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ISSUE:  Accurate  vessel  positioning,  survey 
controls,  and  dredge-travel  monitoring  are  im¬ 
portant  to  many  facets  of  a  dredging  program, 
ranging  from  the  planning  phase  to  contract 
payment  as  well  as  including  satisfaction  of 
environmental  concerns.  Presently,  the  sys¬ 
tems  used  for  horizontal  positioning  or  hydro- 
graphic  survey  vessels  and  dredges  require 
daily  calibration  with  a  known  point  (shore 
station).  The  shore  stations  are  extremely  ex¬ 
pensive  and  labor  intensive  to  calibrate  and 
maintain.  Further,  dredging  and  survey  opera¬ 
tions  are  vertically  referenced  to  the  vessel 
performing  the  work.  This  reference  is  contin¬ 
uously  in  error  by  short-  and  long-term  sea 
surface  action.  Thus  surveyors  cannot  accu¬ 
rately  define  the  datum  at  a  job  site. 

RESEARCH:  A  kinematic  positioning  sys¬ 
tem  was  designed  that  would  support  the 
Corps  dredging  and  hydrographic  survey  oper¬ 
ations.  Because  of  accuracy  contraints  ( 10  cm 
vertical  positioning),  only  survey-quality  GPS 
receivers  were  considered;  use  of  auxiliary 
.sensors  was  studied  to  improve  navigation  per¬ 
formance  when  GPS  was  not  available.  The 
analysis  included  simulation  of  three  import¬ 
ant  areas  of  .system  performance; 


•  Ability  to  resolve  GPS  carrier-lane 
ambiguities. 

•  Ability  to  recover  from  loss  of  GPS  signal. 

•  Navigation  .system  accuracy  operating  with 
one  (or  more)  differential  stations. 

Based  on  the  analysis  and  examination  of  rele¬ 
vant  costs  and  risk  factors,  a  decision  matrix 
was  prepared. 

SUMMARY:  From  the  decision  matrix,  sys¬ 
tems  emerged  as  strong  candidates  as  a  GPS 
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for  maintenance  of  navigational  accuracy  dur¬ 
ing  complete  GPS  outages.  The  system  analy¬ 
sis  and  decision  matrix  will  be  used  to  de¬ 
velop  a  prototype  GPS  that  will  be  field  tested. 
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1 . 0  Summary 


This  report  documents  the  results  of  a  systems  analysis  of  a  Global 
Positioning  System  (GPS)  based  navigation  system  to  support  U.S.  Army  Corps  of 
Engineers  (USAGE)  dredging  and  hydrographic  surveying.  A  horizontal 
positioning  accuracy  of  2  meters  is  required;  however,  decimeter  level 
accuracy  is  highly  desired.  In  the  vertical,  decimeter  level  accuracy  is 
required.  Due  to  these  stringent  accuracy  constraints,  a  significant  portion 
of  this  effort  was  dedicated  to  an  accuracy  prediction  based  upon  extensive 
error  modeling  and  simulation.  These  efforts  are  summarized  in  the  next 
section.  Based  upon  the  results  of  the  analysis,  and  examination  of  relevant 
cost  and  risk  factors,  a  decision  matrix  was  prepared  (see  Table  9.0-1). 

The  analysis  included  simulation  of  three  important  areas  of  system 
performance ; 

•  the  ability  to  resolve  GPS  carrier  lane  ambiguities 

•  the  ability  to  recover  from  loss  of  GPS  signals 

•  navigation  system  accuracy  operating  with  one  (or  more)  differential 

reference  stations. 

Decimeter  level  positioning,  using  GPS,  requires  use  of  carrier  phase 
measurements  which  have  the  integer  cycles,  or  lanes,  properly  resolved. 

The  first  simulation  investigated  the  ability  and  time  required  to 
resolve  the  lanes,  considering  the  possible  receiver  types  (i.e.,  P-code, 
single  frequency  C/A-code,  and  dual  frequency  C/A-code  types).  The  results 
show  that  the  so-called  widelanes  can  be  resolved  within  7  minutes  using  P- 
code  receivers  at  distances  up  to  100  km.  Widelanes  could  be  resolved  only 
for  limited  separation  distances  using  dual  frequency  "codeless"  type 
receivers  (which  derive  second  frequency  information  through  "squaring"  or 
"cross  -  correlating" ,  as  explained  later  in  the  report).  Section  4  describes 
the  simulation  and  results  in  detail. 

In  the  second  simulation,  the  use  of  an  Inertial  Navigation  System  (INS), 
high  quality  clocks,  and  a  barometric  altimeter  were  examined  for  their  use  in 
maintaining  integers  during  gaps  in  GPS  coverage.  Different  quality  INSs  were 
considered  with  both  partial  and  complete  obscuration  of  GPS.  Maximum 
tolerable  gap  lengths  (without  loss  of  the  integer  ambiguities)  were  deter¬ 
mined  as  a  function  of  the  separation  distance  from  the  differential  reference 
station.  The  maximum  tolerable  gap  for  total  loss  of  GPS  and  a  high  quality 
INS  was  found  to  be  about  a  minute.  A  detailed  reporting  of  these  efforts 
appears  in  Section  5. 
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The  third  simulation  predicted  the  resultant  navigation  performance  wiser, 
operating  with  both  single  and  multiple  reference  stations  operating  as  a 
network.  Separation  distances  of  up  to  120  km  are  considered.  .Hulciple 
reference  stations  were  necessary  to  maintain  decimeter  vertical  accuracy  at 
longer  separation  distances.  The  results  are  reported  in  detail  in  Section  6. 
The  conclusions  of  the  study  are  as  follows; 

•  A  decimeter  level  kinematic  positioning  system  employing  differential  GPS 

is  feasible. 

t  The  maximum  operating  range  requires  use  of  P-code  receivers.  Lesser 

receivers  provide  only  limited  range.  A  suitable  receiver  is  presently 
on  the  market,  though  it  is  expensive. 

•  A  single  reference  station  can  supply  decimeter- level  accuracy  out  to 

about  20  km  once  single  frequency  (LI)  lanes  have  been  determined. 

•  Multiple  reference  stations  are  needed  to  maintain  decimeter  vertical 

accuracy  distances  greater  than  20  km.  The  user  must  be  surrounded 
by  the  network  stations  (interpolation) . 

•  Wide -lanes  can  be  resolved  and  the  system  initialized  within  ss?  minutes 

using  six  available  satellites  from  the  full  GPS  constellation.  Use  of 
GLONASS  and/or  geostationary  satellites  can  speed  up  the  process. 

•  Integration  of  a  high  quality  INS  allows  complete  loss  of  GPS  signals  for 

periods  up  to  a  minute  without  having  to  reacquire  the  lanes. 
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2.0  Introduction 


2.1  Review  of  Current  Dredging  Operations 

A  major  function  of  USAGE  is  to  maintain  the  Nation's  ports,  harbors  and 
waterways.  Extensive  hydrographic  surveys  are  required  to  support  the 
planning,  engineering,  design,  and  construction  phases  of  this  effort. 
Hydrographic  surveying  and  vessel  positioning  is  an  essential  engineering 
function  to  prepare  for  and  administer  dredging  projects.  Presently, 
horizontal  positions  of  dredges  and  hydrographic  survey  vessels  a  a  determined 
using  systems  that  electronically  measure  multiple  ranges  or  ranges  and  angles 
from  previously  established  control  points  on  shore.  Most  of  these  systems 
require  the  vessel  to  occupy  a  calibration  point  installed  near  the  job  site 
each  work  day  to  "initialize"  the  system.  Furtnermore,  they  all  require 
establishing  a  series  of  receiver/transmitters  on  control  stations  on  shore. 
Maintaining  these  control  stations,  moving  receiver/transmitters  about,  and 
performing  the  calibration  process  is  extremely  expensive  and  labor  intensive. 
In  addition,  all  dredging  and  survey  operations  are  vertically  referenced  to 
the  vessel  performing  the  work.  This  reference  is  generally  related  to  a  tide 
or  river  gage  to  reduce  the  depth  readings  to  some  datum,  for  instance  mean 
lower  lov?  water  (MLLW) .  This  method  assumes  that  water  surface  elevations  at 
the  gage  site  accurately  represent  the  surface  elevations  at  the  survey  site. 
However,  this  is  not  the  case  and  the  surface  elevations  vary  significantly 
between  the  gage  and  the  survey  site.  Offshore  tide  gages  have  been 
introduced  as  a  means  to  produce  mathematical  models  of  the  surface 
characteristics  of  a  body  of  water.  However,  these  are  expensive  to  install, 
operate  and  maintain.  Furthermore,  the  models  produced  are  limited  in 
accuracy  by  various  tidal  characteristics,  and  metrological,  oceanographic, 
and  hydrological  effects.  Current  technology  does  not  allow  surveyors  to 
accurately  and  efficiently  define  the  vertical  datum  at  a  job  site.  With  GPS, 
vessel  altitude  is  determined  relative  to  the  WGS-84  ellipsoid.  Once  an 
onshore  tide  gauge  reading  has  been  referenced  to  WGS-S^  ellipsoidal  elevation 
(by  collocating  a  GPS  receiver  with  the  tide  gauge),  vessel  depth  gage 
readings  can  be  determined  relative  to  MLLW  without  the  use  of  water  surface 
modelling . 

A  system  requiring  no  onsite  calibration  and  only  one  reference  station 
would  significantly  increase  the  efficiency  and  productivity  of  the  USAGE 
hydrographic  surveyor.  If  this  system  could  simultaneously  provide  sufficient 
vertical  accuracy  it  could  also  serve  a*:  a  means  to  rectify  the  tidal  datura. 
For  such  a  system  to  be  accepted  by  the  community  it  would  have  to  provide  a 
horizontal  positioning  accuracy  of  ±  2  meters,  and  a  vertical  accuracy  of  ± 

0.1  t..eters.  Even  though  a  GPS  kinematic  positioning  system  may  not  require 
calibration  to  its  reference  station,  it  may  have  to  be  calibrated  or  aligned 
to  the  local  project  datum. 

With  the  deployment  of  GPS,  a  possibility  for  a  kinematic  positioning 
system  meeting  these  requirements  has  materialized.  Relative  positions  of  a 
moving  antenna  meeting  the  above  criteria  are  possible  using  GPS. 

Differential  positioning  using  carrier  smoothed  pseudo- ranges  has  demonstrated 
accuracies  of  a  few  meters.  Differential  kinematic  surveys  using  carrier 
phase  measurements  have  demonstrated  accuracies  o£  a  few  centimeters. 
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However,  €'xpe  r  imetu  s  using  cb.e  later  approacti  emploved  niethods  to  initiali.-i 
carrier  lane  ambiguLf.'  which  are  operationailv  iraippropr  i  ate  in  the  iraritc.- 
envi  ronjiient  (occupation  ot  known  baseline,  antetir.a  swaps,  .static:  siiivevs  , 
etc.').  Acidic  Lon;:tl  operacions  may  be  needed  to  recover  tul  1  pt- r !  onr-.itu'e  a::.,- 
loss  ot  signal  and  cvcie  slips.  Data  t  ransiEi  ss  i  on ,  (landling,  .iiui  r..  dia:  t  ;  i,;; 
iiave  nor  been  fullv  developed  for  operational  use. 

2.2  Review  of  US.ACL  Study  Objectives 

Currently,  tlu're  is  no  kinematic  GPS  positioning  sy.stem  which  na-e's  tiic 
operctional  needs  previously  outlined.  The  purpose  of  this  studv  is  t.u  Ccii 
and  analyze  possible  candidate  system  approaches,  determine  their  pt-r for.'tar.c 
and  estimate  costs  to  implement  and  operate. 

2  2.1  Kinematic  Po.sitioning  System  Requirements 

For  the  systems  under  study  to  be  successful,  they  mu.iC  meet  tiu* 
following  requirements : 


They  should  provide  horizontal  positional  accuracies  equivalent  to  tn.i; 
obtained  with  existing  systems,  approximately  i  2  meters.  Decimeter 
horizontal  accuracy  is  highly  desired. 

Existing  water  surface  modeling  techniques  using  on.site  gages  can  ptovi 
vertical  accuracies  of  about  0.1  meter,  however  these  are  not  repeatable  it 
this  level  due  to  metrological  and  hydrological  conditions,  etc.  Therefore, 
the  proposed  systems  should  be  capable  of  redefining  the  vertical  datum  ar  t 
0.1  meter  level  cr  better. 

USAGE  uses  automated  hydrographic  survey  systems  to  produce  engineering 
drawings  for  the  planning,  design  and  construction  of  major  engineering  work 
These  principally  include  navigation  channels  and  dredged  material  di.spo.sal 
areas.  Most  importantly,  the  results  of  these  survevs  are  used  for  bcutiget 
estimates  and  contractor  payments.  As  a  result,  the  scale  of  these  survov.s 
requires  the  positions  be  computed  at  least  once  per  second. 

It  is  desired  that  positions  be  computed  in  "real  time"  so  that  tia'v 
be  provided  to  other  on-board  systems  when  the  event  actually  occurred. 
Alternately,  a  time  tag,  accurate  to  1  millisecond,  must  be  provided  for 
post-mission  matching  of  soundings  and  positions. 

Survey  vessels  and  dredges  must  be  accuratelv  piloted  along  a  pre¬ 
determined  course  or  to  a  specified  position.  Horizontal  position.s  rrufjt:  :>»• 
provided  in  "real  time"  to  an  accuracy  of  ±  2  meters,  for  mission  rwiviga'ior. 
Height  in  real  time  is  required.  Other  requirements  and  desirable  svsifr. 
features  relative  to  vessel  positioning  are  .summarized  below: 

•  The  system  must  operate  in  the  dynamic  conditions  exnerienced  in  small 
craft  (20  to  65  feet)  under  all  weather  conditic)n.s  in  whicii 
hydrographic  survey  operations  are  performed 
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•  The  system  must  determine  positions  of  other  on-board  systems,  such  as 

depth  sounder  transducers,  which  are  not  coincident  with  the  GPS 

antenna(s) . 

•  The  system  must  minimize  the  need  for  shore  based  control  and/or 

reference  stations.  A  single  reference  station  should  have  a  broad 
coverage  area  and  be  capable  of  servicing  multiple  users 
simultaneously . 

•  It  is  desired  that  no  special  calibration  or  operational  procedures  be 

required  to  initialize  the  system  to  conduct  a  mission.  Any  needed 
pre-mission  operations  must  be  within  the  capability  of  the  vessel  crew 
to  perform. 

•  The  system  must  be  robust  in  operation.  It  should  be  designed  to 

minimize  loss  of  GPS  data.  It  is  desired  that  performance  degrade 

gracefully  after  signal  loss.  Recovery  of  full  performance  should  be 
rapid  when  sufficient  signals  are  reacquired.  It  is  desired  that  nj 
special  operational  procedures  be  required  to  restore  full  performance. 

•  The  system  must  be  reliable  with  little  downtime.  When  troubles  do 

arise,  general  repairs  will  normally  be  made  by  survey  technicians. 

•  The  system  must  be  relatively  small  and  light  weight.  Many  USAGE  survey 

vessels  are  in  the  20  to  25  foot  class  with  limited  space  and  weight 
capacity.  A  reasonable  weight  limit  would  be  120  pounds  or  less. 

•  USAGE  is  anticipating  the  development  of  a  helicopter  mounted  LIDAR 

Bathymeter.  A  kinematic  positioning  system  developed  under  this  work 
unit  should  maintain  a  goal  of  interfacing  to  this  helicopter  system. 
This  system  will  have  very  rigid  mission  requirements  in  regards  to 
dynamics,  size,  and  weight;  although  they  are  yet  to  be  defined. 

•  The  system  should  be  capable  of  installation  on  a  buoy  to  be  used  as  a 

reference  point  for  establishing  off  shore  tidal  datums.  Again,  the 
specific  requirements  for  this  application  have  not  been  defined,  but 
they  would  certainly  be  the  most  rigorous  for  size  and  weight  limits. 
However,  the  data  processing  application  may  be  the  simplest  because  of 
the  limited  motions, 

•  The  system  must  operate  with  Block  II  and  Block  HR  satellites  as  well  as 

with  the  present  constellation.  It  Is  fully  expected  that  Selective 
Availability  (S/A)  and  Anti -Spoof ing  (AS)  will  be  implemented  in  the 
Block  II  constellation.  It  Is  desired  that  the  system  not  require 
classified  access  for  full  performance. 

•  The  system  should  have  lower  cost  and/or  better  performance  than  the 

present  methods  for  hydrographic  surveying. 
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Summary  of  Study  Activities 
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The  study  activities  are  summarized  in  the  following  paragraphs: 

•  Develop  alternate  system  implementations  for  analysis. 

•  Prepare  detailed  descriptions  of  selected  systems,  including  hardware, 

software,  and  interface  requirements,  error  budgets,  reference  station 
location(s) ,  data  storage,  communications,  processing  requirements, 
calibration/initialization,  operational  procedures  and  constraints. 

t  Analyze  and  simulate  system  performance  for  representative  operational 
scenarios.  The  analysis  shall  include:  accuracy;  sensitivity  to 
calibration/  initialization  procedures;  sensitivity  to  signal  loss;  and 
accuracy  degradation  and  recovery  from  signal  loss  and  carrier  cycle 
slips . 

•  For  operationally  viable  system  implementations,  estimate  the  system 

costs  including:  development;  procurement;  operation;  and  data 
reduction.  Identify  sources  for  all  components  available 
off-the-shelf.  Scope  the  development  effort  and  assess  the  risk  for 
items  needing  further  development. 

•  Identify  areas  of  potential  growth  and  barriers  to  growth  (e.g., 

helicopter  operation) . 

•  Evaluate  the  tradeoffs  between  system  performance,  risk,  cost  and  growth 

potential . 

•  Prepare  a  decision  matrix  comparing  the  performance,  risk,  cost,  and 

growth  potential  of  the  alternative  system  implementations. 

2.3  Discussion  of  Technical  Approach 

2.3.1  Identification  of  Critical  Issues 

Based  upon  the  stated  navigation  system  requirements,  four  critical 
technical  issues  affecting  performance  were  identified: 

1.  Rapid  and  reliable  acquisition  of  the  lane  ambiguities  associated 
with  the  carrier-phase-derived  pseudorange  measurements  in  the 
expected  dynamic  environment. 

2.  Maintenance  of  positioning  capability  during  partial  or  complete  GPS 
obscuration  and  reconvergence  to  the  correct  set  of  carrier-phase 
integer  ambiguities  upon  reacquisition  of  GPS  signals. 

3.  Estimation  and  removal  of  the  correlated  errors  (i.e.,  S/A  and 
tropospheric  and  ionospheric  delay)  in  the  GPS  measurements  over  the 
expected  range  of  separation  distances  from  the  reference  stations. 
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4.  Maintenance  of  carrier  lock  and/or  detection  and  correction  of  cvcie 
slips  for  the  expected  signal  dynamics  and  noise  envirorutient  . 

The  ability  of  the  system  to  satisfy  the  stringent  accuracy  requirenu-n’s 
is  critically  dependent  on  establishing  the  correct  set  of  lane  ambiguities, 
since  each  lane  error  contributes  a  minimum  of  roughly  20  cm  of  ranging  error 
(for  the  nominal  LI  wavelength).  These  ambiguities  must  be  determined  in  a 
dynamic,  rather  than  a  static  (i.e.,  survey)  environment.  The  time  to  acquire 
che  correct  set  of  ambiguities  can  adversely  affect  system  operation.s  if  it  i;; 
coo  long.  For  example,  requiring  more  than  20  minutes  for  initial  convergence 
(and  perhaps  reconvergence  if  the  correct  set  is  lost)  can  make  the  system 
unusable  for  some  applications. 

The  ability  to  maintain  integer  count  during  GPS  obscurations  requires 
the  use  of  external  sensors,  including  (generally)  a  high  quality  (e.g.. 
rubidium  standard)  clock  and  an  INS.  The  required  sensor  quality  will  depend 
upon  the  length  of  the  measurement  gap  and  the  extent  of  GPS  coverage  loss 
(i.e.,  how  many  of  the  minimum  required  satellites  are  obscured). 

Residual  correlated  errors  in  the  carrier  pseudorange  measurements  made 
available  for  navigation  limit  the  achievable  positioning  accuracy,  since  they 
cannot  be  filtered  (at  least  in  the  short  term),  and  are  difficult  to  estimate 
reliably.  These  correlated  errors  can  also  adversely  affect  both  the  speed 
and  reliability  of  the  ambiguity  resolution  process.  It  is  therefore  critical 
that  the  ability  of  single  or  multiple  reference  stations  to  remove  these 
error  contributions  be  quantified. 

Loss  of  carrier  lock  is  readily  detected  by  most  GPS  receivers  within  20 
msecs.  This  detection  can  be  performed  reliably  by  each  receiver  channel, 
independent  of  the  other  channels,  without  use  of  any  externally- supplied 
information.  Once  a  (momentary)  loss  of  track  is  detected  on  one  channel,  a 
multichannel,  state-of-the-art  GPS  receiver  [2,4,36]  should  be  able  to 
determine  the  number  of  carrier  cycles  which  were  lost  when  track  has  again 
been  established.  The  redundancy  available  in  the  additional  satellites 
should  permit  reconstruction  of  the  (missing)  carrier  cycles.  In  addition,  if 
separate  tracking  of  both  LI  and  L2  carriers  is  performed  within  the  receiver, 
this  information  can  be  used  to  determine  the  number  of  slipped  cycles  on  one 
carrier,  as  long  as  both  carriers  do  not  lose  lock  over  the  same  time  period. 
Thus,  an  examination  of  recovery  from  carrier  loop  loss  of  lock  will  not  be 
performed  as  part  of  this  study,  since  it  can  and  will  be  done  reliably  within 
the  recommended  GPS  receiver  architectures  [1,8,38],  Current  experience  with 
testing  Trimble  Navigation,  Ltd.  (TNL)  survey  receivers  indicates  that  cycle 
slips  are  generally  induced  by  interfering  sources  or  poor  signal  strength; 
i.e.,  caused  by  external  radio  frequency  sources  or  partial  obscuration  of  the 
GPS  signals.  Independent  of  the  cause,  methods  are  used  which  can  detect 
cycle  slips  very  reliably.  Although  no  comprehensive  statistical  analysis  has 
been  performed  for  TNL  receivers,  probabilities  of  undetected  cycle  slips  of 
less  than  one  tenth  of  a  percent  have  been  observed  in  tests  of  simulated 
cycle  slips  over  extended  time  periods  (i.e.,  many  hours). 
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In  addition  to  the  critical  issues  identified  and  discussed  above, 
certain  other  technical  issues  listed  below,  are  addressed  in  the  report 

•  Signal  reception  in  the  expected  multipath  and  dynamic  environment. 

•  Antenna  design,  location  and  mounting. 

•  System  initialization/calibration  requirements. 

•  Referring  positions  to  and  time  synchronization  with  other  on-board 

systems . 

•  Impact  of  S/A  and  AS, 

2.3.2  Simulation  Requirements 

The  preceding  discussion  of  critical  technical  issues  motivates  three 
separate  but  related  simulation  activities: 

1.  A  simulation  of  the  integer  ambiguity  acquisition  (and  reacquisition) 
process . 

2.  A  simulation  of  the  use  of  external  sensors  in  filling  in  GPS 
obscuration  periods. 

3.  A  simulation  of  the  use  of  single  and  multiple  reference  stations  in 
removing  correlated  error  components  in  the  GPS  measurements. 

These  three  simulations  have  been  used  to  address  the  critical  issues 
mentioned  above,  and  are  reported  on  in  the  following  sections  of  the  report. 
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3.0  System  Description  and  Design  Alternatives 

3.1  System  Overview 

The  major  components  of  the  recommended  navigation  system  are  i  1  lu.sr rated 
in  Figure  3.1-1.  Note  chat  it  is  comprised  of  two  major  subsystems:  a 
ship-based  unit,  with  an  associated  communications  link  to  a  master  reference 
station;  and  up  to  three  reference  stations,  including  the  master  and  (up  to) 
two  monitors.  The  number  of  monitor  stations  required  is  a  function  of  thie 
desired  accuracy  and  geographic  coverage;  over  reasonably  short  distances 
;e.g.,  10  km),  a  single  reference  station  should  meet  USAGE  constraints  on 
accuracy . 

3.1.1  Ship-Based  System 

As  illustrated  in  the  figure,  the  ship-based  system  consists  of  a  GPS 
receiver  with  a  rubidium  frequency  standard  input,  an  INS,  a  barometric 
altimeter,  a  communications  subsystem  and  a  display.  Each  of  these  subsystems 
interfaces  with  the  central  processor,  at  the  heart  of  the  system.  The  GPS 
receiver  is  the  primary  source  of  navigation  information  for  the  ship:  the 
ultimate  accuracy  achievable  is  largely  dependent  upon  the  receiver  quality. 
The  other  navigation  sensors  (including  the  rubidium  clock)  serve  primarily  to 
aid  GPS,  when  fewer  chan  the  required  number  of  satellites  are  received.  INS 
provides  a  full  three-dimensional  backup  navigation  capability,  and  is  in  fact 
calibrated  during  periods  of  full  GPS  coverage,  while  the  barometric  altimeter 
and  rubidium  clock  add  stability  to  the  vertical  and  time  channels, 
respectively,  when  full  coverage  is  unavailable.  Note  that,  even  when 
operating  with  both  monitor  stations,  the  ship-based  system  need  only 
communicate  with  the  master  reference  station.  The  major  functions  to  be 
performed  by  the  central  processor  include: 

•  Providing  the  interface,  in  real  time,  with  the  navigation  sensors,  the 

display  and  the  communications  subsystem. 

•  Implementing  the  Kalman  filter  software  for  integrating  the  INS, 

measurements  from  the  GPS  receiver,  clock  and  barometrically  derived 

altitude . 

•  Implementing  the  differential  correction  software,  based  upon  data 

received  from  the  master  reference  station. 

•  Implementing  the  ambiguity  resolution  algorithms,  for  initial  convergence 

and  needed  reconvergence  to  the  correct  set  of  integer  ambiguities. 

•  Providing  the  interface,  in  real  time,  to  other  on  board  systems,  as 

required. 

The  major  modules  of  the  real-time  kinematic  positioning  software  are 
illustrated  in  Figure  3.1.2.  As  indicated,  the  operator  interface  software 
allows  for  system  control  (i.e.,  starting  and  stopping)  and  display,  ,and 
changing  the  values  of  algorithm  parameters.  The  Moding  and  Sequencing  Logic 
controls  the  execution  of  the  Ambiguity  Resolution  Software  and  the  GPS/INS 


9 


DISPLAY  DISPLAY 


LO 


Figure  3.1-1  Navigation  System  Major  Components 


Integration  Filter  to  ensure  that  a  correct  sec  of  integer  ambiguities  is 
maintained:  initially,  and  each  time  a  new  GPS  satellite  is  acquired,  the 

Ambiguity  Resolution  Software  must  be  executed  to  determine  new  integers; 
given  a  set  of  integers,  the  carrier  pseudorange  measurements  can  be  proct-sseci 
by  the  GPS/INS  Kalman  filter.  This  filter  comprises  17  states,  which 
represent  a  compromise  between  fidelity  (the  truth  model,  as  discussed  in 
Section  5,  includes  56  states)  and  imposed  computational  burden.  The  model l.c 
states  include  3  INS  position,  3  INS  velocity  and  3  INS  attitude  errors,  i 
gyro  drifts,  3  accelerometer  biases  and  2  GPS  receiver  clock  error  scate.s. 
Thus,  the  Kalman  filter  and  the  Ambiguity  Resolution  Software  are  running 
concurrently,  with  the  filter  processing  carrier  pseudorange  measuremenc,s 
derived  from  all  satellites  for  which  lane  ambiguities  have  been  resolved. 

The  resolved  carrier  pseudoranges  are  continuously  monitored  by  the  Kalman 
filter.  The  cycle  slip  information  output  by  the  GPS  receiver  can  be  utilised 
to  resolve  any  momentary  losses  of  track,  and  filter  residual  testing  can  he 
used  to  augment  this  process.  Should  a  cycle  slip  be  detected  by  the  GPS 
receiver,  the  Kalman  filter  state  vector  could  be  augmented  in  an  attempt  to 
solve  for  the  slipped  cycles.  Thus,  "false  alarms”  will  not  degrade 
performance,  since  the  filter's  estimate  for  the  unresolvt  d  lane  ambiguity 
should  not  be  close  to  an  integer  value.  For  cycle  slips  /hich  have  not  been 
detected  by  the  GPS  receiver  (a  very  rare  occurrence  for  most  designs),  Kalman 
filter  residual  (parity)  testing  can  be  used  to  isolate  one  or  more  slipped 
channels  if  sufficient  redundancy  exists  (six  satellite  tracking  permits 
isolation  of  a  single  cycle  slip).  When  integers  have  been  lost  on  a 
satellite  which  is  still  in  view,  (e.g.,  as  induced  by  obscuration)  covariance 
and  position  and  velocity  information  is  passed  to  the  Ambiguity  Resolution 
Software,  to  enable  an  accelerated  search  process. 

Depending  upon  the  outputs  of  INS  (i.e.,  position  and  velocity,  vs.  delta 
velocity  and  delta  angle  increments),  INS  interface  software  changes 
significantly.  For  the  Modular  Azimuth  Positioning  System  (MAPS)  Dynamic 
Reference  Unit  (DRU)  INS,  position  and  velocity  outputs  are  available  at  a  5 
Hz  rate,  which  is  adequate  for  this  application.  However,  there  is  a  timing 
uncertainty  of  80  msecs  associated  with  the  position  output;  since  this  cannot 
be  compensated,  it  can  only  be  modelled  as  an  additional  "noise"  source  by  clie 
Kalman  filter,  with  magnitude  computed  as  a  function  of  the  expected  velocitv 
and  acceleration.  Interpolation  of  INS  outputs  may  be  needed  because  the 
depth  sounder  can  output  6-10  readings/sec.  The  barometric  altimeter 
interface  software  is  optional,  depending  on  the  availability  of  the 
barometric  altimeter,  and  supplies  altitude  information  to  the  filter  at  a  1 
Hz  rate.  The  outputs  of  the  GPS/INS  Kalman  filter  which  represent  the  best 
estimates  of  position,  velocity,  attitude,  attitude  rate  and  time,  are  output 
for  transmission  to  other  shipboard  systems;  the  velocity  and  attitude  race 
information  can  be  used  to  synchronize  the  supplied  position  and  attitude 
information . 

A  secondary  usage  for  INS  is  as  an  attitude  reference,  required  for  lever 
arm  compensation  from  the  GPS  antenna  to  other  points  on  the  ship.  Assuming  a 
20  meter  lever  arm  magnitude,  an  attitude  measurement  accuracy  of  a  fc?w 
milliradians  (mrads)  is  necessary  to  reduce  errors  induced  by  the  lever  arm 
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correction  to  the  centiroecer  level.  Such  an  accuracy  is  readily  achievabU- 
using  a  high  quality  INS  which  is  periodically  calibrated  by  GPS.  If  INS  is 
not  present,  an  alternate  attitude  system  is  required,  as  illustrated  ir; 

Figure  3.1-1.  Such  a  system  could  be  based  upon  the  use  of  GPS  in  a  mu  L  t  i ;;  i  i- 
t.i.e.,  at  least  3)  antenna  system  for  three-axis  attitude  information. 
.Accuracies  approaching  0.5  mrad  are  possible  with  sufficient  antenna 
separation  and  an  environment  relatively  free  from  multipath.  Given  thcit  t  he 
relevant  integer  ambiguities  associated  with  the  carrier  difference 
measurements  can  be  reliably  resolved,  a  differencial  phase  measurement 
accuracy  of  0.5  cm  (which  includes  the  effects  of  both  multipath  ar;d  noi.se  .> 
implies  an  antenna  baseline  requirement  of  5  meters  to  achieve  1.5  mrad.s  ot 
accuracy.  A  relatively  low  cost  rate  gyro  (which  can  be  calibrated  bv  the 
GPS-derived  attitude)  may  be  used  to  augment  the  GPS  determined  attitude  to 
provide  backup  attitude  information  during  obscuration  periods.  Such  a  backup 
system  is  required  only  for  maintenance  of  continuous  attitude,  since  the 
integers  can  be  reacquired  fairly  rapidly  for  reasonable  baseline  lengths 
i.e.g.  ,  acquisition  over  a  1  meter  baseline  can  generally  be  performed  in 
seconds ) , 

The  GPS  Receiver  Interface  Software  supplies  the  raw  pseudorange  and 
carrier  phase  measurements  to  the  Ambiguity  Resolution  and  Kalman  filter 
software,  along  with  the  required  ephemeris  data.  The  Communications 
Subsystem  Interface  Software  provides  the  required  data  received  from  the 
master  reference  station  to  the  Ambiguity  Resolution  Software  and  the 
Differential  Correction  Software,  which  computes  the  differential  corrections 
to  be  applied  to  the  measurements  processed  by  the  Kalman  filter. 

3.1.2  Reference  Stations 

The  master  reference  station  maintains  the  required  high  rate 
communication  with  the  ship,  making  possible  use  of  data  from  two  (monitor) 
reference  stations.  Data  transferred  by  the  master  includes  the  pseudorange 
and  carrier  phase  measurements  on  both  LI  and  L2  for  all  satellites  in  view, 
and  spatial  derivatives  of  the  differential  corrections  derived  from  the 
monitor  station  data.  This  spatial  derivative  information  enables  the  user  to 
adjust  the  differential  corrections  computed  by  the  master  station  to  account 
for  the  spatial  decorrelation  of  the  ionosphere  and  SA.  Thus,  the  master 
station  software  includes  an  algorithm  for  computing  these  derivatives  from 
the  raw  measured  data.  Of  course,  when  operating  without  the  monitor 
stations,  the  partial  derivative  information  is  absent,  and  so  this  software 
is  unnecessary. 

Transmissions  from  the  monitor  stations  to  the  master  occur  relatively 
infrequently  (i.e.,  once  per  minute)  and  include  the  pseudorange  and  carrier 
phase  measurements  on  both  LI  and  L2  for  all  satellites  in  view.  No  special 
computations  are  performed  by  the  monitor  stations. 

3.2  GPS  Receiver  Designs 

Each  GPS  receiver  considered  for  use  in  the  recommended  navigation  system 
must  support  the  following  minimum  set  of  requirements: 
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•  Track  the  worst  case  expected  signal  dynamics  on  the  Ll  C/A  code  and 

carrier  with  minimum  probability  of  loss  of  track  for  the  expected 
minimum  signal  to  noise  ratio. 

•  For  dual  frequency  designs,  track  the  L2  signal  with  minimum  probability 

of  loss  of  track  under  the  maximum  expected  dynamics  and  worst  case 
signal  to  noise  ratio. 

•  Provide  for  output  of  measured  pseudorange  and  carrier  phase  on  Ll  (and 

L2 ,  for  dual  frequency  designs)  at  1  Hz  for  all  satellites  in  view;  the 
pseudorange  and  carrier  phase  measurement  noise  variances  should  not 
exceed  1  m^  and  1  cra^,  respectively,  over  the  full  range  of  signal 
dynamics  and  signal  to  noise  ratios. 

•  Allow  for  use  of  an  external  high  quality  (i.e.,  rubidium  frequency 

standard)  clock. 

Given  these  minimum  set  of  requirements,  four  generic  receiver  designs 
were  considered  in  this  study; 

•  A  P-code  dual  frequency  receiver. 

•  A  C/A-code  single  frequency  receiver. 

•  A  C/A-code  dual  frequency  receiver  with  codeless  extraction  of  L2  by 

squaring. 

•  A  C/A-code  dual  frequency  receiver  with  codeless  extraction  of  L2  by 

cross  correlation  with  Ll. 

These  two  types  of  codeless  L2  tracking  receivers  deserve  further 
discussion.  Since  neither  design  derives  L2  information  through  use  of  the 
P-code,  they  both  incur  a  loss  in  noise  performance  relative  to  a  P-code 
receiver.  In  addition,  both  designs  require  L2  tracking  loop  aiding  by  the  Li 
tracking  loop  to  enable  discrimination  of  the  separate  satellite  signals 
around  the  L2  frequency.  The  squaring  process  removes  the  code  modulation  on 
each  satellite  carrier;  in  so  doing,  it  also  removes  the  ability  to  isolate 
specific  satellites  through  code  correlation.  Similarly,  the  receiver  which 
cross  correlates  band-pass-filtered  versions  of  the  received  signal  about  Ll 
and  L2  cannot  uniquely  identify  each  satellite's  signal  on  L2 .  The  major 
performance  differentiating  feature  of  the  two  designs  relates  to  the  ability 
to  perform  wide  laning.  The  receiver  which  relies  upon  signal  squaring 
effectively  doubles  the  incoming  signal  frequency,  so  can,  at  best,  achieve 
one-half  of  the  wide  lane  by  combining  Ll  and  L2  carrier  phase  information. 

In  addition,  optional  operation  with  GLONASS  and  a  geosynchronous 
satellite  was  considered  to  increase  the  available  number  of  signals  for 
reducing  the  convergence  time  of  the  ambiguity  resolution  algorithms  under 
study . 
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3.3  Use  of  Auxiliary  Sensors 

As  illustrated  in  Figure  3.1-1,  three  auxiliary  sensors  are  considered 
for  dynamic  aiding  during  GPS  outages;  an  INS.  a  barometric  altimeter,  and  a 
rubidium  frequency  standard.  Three  generic  INSs  were  considered  in  this 
study,  representing  high,  medium,  and  low  quality  designs.  The  basic  concept 
in  using  INS  aiding  is  chat,  once  the  INS  has  been  calibrated  during  full  GPS 
coverage,  it  should  enable  "coasting"  during  outages,  i.e.,  maintenance  of  the 
correct  integer  counts  for  a  period  of  time  so  that  an  integer  reacquisition 
may  not  be  necessary  when  the  signals  return.  Use  of  either  the  barometric 
altimeter  or  the  rubidium  frequency  standard  bounds  the  error  growth  during 
GPS  outages  in  a  single  direction;  for  the  barometric  altimeter,  INS  altitude 
error  is  limited  Co  the  residual  bias  error  in  the  instrument  (which  has  been 
calibrated  by  GPS  prior  to  the  outage) ,  while  the  rubidium  clock  limits  the 
temporal  error  growth.  The  rubidium  frequency  is  inherently  very  stable  and 
its  frequency  drift  has  been  estimated  by  the  Kalman  filter  during  periods  of 
full  GPS  coverage. 

As  an  alternative  to  the  use  of  the  barometric  altimeter,  zero  mean 
vertical  velocity  damping  was  considered.  This  approach  relies  on  the  fact 
that,  over  many  cycles  of  ocean  wave  induced  oscillation,  the  average  vertical 
velocity  of  the  ship  should  be  close  to  zero;  thus,  the  average  of  INS 
vertical  velocity  over  an  equally  long  period  should  be  an  accurate  measure  of 
INS  vertical  velocity  error  and  can  be  used  to  bound  the  INS  error  growth  in 
the  vertical  channel. 

3.4  Gomraunications  Subsystems 

As  previously  noted,  there  are  (at  most)  three  separate  communication 
links  which  need  to  be  implemented:  as  a  minimum,  data  must  be  transferred 
from  the  master  reference  station  to  the  ship  at  a  fairly  high  rate  (i.e.,  1 
Hz);  optionally,  the  two  monitor  stations  communicate  with  the  master  at  a 
relatively  low  frequency  (e.g.,  once  per  minute).  The  data  requirements  for 
the  master  reference  communications  are  presented  in  Table  3.4-1.  As 
indicated,  approximately  1.0  kbaud  is  needed  to  transmit  the  data  from  up  to 
10  satellites  at  a  1  Hz  rate.  Note  that  several  special  measures  have  been 
taken  to  reduce  the  required  baud  rate: 

•  The  ranges  for  the  LI  code  and  carrier  phase  and  time  have  been  reduced 

below  their  maximum  possible  values;  it  is  therefore  necessary  to 
detect  "rollovers"  of  these  quantities  in  the  ship  computer  and  make 
necessary  adjustments  before  these  quantities  can  be  used. 

•  Only  the  differences  between  the  LI  and  L2  code  and  carrier  phases, 

rather  than  the  whole  valued  L2  phases  are  transmitted;  the  difference 
variable  must  therefore  only  accommodate  twice  the  worst  case 
ionospheric  delay,  multipath  and  noise. 

Data  requirements  for  the  monitor  station  communications  are  indicated  by 
Table  3.4-2,  Similar  measures  have  been  taken  to  reduce  the  baud  rate  of  the 
link,  resulting  in  a  requirement  of  0.02  kbaud. 
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Table  3.4-2 
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Assume  {v,orst  case)  10  satellites  1530  btts/60  sec; 


Mechanizac ion  of  the  coimiiunicaclons  can.  b»-  accoiap I  i  •ha.-a  m,  in.: 

modems  and  standard  Leleplione  lines  for  data  transfer  bt^twet-n  !  ra-  iruini't.r 
stacions  and  th(?  master  reference  station,  .and  a  V.-rv  Hii'.h  Kreqe.i-ne  v  ’.'i;;- 
Ultra  HigVi  Frequency  (1,'HF)  line  of  sight  radio  litik  Ijetveen  the  it.istt-r 
and  the  ship.  Lhse  of  telephone  lines  for  the  monitor  stations  is  -to :  i  •.•.ar  .-d  t  ■, 
the  relatively  low  baud  rate;  of  course,  a  rer'.tal  fee  would  he  as.soc  iati-d  wit,'. 
tf;e  use  of  the  telephone  lines.  The  line  of  sight;  radio  linV;  bt-twci-n  the 
master  reference  station  and  the  ship  can  only  accosmnodate  the  largest 
e.tpected  separation  distaitce  (100  kin)  if  a  sufficiently  Large  r.idlo  tow.-!  it, 
utilized  and  power  is  sufficiently  large  Unless  such  a  tower  i  .s  .■.ire.idv 
available,  its  construction  could  substantially  ii-.crease  the  cost  of  the 
co.TJiiunication  subsystem  required. 
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4.0  Ambiguity  Resolution  Simulation 

This  section  describes  the  simulation  of  the  resolution  of  widelane 
carrier  phase  ambiguities  in  differential  kinematic  positioning.  The  abiliiv 
to  resolve  the  LI  lane  ambiguity  is  not  addressed.  This  limitation  of  the 
analysis  will  be  further  discussed  in  Section  6,  when  resultant  navigation 
performance  predictions  are  presented  based  upon  both  widelaning  and  knowledge 
of  the  LI  lane.  Generally  speaking,  acquisition  of  the  LI  lane,  given  that 
the  widelane  has  been  acquired,  is  not  straightforward.  Determination  and 
fixing  of  the  correct  integer  values  of  the  carrier  phase  ambiguities  is  a 
requirement  for  sub-decimeter  kinematic  positioning  with  GPS.  Some  useful 
methods  for  ambiguity  resolution  in  land  kinematic  positioning  have  been 
developed  over  the  past  few  years.  These  methods  require  a  static 
initialization  and  are  therefore  not  applicable  in  an  operational  marine 
environment . 

Three  methods  can  be  used  for  ambiguity  resolution  "on-the-fly"  (OTF) , 
i.e.,  while  the  remote  receiver  is  in  motion.  The  first  method  resolves 
ambiguities  by  combining  pseudorange  and  carrier  phase  measurements;  this  has 
been  proved  successfully  in  the  field  with  dual  frequency  P-code  receivers 
[34],  and  is  chosen  for  further  simulation  in  this  study. 

A  second  method  is  based  solely  on  redundancy  in  integrated  Doppler 
measurements  to  multiple  satellites  [24,29].  It  is  analogous  to  current 
Doppler  survey  techniques,  but  uses  three  extra  satellite  signals  to  track  the 
motion  of  the  vehicle.  This  method  has  never  been  tested  in  the  field; 
simulation  has  demonstrated  very  long  convergence  times  requiring 
uninterrupted  carrier  tracking,  especially  under  the  effect  of  differential 
atmospheric  delay  and  multipath.  Consequently,  this  technique  is  deemed 
impractical  for  real-time  kinematic  positioning  in  the  context  of  dredging. 

The  third  method  of  ambiguity  resolution  was  also  considered  but  rejected 
because  of  the  perceived  risk  associated  with  its  use.  This  method,  the 
ambiguity  function  method,  searches  for  the  most  likely  set(s)  of 
double -difference  integers  within  the  meter-level  "coarse"  bounds  established 
by  double -difference  code.  With  reasonably  good  code  measurements  these 
bounds  can  be  as  small  as  twenty  wavelengths  (integers)  in  diameter.  To 
illustrate  this  method,  examine  the  situation  where  only  four  satellites  are 
available.  Then  three  double -difference  integers  (an  integer  triplet)  are 
necessary  to  resolve  the  ambiguity;  there  are  roughly  6000  triplets  within  the 
sphere  described  by  the  coarse  bounds  and  any  one  of  them  is  equally 
acceptable.  If  a  fifth  satellite  is  added,  a  fourth  double-difference  integer 
that  is  linearly  related  to  the  other  three  is  now  required  for  a  full 
solution.  Theory  indicates  that  very  few  integer  triplets  will  allow  a  fourth 
integer  that  satisfies  the  linear  relationship  within  the  noise  bounds  of  the 
measurement.  This  linear  relationship  can  be  thought  of  as  a  family  of 
planes,  a  wavelength  apart  and  twice  thickness  of  the  carrier  phase  error, 
that  cuts  through  the  lattice  of  possible  points  indicated  by  the  triplets. 
With  a  typical  carrier  noise  of  0.05-0.1  wavelength,  the  percentage  of  volume 
in  the  planes  (i.e.,  satisfying  the  relationship)  is  only  1%  of  the  total 
volume  of  the  space.  Thu;j ,  the  extra  restriction  eliminates  99%  of  the 
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possible  triplets  of  the  first  three  satellite  pairs.  Once  the  full  GPS 
constellation  is  up,  one  can  depend  on  six  satellites  with  short  outages;  this 
produces  two  extra  restrictions,  eliminating  99.99%  of  the  possible  criplet.s 
With  only  6000  possible  triplets,  six  satellites  will  usually  discriminate  the 
true  solution  immediately.  Similar  methods  have  been  used  in  static  survey 
applications  for  many  years  (ambiguity  function  technique,  [5,30]).  Recent 
advances  in  algorithms  and  faster  microprocessors  are  combining  to  make  this 
technique  much  more  practical  than  previously  thought  possible  (ambiguity 
search  technique ,[ 15 , 9 ]) . 

The  difficulty  with  these  ambiguity  function  and  ambiguity  search 
techniques  is  that  in  field  conditions  there  is  a  possibility  of  large  carrier 
measurement  bias,  on  the  order  of  a  centimeter.  These  biases  can  last  a 
minute  (multipath  at  the  reference  station)  or  many  minutes  (atmospheric 
delays  over  twenty  kilometers)  or  virtually  permanent  (epsilon  orbit  error 
over  twenty  kilometers).  A  static  receiver  can  average  through  these  effects, 
but  a  kinematic  receiver  cannot  The  response  of  the  ambiguity  search  methods 
to  such  biases  is  nonlinear;  although  perturbing  the  measurements  very 
slightly  does  not  change  the  integer  solution  at  all,  a  substantial 
perturbation  of  a  centimeter  may  cause  convergence  on  a  set  of  integers  far 
from  the  true  solution.  The  effect  that  incorporation  of  additional 
satellites  (e.g.,  through  the  Inclusion  of  GLONASS)  has  on  rectifying  this 
problem  is  a  subject  of  current  research  activities.  This  differs  from  the 
code/carrier  convergence  technique  chosen  for  study,  which  always  converges  on 
the  correct  wave  integer.  Estimating  the  biases  is  possible;  if  no  a  priori 
information  on  the  size  of  these  biases  is  available,  the  estimation  problem 
degrades  to  the  multiple  satellite  integrated  Doppler  approach  mentioned  above 
which  converges  very  slowly  unless  many  satellites  are  available.  The 
ambiguity  techniques  do  work  well  in  short-baseline  static  surveys  because 
multipath  errors  are  averaged  over  a  substantial  period.  The  ambiguity 
function  and  ambiguity  search  techniques  were  not  considered  as  viable 
real-time  kinematic  techniques  for  the  following  two  reasons:  OTF  performance 
in  field  conditions  with  reference  station  multipath  and  atmospheric  biases 
has  not  been  demonstrated;  and  quality  assurance  methods  have  not  yet  matured 
to  the  point  that  they  can  be  considered  reliable  for  OTF  ambiguity 
resolution.  The  ambiguity  techniques  work  well  in  static  survey  because  the 
data  is  averaged  over  a  substantial  period. 

The  simulation  is  performed  in  two  steps  for  various  scenarios.  The 
first  step  consists  of  generating  GPS  measurements  based  on  a  particular 
satellite  constellation,  a  model  for  the  measurement  errors  and  a  model  for 
the  motion  of  the  remote  receiver.  In  the  second  step  a  position  and 
ambiguity  determination  algorithm  tries  to  recover  the  carrier  phase 
ambiguities  together  with  the  differential  position.  For  the  present  study, 
the  differential  position  is  a  by-product  of  the  algorithm  and  will  not  be 
further  analyzed  or  discussed.  Section  4.1  describes  the  models  used  in 
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generating  the  simulated  measurements.  Section  4.2  describes  the  solution 
algorithm  employed  to  recover  the  carrier  phase  ambiguities.  Finally,  Section 
4.3  summarizes  the  numerical  results  of  the  study. 

4.1  Measurement  Simulation  Program 

The  measurement  simulation  program  generates  a  file  containing  satellite 
coordinates,  carrier  phase  measurements  and  pseudorange  measurements  for  all 
visible  satellites  above  15  degrees  of  elevation  at  a  one  second  data  rate. 

The  simulated  measurements  are  based  on  a  set  of  satellite  ephemerides ,  a 
model  for  the  receiver  motion  and  models  for  measurement  errors.  Ephemerides 
reference  epoch  is  00^  00"'  00®  of  day  364,  1990. 

4.1.1  Satellite  Ephemerides 

The  ephemeris  files  have  identical  format  for  all  satellites  used  in  this 
study.  The  standard  University  of  New  Brunswick  GPS  ephemeris  file  format  has 
been  used,  containing  seven  records  per  satellite  with  the  GPS  broadcast 
ephemerides  parameters  organized  according  to  Table  4.1.1. 

For  Che  purpose  of  this  study  only  nominal  orbital  parameters  are 
required.  These  nominal  orbital  parameters  describe  circular  orbits  with 
constant  rate  of  right  ascension,  and  the  resulting  ephemeris  files  contain 
many  zeros.  However,  the  general  file  format  has  been  retained  to  enable  the 
processing  of  real  GPS  broadcast  ephemeris  data,  if  required.  The  ephemeris 
files  used  in  this  study  typically  will  have  the  structure  as  indicated  in 
Table  4.1-2  (xx.xx  means  non-zero  value). 

4. 1.1.1  GPS  Ephemerides 

The  full  GPS  constellation  used  in  this  study  is  described  in  [4-3).  The 
24  GPS  satellites  are  irregularly  spaced  in  six  orbital  planes  with  55° 
inclination,  separated  by  60°  in  longitude  of  the  ascending  node.  The  orbits 
are  nominally  circular,  with  a  radius  of  26,609  km,  which  leads  to  an  orbital 
period  of  1/2  of  a  sidereal  day.  The  eccentricity  and  the  argument  of  the 
perigee  are  set  to  zero.  Thus  the  ephemerides  are  given  in  terms  of  Keplerian 
elements  semi-major  axis  a,  inclination  i,  right  ascension  of  the  ascending 
node  Q  and  mean  anomaly  M  for  the  epoch  00:00:00  of  November  26,  1989,  which 
is  day  number  330  of  1989. 

For  use  in  this  study,  the  ephemerides  must  be  extrapolated  to  the 
simulation  epoch  on  day  number  364  of  the  year  1990.  For  this  ephemerides 
prediction,  the  semi -major  axis  and  the  inclination  are  assumed  to  be 
constant.  The  mean  anomaly  of  the  nominal  orbits  repeats  every  12  (solar) 
hours.  Therefore,  at  the  epoch  00:00:00  of  day  364  of  1990,  the  mean 
anomalies  of  the  nominal  GPS  constellation  are  identical  to  the  ones  listed  in 
[4-5].  Right  ascension  of  the  ascending  node  is  subject  to  a  linear  drift  in 
first  order  approximation.  Since  it  is  proportional  to  the  cosine  of  the 
satellite  orbit  inclination: 

dn/dt  -  constant  •  cos  i,  (4.1) 
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Table  4.1-1:  Ephemerls  File  Format 


record  wl: 

14  -  S.V.  number 

N/A 

14  -  GPS  week  Weeks  from  6  Jan  1990 

14  -  User  range  accuracy 

Meters 

14  -  S.V.  health 

N/A 

14  -  Ephemerls  fit  interval 

0-4  hrs , 

1  “  6hrs 

F10.2  -  Age  of  data  ephemerls 

Seconds 

F10.2  -  Ephemerls  reference  time 

Seconds 

F22.15  -  aO  coefficient  for  S.V.  clock 

Seconds 

record  ff2 ; 

F22.15  -  al  coefficient  for  S.V.  clock 

Sec . /Sec . 

F22.15  -  a2  coefficient  for  S.V.  clock 

Sec . /Sec . **2 

F22.15  -  MO  (mean  anolaly  at  reference  time) 

Rad. 

record  #3 : 

F22.15  -  Del  N  (correction  to  mean  motion) 

Rad. /Sec. 

F22.15  -  Eccentricity 

N/A 

F22.15  -  A** . 5(sq . root  of  semimajor  axis) 

Meters  **.5 

record  #4: 

F22.15  -  w  (argument  of  prerigee) 

Rad. 

F22.15  -  10  (inclination  angle  at  ref.  time) 

Rad. 

F22.15  -  I  dot  (rate  of  inclination  angle) 

Rad. /Sec . 

record  #5 : 

F22.15  -  Omega  (right  ascension  of  ref.  time) 

Rad. 

F22.15  -  Omega  dot  (rate  aof  Omega) 

Rad . /Sec . 

F22.15  -  Cue  (cos  harm,  term  to  arg  of  lat) 

Rad. 

record  #6 : 

F22.15  -  Cus  (sin  harm,  term  to  arg  of  lat) 

Rad 

F22.15  -  Cre  (cos  harm,  term  to  orbit  rad) 

Meters 

F22.15  -  Crs  (sin  harm,  term  to  orbit  rad) 

Meters 

record  #7 : 

F22.15  -  Cic  (cos  harm,  term  to  ang  of  incl) 

Rad. 

F22.15  -  Cis  (sin  harm,  term  to  ang  of  incl) 

Rad. 

F22.15  -  Group  delay  differential 

Sec . 

Table  4.1-2:  Ephemerides  File  Structure  for  Nominal  Orbits 

xxxxx  000  0.00  xxxxxx.xx  0 . OOOOOOOOOOOOOOOD-00 

0 . OOOOOOOOOOOOOOOD- 00  0 . OOOOOOOOOOOOOOOD-00  X . xxxxxxxxxxxxxxxD-xx 
0 . OOOOOOOOOOOOOOOD -00  0 . OOOOOOOOOOOOOOOD- 00  X . xxxxxxxxxxxxxxxD-xx 
0 . OOOOOOOOOOOOOOOD -00  X , xxxxxxxxxxxxxxxD-xx  0 . OOOOOOOOOOOOOOOD- 00 
X . xxxxxxxxxxxxxxxD -  XX  X . xxxxxxxxxxxxxxxD - XX  0 . OOOOOOOOOOOOOOOD - 00 
0 . OOOOOOOOOOOOOOOD-00  0 . OOOOOOOOOOOOOOOD-00  0 . OOOOOOOOOOOOOOOD-OO 
0 . OOOOOOOOOOOOOOOD-OO  0 . OOOOOOOOOOOOOOOD-OO  0 . OOOOOOOOOOOOOOOD-OO 


the  rate  of  right  ascension  for  the  GPS  satellites  in  the  nominal  orbits  with 
55°  inclination  can  be  computed  from  the  broadcast  rate  of  right  ascension  of 
an  existing  GPS  satellite  by  using  the  equation: 

dO/dt  jj.  -  4>(,cos  55“  cos  i)  •  dfl/dt  -  (4.2) 

The  resulting  rate  of  right  ascension  for  the  nominal  constellation  amounts  to 
-7.914  •  10"’  rad/sec,  yielding  a  change  in  right  ascension  of  -15.5926“  for 
day  364  of  1990  according  to: 

n  (day  364,  1990)  -  U  (day  330,  1989)  -  15.5926“.  (4,3) 

The  value  Cl^  transmitted  in  the  GPS  broadcast  message  is  not  the  right 
ascension  in  the  sense  of  a  Keplerian  element.  As  pointed  out  in  [4-5], 
is  referenced  to  the  zero  meridian  at  the  beginning  of  the  GPS  week. 

The  relation  between  fl  and  0^  is  given  by: 

n„  -  n  -  GAST„^^,  (4.4) 

where  GASTy^^l^  is  the  Greenwich  Apparent  Sidereal  Time  at  the  beginning  of  the 
GPS  week.  For  the  purpose  of  this  simulation,  CAST  can  be  replaced  by  GMST 
(Greenwich  Mean  Sidereal  Time)  which  is  given  by  the  Astronomical  Almanac 
(1990)  in  units  of  hours: 

GMST  -  6.6265313  +  0.065709824  d  +  1.00273791  t  (4,5) 

where  d  is  the  day  of  the  year  1990,  and  t  is  the  time  of  the  day  in  hours. 

For  the  epoch  00:00:00  of  day  364  of  1990,  this  equation  yields  in  units  of 
degrees : 

GMST  -  98.1736“.  (4,6) 

With  equations  (4.3),  (4.4),  and  (4,6)  we  obtain  finally  the  broadcast  fl^ 
for  the  simulation  epoch: 

»  n  -  15.5926“  -  98.1736*  (4.7) 

The  resulting  broadcast  epheraerides  for  the  24  GPS  satellites  are  listed  in 
Appendix  1.  The  pseudo  random  noise  (PRN)  numbers  assigned  to  the  GPS 
satellites  are  1  through  24. 

4. 1.1, 2  GLONASS  Epheraerides 

The  full  GLONASS  constellation  used  in  this  study  is  described  by  P.Daly 
[6].  The  24  GPS  satellites  are  regularly  spaced  in  three  orbital  planes  with 
64.8“  inclination,  separated  by  120“  in  longitude  of  the  ascending  node.  The 
orbits  are  nominally  circular  at  an  altitude  of  25,507.  km  which  leads  to  an 
orbital  period  of  675.73  minutes.  The  epheraerides  are  given  in  terms  of 
Keplerian  elements  semi-major  axis  a,  inclination  i,  right  ascension  of  the 

ascending  node  Q  and  mean  anomaly  M  for  the  epoch  00:00:00  of  September  7, 

1987,  which  is  day  niimber  250  of  1987. 
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As  for  the  GPS  constellation,  the  ephemer ides  must  be  extrapolated  to  the 
simulation  epoch  on  day  number  364  of  the  year  1990.  The  semi-major  axis  and 
the  inclination  are  assumed  to  be  constant.  The  mean  anomaly  of  the  nominal 
orbits  is  extrapolated  using  the  mean  motion  defined  by: 

dM/dt  -  27r  /  675.73  minutes  «■  0.1549728  •  10'^  rad/sec.  (4.8) 

This  mean  motion  yields  a  phase  shift  of  196.06“  for  the  mean  anomaly  between 
the  epheraerides  epoch  in  [6]  and  the  ephemerides  epoch  for  this  simulation 
s  cudy , 

M  (00:00:00,  day  364,  1990)  -  M  (00:00:00,  day  250,  1987)  +  196.06“  (4.9) 

The  change  in  right  ascension  of  the  ascending  node  is  computed  from  equation 
(42)  for  the  nominal  GLONASS  constellation  and  amounts  to  -6.102  •  10'^ 
rad/sec .  yielding  a  change  in  right  ascension  of  -35.1906°  for  day  364  of  1990 
according  to: 

n  (day  364,  1990)  =  (day  250,  1987)  -  35.1906°.  (4,10) 

■As  discussed  in  the  previous  section,  CAST  must  be  subtracted  from  the 
right  ascension  to  obtain  the  value  fig  transmitted  in  the  GPS  broadcast 
message.  With  equations  (4.3),  (4.4)  and  (4.10)  we  obtain  the  broadcast  fig 
for  the  simulation  epoch: 

-  fl  -  35.1906°  -  98.1736°  (4.11) 

The  resulting  broadcast  ephemerides  file  for  the  24  GLONASS  satellites  is 
listed  in  Appendix  2.  The  PRN  numbers  assigned  to  the  GLONASS  satellites  are 
25  through  48. 

4. 1.1. 3  Geostationary  Satellite  Ephemeris 

By  definition,  a  geostationary  satellite  is  in  an  equatorial  orbit  (i.e. 
inclination  -  0)  with  a  radius  of  about  42,000  km.  It’s  elevation  and  azimuth 
does  not  change  in  time  for  a  fixed  receiver  location.  Since  the 
geostationary  satellite  to  be  used  in  this  study  is  of  fictional  nature  only, 
we  are  free  to  select  its  position  to  ensure  for  the  areas  of  interest  a  high 
satellite  elevation  angle. 

The  satellite  elevation  depends  primarily  on  the  spherical  distance 
between  the  receiver  position  and  the  satellite  footpoint  at  the  equator  as 
shown  in  Figure  4.1.  In  general,  if  the  spherical  distance,  S,  increases, 
the  elevation  angle,  E,  will  decrease.  In  order  to  maximize  the  elevation 
angle  of  the  geostationary  satellite  as  seen  from  different  areas  of  interest, 
Che  corresponding  spherical  distances  must  be  minimized. 

The  relation  between  spherical  distances  and  spherical  coordinates  ,  A 
is  shown  in  Figure  4,2.  To  yield  maximum  satellite  elevation  at  the  U . S .  East 
Coast  (Washington,  D.C.,  039°,  A-77°)  and  in  Alaska  (Anchorage,  061°,  A-150°), 
the  fictional  geostationary  satellite  is  placed  at  a  longitude  A--130°.  This 
satellite  position  results  in  spherical  distances  of  62°  for  both  locations, 
yielding  a  maximum  satellite  elevation  angle  of  19°. 
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Satellite 


Figure  4  1 ;  Satellite  elevation  as  function  of  spherical  distance  between  receiver  and  satellite 


Figure  4,2:  Spherical  distance  between  receiver  and  Geostationary  satellite 
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The  ephemerides  file  (broadcast  GPS  ephemerides  format)  for  the  nomiisai 
orbit  of  a  geostationary  satellite  is  of  particular  simple  structure.  The 
inclination  is  zero  by  definition,  and  the  mean  anomaly  cannot  be  seperated 
from  the  right  ascension.  The  only  remaining  non-zero  elements  are  the  right 
ascension  of  the  ascending  node  and  the  semi-major  axis  of  the  orbit.  Since 
the  broadcast  fl^  is  the  right  ascension  referred  to  the  zero  meridian  at  the 
beginning  of  the  GPS  week,  it  is  identical  to  the  longitude  of  the  a.scending 
node  at  the  beginning  of  the  GPS  week.  This  longitude  has  been  selected  above 
to  maximize  the  satellite  elevation  angle.  Thus  we  have; 

Qq  -  -130.00°.  (‘i.12) 

The  serai-major  axis,  a,  of  a  Keplerian  satellite  orbit  can  be  computed  from 
the  orbital  period,  T,  and  the  gravitational  constant,  using  Kepler's  third 
law : 


a^  -  T^  .  /  (4  3r2)  .  {4.13) 

With  /i  -  3.986008  •  lO'^*  mVsec^  and  T  -  86400  sec,  we  obtain  from  equation 
(4.13)  for  the  semi-major  axis  of  a  geostationary  satellite: 

a  -  42.241  km.  (4.14) 

The  resulting  broadcast  ephemerides  file  for  the  geostationary  satellite  is 
listed  in  Appendix  3.  The  PRN  number  assigned  to  the  geostationary  satellite 
is  25 . 

4.1.2  Receiver  Position  Consideration 

The  coordinates  of  Che  monitor  station  and  the  moving  vessel  need  to  be 
known  for  the  generation  of  the  simulated  GPS  measurements.  Whereas  the 
monitor  station  coordinates  are  time  invariant,  a  new  position  has  to  be 
supplied  for  the  moving  receiver  for  every  measurement  epoch. 

4. 1.2.1  Monitor  Station 

The  ellipsoidal  (geodetic)  coordinates  longitude,  latitude  and  height  of 
the  monitor  station  are  considered  to  be  known  with  sufficient  accuracy.  No 
attempt  is  made  to  include  models  for  differential  positioning  errors 
resulting  from  inaccurate  monitor  station  coordinates. 

4. 1.2. 2  Remote  Station 

For  the  generation  of  simulated  measurements,  a  time  series  of  positions 
of  the  moving  receiver  has  to  be  supplied.  For  the  present  simulation,  this 
position  time  series  should  adequately  represent  the  motion  of  a  dredging 
vessel.  The  following  equations  have  been  used  to  generate  horizontal 
coordinates  latitude  <j>  and  longitude  A  for  the  moving  receiver; 

+  A^(tj)  /  R  (4.15) 

X(t-)  -  Ap  +  AA(tj)  /  (R  cos"V^(Cj)).  (4.16) 
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The  initial  position  (i^g,  Ag)  in  relation  to  the  coordinates  of  the  monitor 
station  determines  the  separation  between  monitor  station  and  user  receiver 
for  a  particular  simulation  run.  The  changes  in  horizontal  coordinates  (units 
of  meters)  are  computed  from; 


A^(t.) 

“  A(^(tj.T) 

+ 

< 

ft 

• 

cos  A(tj) 

(4.17; 

AA(tj) 

-  AA(ti..,) 

+  V(tj)* 

sin  A( tj ) 

(4.18) 

in  terms  of  a  continuously  changing  vessel  speed,  v(tj),  and  a  continuously 
changing  vessel  heading,  A(t.).  Both,  speed  and  azimuth  are  based  on  an 
average  value  and  a  series  of  harmonic  variations  according  to; 


V(tj) 

=  Vg  +  i 

5  V.  .  cos  (v.Atj  +  Mj) 

(4.19) 

j-1 

A(tj) 

1 

> 

o 

+ 

5  A.  •  cos(qjAtj  4-  (Sj) 

(4.20) 

j-1 

with;  Atj  -  tj  -  tg.  The  mean  values  chosen  are; 

Vg  -  1  m/sec  (4.21) 

Ag  -  45’  (4.22) 

and  the  amplitudes,  frequencies,  and  initial  phases  of  the  harmonic  variations 
are  given  by; 


V.  -  (0.055,  0.050,  0.045,  0.040,  0.035,  0.030)  m/sec  (4.23) 

Vj.  -  (27r/10,  2x/13,  2jr/17,  2?r/23,  27r/31,  2;r/36)  (4.24) 

li.  -  (20.07r,  16.97r,  14. Ir,  II.Itt,  9.07r,  8 . 37r)  (4.25) 

A  -  (8°,  r,  20°)  (4.26) 

■7  -  (23r/19,  27r/37,  2t/330)  (4.27) 

5  -  (0.0,  16.27r,  1.2t).  (4.28) 


The  horizontal  components  of  the  vessel  trajectory  contain  constituents 
of  frequencies  between  0.1  Hz  and  0.0003  Hz.  The  deviations  of  the  resulting 
latitude  and  longitude  coordinates  from  a  mean  straight  line  motion  (45 
degrees  of  azimuth)  is  shown  in  Figure  4.3.  The  detailed  horizontal  position 
variations  in  the  first  two  minutes  of  the  simulation  interval  (0  <  <  120 

sec)  are  shown  in  Figure  4.4. 


The  vertical  motion  of  the  vessel  is  considered  to  consist  of  an  initial 
height,  hg,  and  a  time  varying  part,  Dh(tj),  according  to; 

h(tj)  -  hg  +  Ah(tj)  .  (4.2'  ) 

The  time  varying  part  consists  of  a  linear  trend  and  a  series  of  harmonic 
constituents  to  simulate  tidal  variations  (linear  over  half  an  hour)  and  more 
or  less  periodic  height  variations  due  to  swell  and  waves.  It  is  computed 
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Horizontal  coordinate  deviation  of  simulated  vessel  motion  from  straight  line  motion 
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from: 


Ah  (tj)  =•  a*  Atj  + 


cos(i/)j  At,. 


(4.30) 


0.5  m  /  3600  sec 
(0.11,  0.10,  0.09, 


0.08,  0.07,  0.06)  m 


rp-  =  (27r/10,  2X/13,  2ir/l7 ,  2ir/23,  27r/31,  25r/36) 
-  (20.0;r,  16. 9^-,  IA.Itt,  II.jjt,  9.07r,  S.Stt) 


(4.31) 

(4.32) 

(4.33) 
(4 . 34  j 


The  linear  term  accounts  for  a  height  change  of  0.5  meters  over  a  time  of  1 
hour,  and  the  harmonic  constituents  generate  maximum  vertical  motion 
amplitudes  of  0.5  meters.  The  periods  of  the  vertical  motion  range  from  10 
second  to  36  seconds.  The  total  vertical  motion  profile  is  shown  in  Figure 
4.5.  The  detailed  vertical  position  variation  in  the  first  two  minutes  of  the 
simulation  time  interval  is  depicted  in  Figure  4.6. 

A  data  file  containing  three  dimensional  position  coordinate  changes  Ah(tj), 
Ai^(t.),  AA(tj)  at  one  second  intervals  for  a  total  time  span  of  one  half  hour 
is  input  to  the  GPS  measurement  simulation  software. 

4.1.3  Measurement  Errors 

GPS  measurements  pseudorange,  P,  and  carrier  phase,  F,  can  be  represented 


where : 


P  -  p  +  dp 

$  -  p  +  dj 

p  -  i  1  r  -  R 


(4.35) 

(4.36) 


(4.37) 


is  the  geometric  distance  between  the  satellite  antenna  position  vector,  r,  at 
signal  transmission  time  (as  represented  by  the  satellite  ephemerides)  and  the 
(unknown)  receiver  antenna  position  vector,  R,  at  signal  reception  time.  The 
terms  dp  and  d^  represent  errors  and  biases  in  the  measurements  and  consist  of 


dp  -  dp  +  dj^,,  +  dj  +  c.(dt  -  dT)  +  A,.  +  »?,.+ 


(4.38) 


where : 


+  ^trop  +  c*(dt 


dT)  +  A  +  r)  + 


(4.39) 


is  the  range  error  resulting  from  satellite  ephemerides  errors, 
is  the  range  error  caused  by  (dispersive)  ionospheric  signal 
delay , 

is  the  range  error  caused  by  tropospheric  signal  delay, 
is  the  speed  of  light  in  vacuum, 

is  the  range  error  caused  by  the  satellite  clock  error, 
is  the  range  error  caused  by  the  receiver  clock  error, 
is  the  code  ambiguity, 
is  the  carrier  phase  ambiguity, 

is  the  range  error  caused  by  code  signal  multipath, 
is  the  range  error  caused  by  carrier  signal  multipath. 
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Figure  4.5:  Simulated  vertical  vessel  motion  for  total  simulation  interval 


Figure  4,6:  Simulated  vertical  vessel  motion  for  first  two  minutes 
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and  £p,£,  are  the  range  errors  in  pseudoranges  and  carrier  phases  caused  by 
measurement  noise.  For  a  realistic  simulation  study,  the  simulated  GPS 
measurements  must  be  contaminated  by  artificially  generated  errors.  The 
following  sections  discuss  the  significance  of  the  individual  errors  and  tiieir 
generation  terra  uy  terra. 

4. 1.3.1  Satellite  Ephemeris  Errors 

Broadcast  satellite  ephemerides  errors  result  from  imperfections  in  the 
prediction  process  of  the  satellite  orbits.  The  size  of  these  errors  has  been 
assessed  recently  in  [4-12]  by  comparing  Broadcast  Ephemerides  with  the  more 
accurate  Precise  Ephemerides.  This  comparison  does  not  yield  the  actual 
broadcast  ephemerides  errors,  but  it  tells  us  about  the  order  of  magnitude  of 
these  errors.  The  authors  of  the  aforementioned  study  concluded,  that  the 
maximum  broadcast  ephemerides  errors  for  GPS  satellites  with  Cesium  clocks 
(which  will  be  assumed  for  th\.s  study)  were  of  the  order  of  five  meters  in  the 
along- track  component,  three  meters  in  cross-track  direction  and  two  meters  in 
radial  direction.  These  results  were  obtained  when  ephemerides  were  uploaded 
three  times  per  day.  For  the  present  study,  we  take  the  pessimistic  view  to 
use  these  numbers  as  standard  deviations  for  the  broadcast  ephemerides  errors. 
The  errors  are  taken  to  be  constant  during  the  simulation  time  interval. 

The  geometric  pseudorange  and  carrier  phase  measurements  are  based  on  the 
true  satellite  position  which  is  computed  according  to; 

^true  “  r  +  dr  (4.40) 

with  r  being  the  satellite  position  from  the  ephemeris  file.  The  ephemeris 
error  is  calculated  according  to; 

dr  -  R(0,  5  tn).e3^„^g  +  R(0.  3  ni).e^^^33  +  R(0.  2  (4.41) 

where:  R(qi,  jS)  is  a  normal  distribution  random  number  with  mean  a  and 
standard  deviation  p.  and  are  unit  vectors  pointing  in 

along-track,  cross-track  and  radial  direction  respectively.  The  resulting 
ephemerides  errors  are  constant  for  a  particular  satellite  during  the 
simulation  interval,  but  are  different  from  satellite  to  satellite.  These 
ephemerides  errors  do  not  include  the  orbital  error  caused  by  S/A,  see  section 
4.1.3.11. 

4. 1.3. 2  Satellite  Clock  Errors 

The  position  determination  program  described  in  section  3  of  this  report 
is  based  on  the  analysis  of  double  differences  of  both  pseudorange  and  carrier 
phase  measurements.  In  this  type  of  processing,  the  satellite  clock  error  is 
effectively  eliminated  if  the  measurements  at  the  monitor  station  and  the 
remoue  receiver  are  collected  simultaneously.  Therefore,  the  inclusion  of 
simulated  satellite  clock  errors  in  the  present  study  is  not  necessary. 
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■  *.1.3.3  Ionospheric  Dispersive  Refraction 

The  range  error  caused  by  the  ir>.teracr  ion  of  th.e  dPS  signal  w;:!;  t;.t- 
ionospherc  depends  on  the  total  electron  content,  TKC ,  along  the  path  a; 
JPS  signal  and  the  frequencv  of  the  carrier  signal,  t,  accordini'  ti' 


TEC  depends  on  the  electron  distribution  in  the  annospiiere  ,  wiiicii  :,ot 
only  chan,ges  with  diurnal,  annual  and  solar  (11  vear)  cycles,  but  i,s  al,s<> 
subject  to  short  periodic  variations.  The  high  frequency  var  i .» t  i  on.s  are  o: 
concern  for  signal  cracking,  and  the  correction  of  carrier  piiase 
discontinuities  [a-dl.  Their  effect  on  positions  tends  to  average  out  it; 
time.  Short  periodic  variations  in  the  ionospiiere  will  not  he  further 
discussed  in  this  study.  The  slowly  varying  part  of  the  TEC  has  been  show:,  to 
lead  Co  biases  in  the  estimated  positions  [4-6]  and  is  more  relevant  as  tar  as 
differencial  position  accuracy  is  concerned. 

To  generate  range  errors  due  to  a  slowly  varying  ionosphere,  we  considt:! 
the  following  simplified  model  for  Che  ionosphere  ia-dj.  The  ionospiiere  i  .s 
represented  by  a  chin  spherical  layer  at  a  fixed  altitude  with  a  surface 
density  equal  to  the  total  vertical  electron  content,  VF.C .  This  surface 
density  is  spatially  constant,  and  does  not  change  for  the  measurement 
simulation  time  .span.  From  Figure  4.7  we  see.  that  in  this  model  for  a 
particular  GPS  signal,  the  relation  between  TEC  and  VEC  is 

TEC  -  VEC  (4.43) 

cos  z ' 

where  z'  is  the  zenith  angle  of  the  GPS  signal  at  the  ionospheric  layer.  ' 
can  be  related  to  the  zenith  angle  at  the  surface  of  the  earth,  z,  through: 

sin  z'  -  R  sin  z, 

R+H 

Equations  (4.42)  through  (4,44)  allow  the  simulation  of  ionospheric  delay 
errors,  if  the  height  of  the  ionospheric  layer,  H,  and  VEC  are  specified. 
Typically,  altitudes  between  300  and  400  km  are  chosen  for  an  ionospheric 
surface  layer  model.  For  this  study,  we  shall  use  H  -  350  km.  The  vertical 
electron  density  varies  between  110^*  and  210^®  m'^  depending  on  solar 

activity  and  season  of  the  year  [11].  We  shall  use  VEC  -  0.510^®  m‘^  for 

our  study,  which  yields  an  ionospheric  range  error  of  about  9  meters  for  a 
ver’‘i'al  incident  GTC  signal  and  22  meters  range  error  for  a  .signal  at  IS" 
elevation . 

Some  GPS  receivers  allow  for  the  measurement  of  the  ionospheric  delav 
difference  between  the  two  GPS  signals  by  cross  correlating  the  P-code  bit 
structures  on  LI  and  L2 .  The  measurement  can  be  uced  to  correct  single 
frequency  C/A-code  pseudoranges  for  ionospheric  delay. 
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4. 1.3.4  Tropospheric  Effects 


Tropospheric  delay  is  caused  by  the  lower  non-dispersive  part  of  the 
atmosphere.  It  can  be  modelled  exactly,  if  the  atmospheric  pressure,  the 
atmospheric  temperature  and  the  humidity  are  known  along  the  signal  path. 

Since  these  quantities  can  not  be  measured  directly  along  the  signal  path, 
they  are  approximated  be  an  external  model  for  the  vertical  profile  of  the 
atmosphere,  driven  by  surface  meteorological  data. 

The  total  tropospheric  delay  is  usually  separated  in  the  "dry"  delay  part 
and  the  "wet"  delay  part.  The  dry  delay  depends  on  pressure  and  temperature 
only  and  amounts  to  about  2  meters  in  the  zenith  direction.  It  can  be  very 
accurately  determined  from  surface  meteorological  data.  The  wet  delay  depends 
also  on  the  atmospheric  humidity  and  amounts  to  about  0.3  meters  in  the  zenith 
direction.  Its  relative  accuracy  as  determined  from  surface  met,  data  is  not 
as  high  as  for  the  dry  component,  primarily  because  of  irregularities  in  the 
vertical  humidity  profile. 

The  accuracy  of  the  vertical  tropospheric  delay  as  calculated  from 
surface  pressure,  temperature  and  humidity  is  of  the  order  of  up  to  a  few 
centimeters.  Its  variations  over  time  depend  on  local  weather  conditions  and 
cannot  be  easily  and  realistically  modelled  in  a  simulation.  Tropospheric 
delay  is  not  included  in  this  simulation  study. 

4. 1.3. 5  Receiver  Clock  Error 

The  double  difference  processing  of  pseudoranges  and  carrier  phases  in 
the  position  determination  program  eliminates  receiver  clock  error  completely 
if  the  measurements  to  different  satellites  are  referenced  to  the  same 
receiver  time  tag.  This  is  the  case  in  all  relevant  types  of  presently 
available  GPS  receivers.  Therefore,  the  inclusion  of  simulated  receiver  clock 
errors  in  the  present  study  is  not  necessary. 

4. 1.3. 6  Code  Ambiguity 

GPS  pseudoranges  are  measured  by  aligning  codes  transmitted  on  the 
carrier  signal  with  similar  codes  generated  in  the  GPS  receiver.  For  P-code 
pseudoranges,  this  alignment  can  be  done  unambiguously.  For  C/A-code 
pseudoranges,  the  code  ambiguity  in  the  measured  pseudoranges  equals  an 
integer  multiple  of  the  C/A-code  length  of  1.10'^  sec,  which  is  equivalent  to 
about  300  km.  This  ambiguity  can  easily  be  resolved  and  needs  no  further 
investigation  or  simulation  in  the  present  study. 

Codeless  pseudoranging  utilizes  the  structure  of  the  code  only  and  not 
the  code  itself.  If  pseudoranges  are  measured  codelessly,  their  ambiguity  is 
equal  to  an  integer  multiple  of  the  code  chip  length.  Since  codeless 
pseudoranging  makes  sense  only  with  the  P-code  (C/A-code  is  freely  available) , 
the  ambiguity  in  codeless  pseudoranging  is  a  multiple  of  the  P-code  chip 
length  of  about  30  meters.  This  ambiguity  can  be  easily  resolved  using 
conventional  C/A-code  pseudoranges,  and  will  not  be  further  investigated  or 
simulated  in  this  study. 
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4. 1.3. 7  Carrier  Cycle  Aiabiguity 

The  carrier  phase  ambiguity,  A,  is  a  constant  integer  multiple  N  of  Che 
carrier  wave  length  A  for  code  tracking  receiver  channels ,  and  an  integer 
multiple  of  half  the  carrier  wave  length  for  codeless  tracking  receiver 
channels . 

A  -  N*A,  A  -  N-A/2  (4.45) 

The  correct  resolution  of  the  carrier  phase  ambiguity  is  of  utmost  importance 
in  the  present  study.  To  enable  an  easy  assessment  of  the  correct  carrier 
phase  ambiguity  resolution  in  the  position  determination  software,  all 
simulated  carrier  phase  measurements  are  generated  with  a  unique  ambiguity 

A  -  PRN.IOOO.A,  A  -  PRN*1000*A/2  (4.46) 

for  code  tracking  and  codeless  tracking  receiver  channels,  PRN  is  the 
satellite's  pseudo  random  noise  number. 

4. 1.3. 8  Carrier  Cycle  Slips 

Carrier  cycle  slips  are  changes  in  the  carrier  phase  ambiguity  by  an 
integer  number  of  carrier  wave  lengths  (half  wave  lengths  for  codeless 
tracking  channels) .  They  are  induced  by  a  temporary  loss  of  lock  of  the 
receiver  channel.  To  retain  the  full  accuracy  potential  of  carrier  phase 
measurements,  cycle  slips  must  be  detected  and  corrected.  The  detection  of 
cycle  slips  can  be  achieved  by  a  variety  of  procedures,  and  will  not  be  dealt 
with  in  this  study. 

The  occurrence  of  cycle  slips  can  lead  to  three  different  situations, 
depending  on  how  many  channels  retain  lock  to  the  GPS  signals.  If  loss  of 
lock  occurs  on  all  channels  simultaneously  and  no  inertial  system  is  used,  or 
the  loss  of  lock  lasts  longer  than  INS  can  maintain  the  phase  ambiguity,  all 
carrier  phase  ambiguities  have  to  be  reestablished  in  the  same  way  as  at  the 
beginning  of  a  survey.  This  case  does  not  require  a  separate  simulation.  The 
use  of  an  inertial  system  to  bridge  periods  of  loss  of  lock  to  the  GPS  signals 
is  evaluated  in  Section  5.  Generally,  if  lock  is  retained  for  at  least  four 
channels  and  Perpendicular  Dilusion  of  Precision  (PDOP)  is  reasonable, 
instantaneous  cycle  slip  fixing  for  the  other  channels  is  trivial  (enough 
instantaneous  redundancy  in  the  measurements).  Again,  this  case  does  not 
require  simulation  studies. 

In  the  third  category,  lock  is  maintained  for  less  than  four  channels, 
and  at  least  for  one  channel.  In  this  case,  cycle  slip  correction  is  not 
trivial  and  needs  simulation  in  the  context  of  this  study.  Obviously,  the 
problem  of  determining  these  cycle  slips  is  identical  to  determining  carrier 
phase  ambiguities,  with  some  of  the  carrier  phase  ambiguities  fixed  already. 

In  this  study,  cycle  slip  recovery  is  simulated  assuming  two  double  difference 
ambiguities  (lock  maintained  on  three  satellites)  are  fixed  at  their  proper 
integer  values,  and  all  other  carrier  phase  ambiguities  are  completely 
unknown . 
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4. 1.3.9  Multipath 


Multipath  errors  result  from  the  interference  of  signals  which  have 
travelled  along  signal  paths  of  different  length.  Usually,  one  of  the  signals 
is  travelling  along  the  line  of  sight,  and  the  other  components  result  from 
reflections  at  conducting  material  in  the  vicinity  of  the  GPS  antenna.  The 
changing  geometry  between  satellite,  reflectors  and  antenna  generate.^  cyclic 
multipath  errors.  The  amplitudes,  phases  and  frequencies  of  these  errors 
depend  on  a  number  of  factors  like: 

•  conductivity  of  the  reflecting  material, 

•  shape  of  the  reflecting  material, 

•  distance  from  the  reflecting  material, 

•  geometry  between  satellite,  reflector  and  antenna, 

•  antenna  sensitivity  pattern, 

•  wave  length  (chip  length)  of  signal. 

A  complete  model  containing  these  parameters  is  rather  complicated  and  beyond 
the  scope  of  this  report.  For  the  present  study,  the  follo'wing  simplified 
multipath  error  model  is  assumed: 

7/  -  D  •  cos  (vt  +  w)  (4.47) 

where:  D  is  the  multipath  error  amplitude,  n  is  the  multipath  error  frequency 

and  w  is  the  multipath  error  initial  phase.  The  parameters  adopted  for  the 
present  study  are: 

D  -  R(0,  3.0  m)  for  C/A-code  pseudoranges. 

D  -  R(0,  0.5  m)  for  P-code  pseudoranges. 

D  “  R(0,  1  cm)  for  carrier  phases. 

V  -  27r/300  seconds  (period  of  five  minutes)  . 

w  -  2ir  •  0.0112  •  PRN  for  LI  pseudoranges. 

w  “  2n  •  0.0173  •  PRN  for  L2  pseudoranges. 

w  -  •  0.0133  •  PRN  for  LI  carrier  phases, 

w  -  271'  •  0.0194  •  PRN  for  L2  carrier  phases. 

Obviously,  these  parameters  together  with  equation  (4.47)  model  multipath 
errors  as  simple  cosine  waves  with  five  minute  periods,  and  with  different 
initial  phases  for  individual  satellites.  No  multipath  model  is  employed  for 
measurements  of  differential  ionospheric  group  delays  as  referenced  in  section 
4. 1.3. 3.  This  somewhat  optimistic  assumption  can  be  better  understood  as  more 
field  experience  is  gained  with  cross-correlating  receivers  operating  in  a 
dynamic  environment. 

4.1.3.10  Measurement  Noise 

All  GPS  measurements  are  contaminated  by  noise  resulting  from  the 
tracking  loop  inaccuracies  and  ocher  limitations.  This  measurement  noise  will 
be  different  in  size  frequency  distribution  for  different  receivers.  For  the 
purpose  of  this  study,  we  will  assume  the  measurement  noise  to  be  a  normal 
distribution  white  noise  with  zero  mean.  For  pseudoranges  measured  with  code 
tracking  and  codeless  cracking  receiver  channels,  the  standard  deviation  of 
the  noise  is  chosen  to  be  one  percent  of  the  chip  length  of  the  code  [4-16]. 
These  numbers  may  appear  large  when  compared  with  the  accuracies  quoted  by  GPS 
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receiver  manufacturers,  since  the  instantaneous  code  measurements  are  often 
averaged  over  limited  time  periods  for  output. 

fc'p  “  R(0,  0,3  m) .  for  P-code  pseudoranges  (4,48) 

Sp  “  R(0,  3  ra) ,  for  C/A-code  pseudoranges  (4.49) 

For  measurements  of  differential  ionospheric  group  delays  from  P-code 
cross  correlation,  we  assume  the  same  noise  level  as  for  P-code  pseudoranges. 
For  carrier  phase  measurements,  the  noise  level  is  of  the  order  of  about  half 
a  centimeter  for  a  code  tracking  receiver  channel,  and  a  few  centimeters  for  a 
codeless  cracking  receiver.  For  this  study,  the  simulated  carrier  phase 
measurements  will  be  contaminated  by  measurement  noise  according  to: 

fc'j  -  R(0,  0.5  cm),  for  code  tracking  receivers  (4.50) 

and 

-  R(0,  2  cm),  for  codeless  cracking  receivers.  (4.51) 

4.1.3.11  Selective  Availability 

S/A  is  the  intentional  degradation  of  the  GPS  capability  to  a  level  of 
about  100  meters  horizontal  position  accuracy.  Assuming  a  horizontal  dilution 
of  precision  (HOOP)  between  3  and  ,,  the  range  measurement  errors  required  to 
achieve  this  level  of  S/A  are  of  the  order  of  20  to  30  meters.  These  range 
errors  can  be  introduced  in  two  different  procedures.  The  first  procedure  is 
to  introduce  satellite  clock  errors  yielding  the  appropriate  range  error 
level.  This  type  of  S/A  error  is  eliminated  in  relative  positioning  using 
double  differences  of  pseudoranges  and/or  carrier  ph«''“s  Therefore,  this  S/A 
error  will  not  be  simulated  in  this  study. 

The  second  type  of  S/A  errors  is  caused  by  changing  the  satellite 
ephemerides.  Thi.s  type  of  S/A  error  does  not  completely  cancel  in  relative 
positioning  and  needs  to  be  investigated.  To  assess  the  impact  on  relative 
positioning,  the  ephemerides  errors  (c.f.  section  4. 1,3.1)  in  the  simulated 
measurements  are  changed  from  equation  (4.41)  to  a  level  of  20  meters  in  all 
three  directions. 

dr  -  R(0.  20  m).e^i^  +  R(0,  20  m).e^^„^^  R(0,  20  m).e^^^.3^  (4.52) 

4.1.4  Simulation  Program  Operation 

The  simulation  program  performs  the  following  steps  for  each  simulation  run; 

•  Read  simulation  parameters  from  command  file. 

•  Do  Chebychev  polynomial  approximation  for  satellite  orbits  (computation 

speed ! ) . 

•  Check  visibility  of  satellites  for  each  simulation  epoch. 

•  Read  time  tag  and  remote  receiver  position  from  file. 

•  Compute  satellite  coordinates  for  all  visible  satellites  at  signal 

transmission  times. 

•  Add  epemeridas  error  to  computed  satellite  position. 
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•  Compute  geometric  range  between  satellites  and  monitor  station,  and 

s.*tellite  and  remote  station. 

•  Add  ionospheric  delay  for  pseudoranges. 

•  Subtract  ionospheric  delay  for  carrier  phases. 

•  Add  multipath  errors  for  both  pseudoranges  and  carrier  phases. 

•  Add  ambiguity  terms  to  carrier  phases. 

•  Add  noise  to  both  carrier  phases  and  pseudoranges. 

•  Form  single  differences  (between  monitor  and  remote  station) . 

•  Form  double  differences  (between  satellites). 

•  Write  simulated  measurements  (double  differences  of  pseudoranges  and 

carrier  phases)  and  satellite  coordinates  to  output  file. 

4.2  Position  Determination  Program 

This  section  describes  the  models  behind  a  computer  program  for  the 
computation  of  a  time  series  of  differential  positions  and  i.  corresponding 
time  series  of  carrier  phase  ambiguity  estimates.  In  the  present  context,  the 
main  purpose  of  this  program  is  to  determine  the  integer  values  of  the  carrier 
phase  ambiguities.  The  following  steps  are  executed  in  each  run  of  the 
position  determination  program; 

•  Read  instructions  from  command  file. 

•  Open  files  and  initialize  remote  receiver  position  with  a  standard 

deviation  of  ±  1,000  km.  This  is  replaced  after  the  last  epoch  by 
the  result  of  pseudorange  positioning. 

•  For  each  measurement  epoch; 

•  Read  one  record  of  simulated  data  (satellite  coordinates,  double 

differences  of  pseudoranges  and  carrier  phases). 

•  Form  linearized  observation  equations. 

•  Combine  with  information  retained  from  previous  epoch. 

•  Solve  for  remote  receiver  position  and  carrier  phase  double  difference 

ambiguities . 

•  If  real  number  estimate  for  ambiguities  allows  proper  determination  of 

integers,  fix  the  ambiguity. 

4.2.1  Observation  Models 

The  measurements  utilized  in  the  position  determination  program  are 
double  differences  of  pseudoranges  and  carrier  phases.  As  mentioned  above, 
these  double  differences  are  assumed  free  of  satellite  and  receiver  clock 
errors.  We  obtain  from  the  observation  equations  for  pseudoranges  and 
carrier  phases  (c.f.  section  4.1.3) 

VAP  -  VAp  +  VA(dp  +  d-^  +  d^^^^  +  \  (4.53) 

VA$  -  VAp  +  VA(dp  -  dj^  +  +  A  +  rj)  +  VAe^  (4.54) 

where  VA  designates  the  double  difference  operator.  The  position  information 

of  the  measurements  is  contained  in  the  double  difference  of  the  geometric 
range,  VAp.  Because  of  positive  spatial  correlation,  the  double  differences 
of  range  errors  caused  by  satellite  ephemerides  errors,  VAdp,  and  caused  by 
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ionospheric  refraccion,  are  greatly  reduced  compared  with  the 

corresponding  undifferenced  values.  On  the  other  hand,  those  errors  that  are 
uncorrelated  between  satellites  and  receivers,  are  increased  by  a  factor  of 
two  (there  are  four  measurements  combined  in  a  double  difference).  Usually, 
this  will  be  the  case  for  multipath  induced  errors,  and  for  the  measurement 
noise . 

4 . 2 . 1 . 1  Pseudoranges 

Pseudoranges  obtained  from  C/A-code  correlation  contain  an  ambiguity 
equal  to  an  integer  multiple  of  the  code  length.  As  mentioned  above,  the 
determination  of  this  ambiguity  is  trivial  and  will  not  be  further  discussed 
here.  (A^.  -  0).  Therefore,  the  model  for  both  C/A-code  and  P-code 
pseudorange  double  differences  reads; 

VAP  -  VAp  +  VAej.  (4.55) 

compared  to  equation  (4.53),  the  meaning  of  the  terms  in  equation  (4.54)  has 
been  changed  in  the  following  manner:  The  observation  double  difference  on 
the  left  hand  side  is  corrected  for  all  correctable  errors,  and  the  error  term 
on  the  right  hand  side  contains  all  those  errors  that  cannot  be  explicitly 
corrected. 

4 . 2 . 1 . 2  Carrier  Phases 

Replacing  the  carrier  phase  ambiguity  term  by  an  integer  multiple  of  the 
carrier  wavelength  according  to: 

VAA  -  AVAN,  (4.56) 

we  obtain  the  observation  model  for  double  difference  carrier  phase 
measurements ; 

VA$  -  VAp  +  AVAN  +  VAe,.  (4.57) 

4.2.2  Adjustment  Model 

Equations  (4.55)  and  (4.57)  are  the  basic  non-linear  description  of  the 
measurements  pseudoranges  and  carrier  phases  (on  the  left  hand  side)  in  terms 
of  unknown  parameters  (on  the  right).  The  unknown  parameters  are  the  remote 
receiver  coordinates  (hidden  in  the  double  difference  of  geometric  ranges, 

VAp)  ,  the  double  difference  of  carrier  phase  ambiguities  VAN,  and  the 
unmodelled  residual  carrier  pseudorange  and  carrier  phase  errors,  VAsj  and 
VAfij  respectively.  The  solution  of  theses  equations  in  a  least  squares  sense 
consists  of  the  following  steps; 

•  Linearization  with  respect  to  receiver  coordinates  (to  enable  the  use 

of  matrix  calculus  in  the  solution  process). 

•  Forming  a  linear  system  of  equations  based  on  the  least  squares 

principle,  i.e,  requiring  the  weighted  sum  of  the  estimated  residuals 
to  attain  a  minimum. 

•  Inverting  the  linear  equation  system  to  obtain  an  explicit 

solution  for  receiver  coordinates  and  carrier  phase  ambiguities. 
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All  calculations  are  being  done  in  cartesian  coordinates  in  a  geocentric 
earth- fixed  coordinate  system.  Input  and  output  of  the  position  determinat ioii 
program  are  in  geodetic  (ellipsoidal)  coordinates. 

4. 2. 2.1  Linearized  Observation  Equations 

The  geometric  range  between  satellite  and  receiver  is  a  non-linear 
function  of  satellite  and  receiver  coordinates  according  to: 

r  -  ij  r  -  R  11  -  ; (X  -  X)2  +  (y  -  Y)^  +  (z  -  Z)2]’/2  (4.58) 

where  we  have  denoted  the  cartesian  coordinates  of  the  satellite  bv  lower  case 
letters  and  the  coordinates  of  the  receiver  by  upper  case  letters. 

Linearizing  this  equation  with  respect  to  receiver  coordinates  at  the 
approximate  position  R  -  we  obtain: 

r  -  11  r  -  R^  11  f(d.dR)  ]  ]  r  -  R  )  |  .dR  (4.59) 

where  dR  -  R  -  R“  is  the  difference  between  the  unknown  true  position  and  the 
known  approximate  position  of  the  receiver.  Squares  and  higher  powers  of  dR 
have  been  neglected  in  equation  (4.59).  Denoting  the  partial  derivatives  of 
the  range  with  respect  to  receiver  coordinates  by; 

ag  -  [(Ax/p)  ,  (Ay/p)  ,  (Az/p))  (4.60) 


we  obtain: 

p  -  P°  +  'dR  (4.61) 

with  the  unknown  coordinate  increment; 

dR  -  [dX,  dY,  dZ]^  (4.62) 

Inserting  equation  (4.61)  into  the  non-linear  pseudorange  observation 
equation  (4,55)  we  get  after  re-arranging  the  terms: 

VAP  -  VAPV  -  (^Sr)  •  dR  +  VAcp.  (4.63) 

Since  the  coordinates  of  the  monitor  station  are  known  (cf.  section  4. 1.2.1), 
no  coordinate  increment  for  the  monitor  station  is  included  in  this  equation. 
For  the  same  reason,  only  a  single  difference  (between  satellites)  of  the 
partial  derivatives,  -a„,  needs  to  be  included. 

For  n  pseudorange  (code  phase)  double  differences  observed  simultaneously 
we  obtain  from  the  previous  equation: 

Ip  “  *  dR  +  Bp  (4.64) 

where  Ip  is  the  (n,l)  vector  of  reduced  pseudorange  measurements  containing 
the  quantities  from  the  left  hand  side  of  equation  (4.63),  is  the  (n,3) 
natrix  of  partial  derivatives  of  the  double  difference  measurements  with 
respect  to  remote  station  coordinates  (design  matrix),  and  is  the  (n,l) 
vector  ^  unmodelled  pseudorange  errors. 


Similarly,  the  n  simultaneously  observed  reduced  carrier  phase  double 
difference  measurements  can  be  organized  according  to; 

”  Ajj  •  dR  +  Ajj  •  N  +  .  (4.65) 

The  accuracy  of  the  pseudorange  and  carrier  phase  double  difference 
measurements  is  specified  by  defining  appropriate  covariance  matrices  and 
Cj.  These  matrices  contain  on  the  diagonal  the  variances  of  Bp  and 
respectively,  and  the  corresponding  covariances  in  the  off-diagonal  elements. 
In  our  simulation  study,  all  covariances  were  neglected. 

If  we  have  some  a  priori  information  N°  about  the  carrier  phase  ambiguity 
unknowns  or  about  the  accuracy  of  the  approximate  coordinates  used  in  the 
linearization  process  described  above,  we  can  write  additional  pseudo¬ 
observation  equations: 

N°  =  N  +  e^j.  (4.66) 

The  pseudo-observation  equation  for  approximate  coordinates 
R°  =  R  can  be  re-arranged  to  read: 

0  -  dR  -H  (4.67) 

The  a  priori  information  is  specified  by  defining  appropriate  covariance 
matrices  and  .  These  matrices  contain  the  variances  and  covariances  of 
Cp  and  respectively. 

4. 2. 2. 2  Combining  Code  and  Carrier 

Combining  all  observation  equations  (4.64)  through  (4.67)  yield  the 
following  linear  equation  system: 


L  -  A  •  q  -t-  e 

with  the  observation  vector: 


L  -  [  lp^  1/.  ir\  oM^ 

the  unknowns : 

q  -  [  dR^ 

and  the  design  matrix: 


A 


The  least  squares  solution  of  (4.68)  is  given  by: 
q  -  (A^C'U)’’  a’’  C’l 


(4.68) 


(4.69) 


(4.70) 


(4.71) 


(4.72) 
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with  C  being  a  block-diagonal  matrix; 

C  -  diag  (Cp.  C,,  C^.  Cg).  (4.73) 

In  numerical  applications,  the  special  structure  of  the  design  matrix  (4.71) 
with  zero  and  identity  sub-matrices  is  utilized  for  a  reduction  of  the 
computational  burden  in  the  inversion  according  to  equation  (4.72).  The 
expected  accuracy  of  the  determined  parameters  is  given  by: 

-  (A^C-’a)-''  (4.74) 

i.e.  the  square  roots  of  the  diagonal  elements  of  Cq  are  the  estimated 
standard  deviations  of  the  errors  of  the  parameters  q.  More  information  about 
the  least  squares  procedure  can  be  found  in  any  textbook  on  least  squares 
adj  us  tment . 

4. 2. 2. 3  Parameter  Transition  Model 

The  a  priori  information  on  carrier  phase  ambiguities  and  approximate 
receiver  coordinates  is  obtained  from  a  parameter  transition  model.  A 
parameter  transition  model  describes  how  a  parameter  changes  from  epoch  to 
epoch.  Since  carrier  phase  ambiguities  are  invariant  in  time,  their  parameter 
transition  model  is  simple. 

-  Nj  (4.75) 

and  if  an  estimate  for  N  and  its  associated  covariance  matrix  exists  from 
the  previous  epoch,  it  can  be  easily  incorporated  in  the  solution  f^r  the  next 
epoch  according  to  equation  (4,66),  The  transition  model  for  the  moving 
receiver  position  has  to  be  derived  from  the  assumed  dynamics  model  of  the 
vessel  in  the  particular  environment  (sea,  state,  etc.).  Since  a  very  general 
linear  prediction  model  assuming  constant  velocity  proved  insufficient,  and  a 
more  detailed  model  was  seen  to  limit  the  usefulness  of  the  parameter 
transition  model  to  particular  applications,  the  idea  of  a  transition  model 
for  the  receiver  position  was  abandoned  altogether.  This,  on  the  other  hand, 
makes  the  position  determination  procedure  independent  of  the  dynamics  of  the 
receiver  motion.  It  can  be  used  without  modification  for  any  type  of  vehicle. 

4. 2. 2. 4  Wide-Laning 

Wide-laning  is  a  procedure  to  obtain  from  two  carrier  signals  at 
different  wave  lengths  a  linear  combination  with  a  longer  wave  length  and  at 
the  same  time  retaining  the  integer  nature  of  the  carrier  ambiguity.  For  GPS 
carrier  signals  with  wavelength  and  X2  this  particular  combination  is  given 
by: 


-  (l/A,  -  l/Aj)-’  /A,  -  $2 


(4.76) 


and 


Ay  -  (1/A.,  -  l/Aj)'^  86  cm  is  the  wide  lane  signal  wave  length. 

If  Che  second  frequency  is  observed  with  a  squaring  type  of  receiver,  its 
effective  wave  ength  is  only  half  of  Aj.  In  this  case,  the  wide -lane 


42 


combination  can  be  formed  according  to: 


-  (2/A2  -  1/A,)-’  [2$2/'^2  -$i/A,]  (4.77) 

and 

Ay  -  (2/A2  -  l/A,)-’  «  34  cm 

is  the  wide-lane  signal  wave  length.  Alternatively,  the  wide-lane  can  be 
formed: 

<I>y  -  (2/A,  -  2/A2)-’  [2  $,/A,  -  2*2/A2]  (4-78) 

and: 

Ay  =  (2/A,  -  2/A2)-''  w  43  cm. 


The  reason  for  wide-laning  dual  frequency  carrier  phase  measurements  is 
to  use  the  effective  longer  wave  length  of  the  signal  for  ambiguity 
resolution.  The  non-zero  elements  of  the  design  matrix  Ay  for  ambiguities 
(cf.  equation  (4.65))  are  filled  with  the  effective  wave  length  of  the  carrier 
signal.  The  larger  this  wave  length  is,  the  easier  can  the  integer  values  of 
the  carrier  phase  ambiguities  N  be  determined  in  the  solution  (4.72)  and 
(4.73).  Obviously,  all  of  the  above  wide-lane  combinations  have  wave  lengths 
larger  than  the  single  frequency  GPS  measurements. 


The  major  drawbacks  in  wide-laning  are  increased  measurement  noise,  and  a 
slightly  increased  sensitivity  to  the  ionosphere.  If  P-code  pseudoranges  are 
observed  on  both  GPS  signal  frequencies,  it  is  possible  to  find  a  particular 
linear  combination  of  these  pseudoranges  that  is  the  same  way  affected  by  the 
ionosphere  as  the  86  cm  wide-lane  combination  of  the  carrier  phases  [40]. 

Using  this  particular  pseudorange  combination  together  with  the  86  cm  carrier 
wide  lane  allows  an  ambiguity  resolution  that  is  not  adversely  affected  by 
ionospheric  delays. 


In  a  similar  fashion,  a  C/A-code  pseudorange  and  the  differential 
ionospheric  delay  measured  through  cross  correlation  of  the  P-code  bit 
structure  can  be  combined  into  a  pseudorange  that  is  the  same  way  affected  by 
the  ionosphere  as  the  wide-lane  carrier  phase  combination. 


The  effect  of  the  ionosphere  on  single  frequency  C/A-code  pseudoranges  is 
different  from  the  ionospheric  delay  in  the  34  cm  or  43  cm  carrier  phase  wide- 
lane  combinations.  Therefore,  the  wide-lane  ambiguity  resolution  through 
C/A-code  pseudoranges  will  be  hampered  by  ionospheric  delay. 

4. 2. 2. 5  Ambiguity  Resolution 


Ambiguity  resolution  means  that  the  integer  value  of  the  carrier  phase 
double  difference  ambiguity  terra  (either  single  frequency  or  wide  lane)  is 
determined.  If  this  is  possible,  the  ambiguity  term  N  is  not  an  unknown 
parameter  anymore  and  can  be  removed  from  the  adjustment  procedure  described 
above.  After  removing  the  ambiguities,  GPS  becomes  effectively  a  carrier 
signal  ranging  system  and  individual  kinematic  positions  can  be  determined 
with  the  accuracy  of  the  carrier  phase  measurement.  If  the  integer  values  of 
the  ambiguities  are  determined  incorrectly,  all  positions  will  be  in  error. 
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Therefore,  the  integrity  of  the  ambiguity  resolution  procedure  has  to  bt- 
guaranteed  through  appropriate  measures. 

For  every  measurement  epoch,  the  least  squares  algor itnra  yields  a  rt-ai 
number  estimate  for  the  ambiguity,  and  an  associated  error  standard  deviatio:; 
(the  square  root  of  the  corresponding  diagonal  element  of  the  covariance 
matrix).  The  following  four  criteria  are  used  to  decide  on  the  quality  of  the- 
ambiguity  resolution.  The  ambiguity  is  fixed  if: 

•  the  real  number  ambiguity  estimate  is  within  ±0.25  cycle  of  an 

integer , 

•  the  estimated  error  standard  deviation  of  the  real  number  ambiguity 

estimate  is  less  than  0.1  cycle, 

•  the  scaled  estimated  error  standard  deviation  of  the  real  number 

ambiguity  estimate  is  less  than  0.1  cycle  and 

•  the  difference  between  the  real  number  estimate  and  the  nearest 

integer  is  within  three  times  the  estimated  and  scaled  error 

standard  deviation  for  the  real  number  ambiguity  estimate. 

Obviously,  the  first  criterion  ensures  that  only  real  number  estimates 
close  to  integers  are  considered  for  ambiguity  fixing.  The  second  criterion 
ensures  that  the  nearest  integer  is  the  only  possible  integer  for  ambiguity 
fixing.  The  third  criterion  needs  a  longer  explanation.  The  estimated 
standard  deviation  of  the  real  number  estimates  is  based  on  the  covariance 
matrix  used  to  describe  the  stochastic  model  for  the  measurements.  Since  the 
measurement  covariance  matrix  described  above  does  not  account  for  temporal 
correlations,  the  estimated  error  standard  deviation  of  the  real  number 
ambiguity  estimate  is  based  on  uncorrelated  input  data.  On  the  other  hand, 
the  error  model  for  multipath  described  in  section  4. 1.3. 9  introduces  strong 
positive  temporal  correlation  into  the  measurement  errors.  Neglecting  these 
correlations  yields  overly  optimistic  ambiguity  error  standard  deviations.  A 
simple  procedure  has  been  derived  to  correct  this  problem.  It  can  be  shown, 
that  for  uncorrclated  input  the  least  squares  ambiguity  estimate  accuracy  is 
roughly  inversely  proportional  to  the  square  root  of  the  number  of  measurement 
epochs  processed.  It  also  can  be  shown,  that  for  errors  with  sinusoidal 
components  with  typical  periods,  T,  the  least  squares  ambiguity  estimate 
accuracy  is  roughly  inversely  proportional  to  the  angular  frequency  2;r/T  times 
elapsed  time.  Based  on  these  findings,  the  estimated  error  standard 
deviations  were  scaled  by  multiplying  them  with: 

scale  =  T  n  (4,79) 

2  TT  t 

where  n  is  the  ntuaber  of  measurement  epochs  since  lock  on,  and  t  is  the 
elapsed  time  since  lock  on.  In  this  study  the  measurement  rate  is  1  Hz ,  and 
the  multipath  period  is  5  minutes.  With  these  numbers,  the  scale  factor  is 
greater  than  1  for  elapsed  time  less  than  one  half  hour.  This  simple  scaling 
restores  to  a  great  extent  the  effect  of  the  neglected  temporal  correlations. 
The  fourth  criterion  listed  above  ensures  that  the  real  number  estimate  of  the 
ambiguity,  its  error  standard  deviation,  and  the  nearest  integer  are 
compatible.  Once  all  ambiguities  are  determined  and  fixed  to  integers,  the 
task  of  the  simulation  program  is  completed.  Any  loss  of  phase  lock 
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afterwards  is  assumed  to  be  detected  by  the  receiver.  Depending  on  the  nundn-r 
of  lose  signals,  the  different  procedures  outlined  in  Section  4. 1.3. 8  wovild 
have  Co  be  invoked  then. 

-^.3  Simulation  Results 

This  section  discusses  the  results  obtained  using  the  measurement 
simulation  and  position  determination  programs  described  above. 

4.3.1  Selection  of  Simulation  Parameters 

The  performance  of  GPS  (and  any  other  satellite  positioning  system;  is 
different  for  different  locations.  This  difference  is  caused  by  the  number  or 
visible  satellites  and  the  geometry  between  these  satellites  and  the  user 
position.  Especially  in  those  periods  when  only  the  minimum  number  of 
satellites  required  for  a  position  fix  is  available,  a  weak  satellite  geometry 
can  lead  to  poor  or  even  useless  positioning  results.  A  comolete  simulation 
study  would  have  to  assess  the  system  performance  at  all  possible  user 
locations  during  2A  hours  a  day.  However,  to  keep  the  number  of  simulations 
as  low  as  possible,  a  representative  simulation  location  was  chosen.  Since 
satellite  visibility  also  changes  in  time,  representative  time  windows  were 
chosen.  The  effect  of  certain  measurement  errors  depends  on  the  separation 
between  monitor  station  and  user  site.  Several  representative  separation 
distances  are  chosen.  These  selection  processes  are  described  below.  For  all 
simulations,  an  elevation  cut-off  angle  of  15  degrees  is  assumed,  i.e. 
measurements  Co  satellites  below  15  degrees  of  elevation  are  not  included. 

4. 3. 1.1  Selection  of  Simulation  Location 

To  assess  the  variability  in  the  number  of  satellites  available  in 
different  parts  of  the  U.S.A.,  satellite  visibility  was  predicted  at  the  three 
locations : 


• 

Washington,  D.C. 

(<5 

39°, 

A  -77°) 

• 

Anchorage ,  AK 

er, 

A  -150°) 

• 

New  Orleans ,  LA 

30°, 

A  -90°). 

These  locations  were  chosen  to  represent  an  east  coast  mid-latitude  site,  a 
west  coast  northern  latitude  site  and  a  site  in  the  Mississippi  delta.  For 
each  site,  the  number  of  visible  GPS  and  GLONASS  satellites  was  predicted 
every  half  hour  on  December  30,  1990,  between  0:00  hours  and  24:00  hours  UT. 
using  the  software  package  MacGEPSAL  [29],  Ephemerides  files  for  GPS 
satellites  were  generated  based  on  the  full  24  satellite  constellation  as 
described  in  [12].  Ephemerides  files  for  GLONASS  satellites  are  based  on  the 
constellation  described  in  [6],  Figures  4.8  through  4.13  show  the  prediction 
results.  It  can  bi,  seen,  that  Che  number  of  visible  GPS  and  GLONASS 
satellites  drops  below  5  occasionally,  for  an  elevation  cut-off  angle  cf  15°. 
Each  data  point  represents  a  visibility  prediction  at  30  minute  intervals. 
Obviously,  there  are  no  two  subsequent  epochs  with  less  than  5  satellites. 
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The  visibility  prediction  results  are  summarized  in  the  following 
statistic : 


Washington,  D.C.;  Minimum  number  of  GPS/GLON'ASS/total  satellites:  4/4/8 

Maximum  number  of  GPS/GLONASS/total  satellites:  8/8/16 

Average  number  of  GPS/GLONASS/total  satellites;  6/6/12 

Anchorage,  AK.:  Minimum  number  of  GPS/GLONASS/total  satellites:  4/6/10 

Maximum  number  of  GPS/GLONASS/total  satellites:  9/9/16 

Average  number  of  GPS/GLONASS/total  satellites;  6/7/13 

New  Orleans,  LA:  Minimum  number  of  GPS/GLONASS/total  satellites;  4/5/9 

Maximum  number  of  GPS/GLONASS/total  satellites:  8/9/16 

Average  number  of  GPS/GLONASS/total  satellites:  6/6/12 


Obviously,  there  are  only  minor  differences  in  the  satellite  visibility  at  the 
three  locations  involved.  In  terms  of  minimum  and  average  number  of  visible 
satellites,  Che  Washington,  D.C..  results  are  slightly  inferior  Co  the  other 
locations.  Therefore,  Washington,  D.C.,  was  chosen  to  evaluate  the  kinematic 
GPS  perfonnance . 

A  geost.ationary  satellite  is  permanently  in  view  at  all  three  sites,  if 
placed  in  an  appropriate  equatorial  orbit.  Therefore,  including  a 
geoscationar>  satellite  in  the  simulation  increases  the  number  of  visible 
satellites  by  one. 

4. 3. 1.2  Selection  of  Simulation  Time  Intervals 

It  can  be  seen  in  Figure  4.8,  chat  in  Washington,  D.C.,  the  average 
number  of  visible  satellites  is  six.  This  average  visibility  occurs  for  an 
extended  period  of  around  midnight.  To  assess  the  kinematic  GPS  ambiguity 
resolution  under  these  average  conditions,  the  time  interval  0:00  -  01:00  UT 
is  chosen  for  simulation.  The  same  time  interval  was  kept  for  the  simulations 
involving  a  geostationary  satellite  and  the  GLONASS -type  satellites. 

4. 3. 1.3  Distance  to  Monitor  Station 

Several  errors  in  differential  GPS  positioning  increase  with  the 
separation  distance  between  monitor  station  and  remote  user.  To  assess  the 
resulting  decrease  in  position  accuracy,  and  the  corresponding  increase  in 
convergence  time  for  the  carrier  phase  ambiguity  resolution,  the  simulations 
were  performed  for  several  representative  separation  distances.  Distances 
chosen  are  10  km.  20  km,  40  km  and  100  km.  If  for  a  particular  satellite 
constellation  and  receiver  configuration  an  ambiguity  resolution  at  a 
particular  distance  was  not  successful,  no  attempt  was  made  to  resolve 
ambiguities  for  greater  distances. 


46 


12 
1 1 
10 
9 
8 
7 
6 
5 
4 
3 
2 
1 

.000  4.800  9.600  14.400  19.200  24.000 

Time  (UTO  in  hours  (Day  354,  1990) 

Washington,  O.C. 

Figure  4.8:  Number  of  GPS  satellites  visible  in  Washington,  D.C. 
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Figure  4.9:  Number  of  Glonass  satellites  visible  in  Washington,  D.C 
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Figure  4,10:  Number  of  GPS  satellites  visible  in  Anchorage,  AK 


Ephemer i des  file:  GLONASS . 24 


12 
1 1 
10 
9 
8 
7 
6 
5 
4 
3 
2 
1 


LATITUDE 

61  d  0  m  Os 

LONGITUDE 
-150  d  0  m  0  s 


MASK  ANCLE 
15.0  deg 

DATA  INTERUAL 
30  min 


.000  4.800  9.600  14.400  19.200  24.000 

Time  <UTC)  in  hours  fDoy  364,  1990 > 

Anchorage 

Fieure  4, 1 1 :  Number  of  Glonasr.  satellites  visible  in  Anchorage,  AK 


48 


\2 
J 1 
10 
9 
8 
7 
6 
5 
4 
3 
2 
1 

.000  4.800  9.600  14.400  19.200  24,000 

Time  <UTC)  in  hours  (Day  364,  1990) 

Neu»  Orleans 

Figure  4.12:  Number  of  GPS  satellites  visible  in  New  Orleans,  LA 
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Figure  4.13:  Number  of  Glonass  satellites  visible  in  New  Orleans,  LA 
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4.3.2  GPS  SaCellLces  Alone 


•  GPS  satellites  are  available  in  the  simulation  time  interval. 

4. 3. 2.1  Single  Frequency  C/A-code 

Carrier  phase  double  difference  ambiguities  did  not  converge  to  their 
integer  values  within  15  minutes  for  10  km  distance  between  monitor  and  remote 
station.  Obviously,  the  same  can  be  expected  for  longer  distances. 

Continuing  the  measurements  longer  than  15  minutes  will  not  improve  the 
results,  because  the  difference  in  ionospheric  delay  effects  on  pseudorange 
and  carrier  signal  leads  to  a  divergence  for  single  frequency  observations. 

The  primary  reason  for  this  failure  is  the  size  of  multipath  introduced 
in  the  C/A-code  pseudorange  measurements  (3  meter  error  amplitude).  Probably, 
a  significant  reduction  of  C/A-code  multipath  error  levels  (less  than  1  meter) 
through  appropriate  measures  would  enable  ambiguity  resolution  with  a  C/A-code 
system . 

An  alternative  to  the  procedure  of  combining  pseudoranges  and  carrier 
phases  is  to  use  carrier  phase  measurements  to  multiple  satellites  alone  [24). 
This  procedure  is  free  from  any  pseudorange  multipath  errors,  but  requires 
long  convergence  intervals. 

If  at  least  two  carrier  phase  double  difference  ambiguities  are  fixed  to 
their  integer  values,  the  remaining  unknown  ambiguities  converge  to  their 
correct  integer  values  within  13  minutes  (±  5  minutes)  for  a  distance  between 
monitor  and  remote  station  of  20  km.  Cycle  slip  recovery  after  partial  phase 
lock  loss  is  possible. 

4. 3. 2. 2  Single  Frequency  C/A-code  Plus  Codeless  Second  Carrier 

For  all  simulations  Involving  C/A-code  receivers  with  codeless  second 
frequency,  the  34  cm  wide-lane  wave  length  was  chosen  (cf.  Section  4. 2. 2. 4). 
For  a  discussion  of  the  possible  improvement  with  the  43  cm  wide  lane  wave 
length,  see  Section  4.4  Conclusions. 

Carrier  phase  double  difference  wide-lane  (squaring)  ambiguities  do  not 
converge  to  their  integer  values  within  15  minutes  for  10  km  distance  between 
monitor  and  remote  station.  Obviously,  the  same  can  be  expected  for  longer 
distances.  Continuing  the  measurements  longer  than  15  minutes  will  not 
improve  the  results,  because  the  difference  in  ionospheric  delay  effects  on 

pseudorange  and  carrier  signal  leads  to  divergence  even  larger  than  for  single 
frequency  measurements.  The  primary  reason  for  non-convergence  is  again 
multipath  introduced  in  the  C/A-code  pseudorange  measurements. 

Cycle  slip  recovery  after  partial  phase  lock  loss  was  not  successful  for 
20  km  separation  distance,  primarily  because  of  the  contamination  by 
ionospheric  delay  effects. 
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4. 3. 2. 3  Code  and  Codeless  Dual  Frequency  Pseudorange  and  Carrier  Phase 

For  both  receiver  types,  full  wide-laning  is  possible.  The  remaining 
main  differences  are  the  larger  noise  and  multipath  contamination  of  the 
pseudoranges  in  the  codeless  system. 

For  the  full  P-code  receiver,  the  wide-lane  carrier  phase  ambiguities 
converge  to  their  correct  integer  values  within  6  minutes  (±  2  minutes)  for 
both  10  km  and  20  km  separation  distance,  and  within  7  minutes  (±  2  minutes) 
for  distances  of  40  km  and  100  km  between  monitor  and  remote  station.  Cycle 
slip  recovery  after  partial  phase  lock  loss  was  possible  after  3  minutes  (±  2 
minutes) . 

For  the  P-codeless  type  of  receiver,  the  wide  lane  carrier  phase 
ambiguities  did  not  converge  to  integers  within  15  minutes  for  any  of  the 
above  distances.  The  prime  reason  for  this  failure  is  the  multipath  in  the  LI 
C/A- code  derived  pseudoranges.  Simulations  with  reduced  multipath  amplitudes 
(1  meter  instead  of  3  meters)  yielded  convergence  to  correct  integers  within 
10  minutes  for  10  km  separation  distance,  and  within  15  minutes  for  100  km 
separation  distance. 

4.3.3  GPS  Plus  GLONASS 

A  total  of  12  satellites  are  utilized  in  this  part  of  the  simulation 
study  (6  GPS  satellites  and  6  GLONASS  satellites). 

4. 3. 3.1  Single  Frequency  C/A-code 

Improvement  compared  to  GPS  alone,  but  still  no  convergence  to  integer 
values  for  Che  carrier  phase  double  difference  ambiguities  within  15  minutes. 

If  at  least  two  carrier  phase  double  difference  ambiguities  are  fixed  to 
their  integer  values,  the  remaining  unknown  ambiguities  converge  to  their 
correct  integer  values  within  8  minutes  (±  5  minutes)  for  a  distance  between 
monitor  and  remote  station  of  20  km.  This  is  a  slight  improvement  compared  to 
GPS  alone. 

4. 3. 3. 2  Single  Frequency  C/A-code  Plus  Codeless  Second  Carrier 

Carrier  phase  double  difference  wide-lane  (squaring)  ambiguities  did 
converge  to  their  integer  values  within  about  15  minutes  for  10  km  separation 
distance . 

Cycle  slip  recovery  after  partial  phase  lock  loss  was  successful  for 
20  km  separation  distance  after  about  7  minutes. 


4. 3. 3. 3  Code  and  Codeless  Dual  Frequency  Pseudorange  and  Carrier  Phase 

For  both  receiver  types,  the  additional  measurements  lead  to  improvements 
in  convergence  time  compared  to  GPS  satellites  alone.  For  the  full  P-code 
receiver,  the  wide-lane  carrier  phase  ambiguities  converge  to  their  co.-rect 
Integer  values  within  4  minutes  (±  2  minutes)  for  all  separation  distances 
between  monitor  and  remote  station. 
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For  Che  P-codeless  type  of  receiver,  the  wide-larie  carrier  phase 
ambiguities  did  converge  to  their  correct  integer  values  within  9  minutes  for 
separation  distances  of  10  km  and  40  km.  and  within  10  minutes  for  100  kin 
separation  distance.  Reducing  Che  rauicipath  amplitude  to  1  meter,  convergence 
was  achieved  within  6  minutes  for  distances  of  10  km  and  100  km  between 
monitor  and  remote  station. 

4.3.4  GPS  Plus  Geostationary  Satellite 

4. 3.4.1  Single  Frequency  C/A-code 

Improvement  compared  to  GPS  alone,  but  still  no  convergence  Co  integer 
values  for  Che  carrier  phase  double  difference  ambiguities  within  15  minutes. 

If  at  least  two  carrier  phase  double  difference  ambiguities  are  fixed  to 
their  integer  values,  the  remaining  unknown  ambiguities  converge  to  their 
correct  integer  values  within  11  minutes  (±  5  minutes)  for  a  distance  between 
monitor  and  remote  station  of  20  km.  This  is  a  slight  improvement  compared  to 
GPS  alone. 

4. 3.4. 2  Single  Frequency  C/A-code  Plus  Codeless  Second  Carrier 

No  significant  difference  compared  to  GPS  alone  for  ambiguity 
convergence.  However,  cycle  slip  recovery  after  partial  phase  lock  loss  was 
successful  for  20  km  separation  distance  after  about  8  minutes. 

4. 3.4. 3  Code  and  Codeless  Dual  Frequency  Pseudorange  and  Carrier  Phase 

For  both  receiver  types,  the  measurements  to  the  additional  geostationary 
satellite  lead  to  slight  improvements  in  convergence  time  compared  with  GPS 
satellites  alone.  For  the  full  P-code  receiver,  the  wide-lane  carrier  phase 
ambiguities  converge  to  their  correct  integer  values  within  4  minutes  (±  2 
minutes)  for  10  km  and  20  km  separation  distance,  and  within  5  minutes  for 
separation  distances  of  40  km  and  100  km  between  monitor  and  remote  station. 

As  for  GPS  satellite  alone,  cycle  slip  recovery  was  possible  after  3  minutes 
(±  2  minutes) . 

For  the  P-codeless  type  of  receiver,  the  wide-lane  carrier  phase 
ambiguities  did  converge  to  tleir  correct  integer  values  within  11  minutes  for 
all  separation  distances.  Reducing  the  multipath  amplitude  to  1  meter, 
convergence  was  achieved  within  7  minutes  for  distances  of  10  km  and  100  km 
between  monitor  and  remote  station. 

4.4  Summary  and  Conclusions 

We  have  simulated  the  capability  of  four  different  GPS  receiver 
configurations  for  carrier  phase  ambiguity  resolution  with  simultaneous 
differential  positioning.  The  four  receiver  configurations  were 

•  Single  frequency  C/A-code  receiver. 

•  Single  frequency  C/A-code  receiver  with  codeless  (squaring)  second 

frequency . 

•  Dual  frequency  codeless  pseudorange  and  carrier  phase  receiver  (Rogue 

type) . 

•  Dual  frequency  P-code  receiver  (pseudorange  and  carrier) . 
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It  was  assumed  that  any  of  these  configurations  can  be  used  in  a 
kinematic  environment.  The  simulated  satellite  constellations  included  the 
complete  24  satellite  GPS  configuration,  the  GPS  constellation  plus  one 
geostationary  satellite  and  a  combined  GPS-GLONASS  constellation  with  a  total 
of  48  satellites. 

The  simulation  results  can  be  summarized  as  follows.  Wide- lane 
ambiguities  were  not  successfully  resolved  with  single  frequency  C/A- code 
receivers,  and  with  C/A-code  receivers  providing  a  second  "squared"  carrier 
phase  measurement.  The  convergence  times  for  full  P-code  receivers  and 
P-codeless  receivers  (Rogue  type)  are  summarized  in  Table  4.2.  Results  for 
the  codeless  type  are  provided  for  the  two  pseudorange  multipath  levels  3 
meter  (M  -  3)  and  1  meter  (M  -  1),  Obviously,  multipath  is  the  dominating 
error  in  the  codeless  receiver  type. 

Obviously,  there  is  some  gain  in  convergence  time  by  including  a 
geostationary  satellite  or  GLONASS  satellites  in  the  solution.  For  C/A-code 
single  frequency  receivers  with  codeless  second  frequency,  ambiguity 
resolution  was  not  possible  for  separation  distances  of  20  km  and  more,  and 
even  for  10  km  separation  distance  the  convergence  was  not  achieved  for  all 
tested  satellite  constellations.  Results  were  even  worse  with  a  C/A-code 
single  frequency  receiver  without  the  codeless  second  frequency. 

The  prime  reason  for  the  failure  of  the  C/A-code  receiver  systems  in  the 
present  context  is  the  pseudorange  multipath  model  used  in  the  simulations. 

The  C/A-code  multipath  model  consists  of  sinusoidal  variations  with  amplitudes 
of  3  meters  and  periods  of  5  minutes.  The  convergence  time  can  be  shown  to 
depend  linearly  on  both  the  amplitude  and  the  period  of  the  multipath. 
Therefore,  a  significant  reduction  in  either  of  these  two  parameters  would  be 
required  to  reduce  the  convergence  time  significantly.  Test  runs  with  a 
reduced  multipath  amplitude  of  1  meter  for  20  km  separation  distances  showed 
the  expected  improvement,  but  did  not  allow  complete  ambiguity  resolution 
within  15  minutes. 

The  big  advantage  of  a  dual  frequency  system  is  best  understood  if  we 
look  at  the  ratio  between  wavelength  and  multipath  error  amplitude.  For  the 
simulation  parameters  chosen  in  this  investigation,  this  ratio  is  about  1.7 
for  wide-laning  with  a  full  dual  frequency  P-code  receiver  and  about  0.3  (M=-3) 
or  0.9  (M-1)  for  the  P-codeless  receiver  type,  compared  to  0.11  for  the  34  cm 
wide -lane  with  the  codeless  second  frequency  carrier  described  above.  A 
change  to  the  other  possible  wide- lane  combination  yielding  a  43  cm  wide -lane 
would  improve  this  factor  to  a  value  of  0.14. 
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Table  4.2:  Wide-Lane  Ambiguity  Convergence  Times  (minutes) 
for  Dual  Frequency  receiver  configurations. 


Distance  between  monitor  and  remote 
station 

10km 

20km 

40km 

100km 

GP§  ^lone  (6  sate).lite^)  I 

Full  P-code  Rx 

6 

6 

7 

P-codeless  Rx,  M-3 

no  convergence  any  distance 

P-codeless  Rx,  M-1 

10 

10 

11 

15 

GPS  and  eeostationarv  sat.  il  satellites) 

Full  P-code  Rx 

4 

4 

5 

5 

P-codeless  Rx,  M-3 

11 

11 

11 

11 

P-codeless  Rx,  M-1 

7 

7 

7 

7 

GPS  and  GLONASS  (12  satell 

Ites)  1 

Full  P-code  Rx 

4 

4 

4 

4 

P-codeless  Rx,  M-3 

9 

9 

9 

10 

P-codeless  Rx,  M-1 

6 

6 

6 

6 
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5.0  Measurement  Gap  Filling  Studies 

5 . 1  Requirements  Review 

The  measurement  gap  filling  simulation  was  developed  to  determine  the 
maximiim  gap  lengths  which  could  be  tolerated  before  it  became  necessary  to 
reinitiate  an  integer  search.  High,  medium,  and  low  quality  INSs  were 
integrated  with  receivers  capable  of  wide-laning,  half-wide-laning  and  normal 
laning.  In  addition,  the  use  of  a  barometric  altimeter  and  a  rubidium  clock 
are  considered  to  enhance  the  coasting  performance.  No  specific  goals  (in 
terms  of  minimum  tolerable  gap  length)  were  established;  however,  it  was 
generally  understood  that  gaps  due  to  signal  obscuration  by  bridges  could  be 
up  to  a  minute  in  length. 

5.2  Simulation  Description 

5.2.1  Overview 

An  optimal  covariance  simulation  was  used  to  first  determine  the  extent 
to  which  the  raw  inertial  drift  could  be  calibrated  using  GPS,  and  then  to 
determine  how  the  calibrated  INS  would  perform  with  partial  or  complete  loss 
of  GPS.  Figure  5.2-1  illustrates  the  major  components  of  the  simulation  and 
their  interaction.  The  vehicle  motion  model  and  INS,  GPS  and  baro  altimeter 
error  models  provide  the  necessary  parameters  for  running  the  optimal 
covariance  equations,  which  generate  the  predicted  one  sigma  error  values  for 
each  of  the  modelled  states  as  a  function  of  time.  The  equations  are  first 
run  for  a  60  second  period  with  complete  (i.e,,  4  satellite)  GPS  coverage 
(this  time  period  permits  the  calibration  of  INS  position  and  velocity  errors 
to  reach  a  near  steady  state  condition);  following  this,  a  measurement  gap  of 
specified  length  is  simulated,  during  which  time  a  specified  number  of  GPS 
satellites  are  assumed  obscured;  at  the  completion  of  the  simulated 
obscuration  period,  a  20  second  period  with  full  coverage  is  simulated,  to 
enable  the  filter  to  reestablish  the  correct  integer  count.  Obviously, 
inclusion  of  more  than  the  minimum  number  of  satellites  will  enhance  INS 
calibration  and  so  permit  longer  coasting  periods.  The  simulations  reported 
on  are  limited  to  4  satellites  to  ensure  conservative  predictions. 

The  predicted  one  sigma  values  output  from  the  simulation  represent  the 
performance  of  a  fully  optimal  Kalman  filter,  with  complete  statistical  and 
dynamic  knowledge  of  all  sources  of  error.  These  predictions,  however,  should 
be  only  slightly  optimistic,  since  reduced  state  filters  can  in  general  be 
designed  to  perform  nearly  as  well  as  optimal  filters  through  appropriate 
choice  of  modelled  process  and  measurement  noise  statistics  [10]. 

5.2.2  Simulation  Models 

Each  of  the  models  used  in  the  covariance  simulation  is  described  in  this 
section;  where  appropriate,  equation- level  descriptions  are  provided.  In 
addition,  references  for  the  error  models  and  their  parameters  are  provided 
wherever  possible. 
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Error  Model  Pars. 
Satellite  Geometry 


5. 2. 2.1  Vehicle  Motion  Model 


The  motion  of  the  ship  is  represented  by  straight  line  segments 
(corresponding  to  a  constant  ship  velocity  at  specified  heading)  separated  by 
a  single  turn,  which  is  specified  by  a  fixed  turn  rate  and  duration. 
Superimposed  upon  this  nominal  trajectory  is  a  sinusoidal  heave  motion  of 
specified  amplitude  and  frequency.  For  the  results  reported  in  this  section, 
the  ship  is  initially  located  at  45  degrees  latitude  with  a  45  degree  heading 
and  a  velocity  of  10  knots.  A  turn  is  simulated  after  80  seconds,  completing 
a  course  reversal  after  110  seconds.  Note  that  the  turn  begins  and  is 
completed  during  the  simulated  gap  in  GPS  coverage,  to  make  INS  error 
projections  somewhat  pessimistic.  The  simulated  heave  motion  is  1  meter  in 
amplitude,  with  a  period  of  20  seconds.  The  equations  provided  below  specify 
the  form  of  the  vehicle  motion  model  (note  that  a  spherical  earth  model  is 
used  in  updating  vehicle  latitude  and  longitude) : 

Vehicle  Motion  Model  Equation  Summary 


Given; 

•  ^manstr’  ^man 

If: 

^■manstr  ^  ^^manstr 

^t-nan)  «  - 

H 

man 

Else : 

H  -  0 

Heading  Update  Equation: 

+  H  At 

Latitude  Update  Equation; 

k  ^ 

Longitude  Update  Equation: 

^  sinH^  /  (Re  cosL,^) 

Altitude  Equation: 

hk,T  -  A^  sin  w^,(^,,) 

In  the  above  equations,  H  and  H  denote  vessel  heading  and  heading  rate, 
respectively,  and  At^^^g^  denote  the  start  and  duration  of  the  specified 

maneuver,  respectively,  v  denotes  vessel  speed  and  represents  the  frequency 
of  the  heave  oscillation. 

5. 2. 2. 2  INS  Error  Model 

The  INS  error  model  is  abstracted  from  [39];  it  represents  the  position, 
velocity  and  attitude  errors  in  a  geographic  frame  (i.e.,  east,  north,  up)  and 
includes  both  Schuler  and  earth  rate  error  dynamics.  Both  local  level  and 
strapdown  mechanizations  are  readily  accommodated  in  this  formulation.  The 
dynamics  matrix  for  the  fundamental  (i.e.,  unforced)  error  dynamics, 
abstracted  from  [39j,  is  repeated  here  in  Tables  5.2-1  through  5.2-3.  Two 
minor  changes  have  been  made  to  the  dynamics  matrix  appearing  in  Table  5.2-1: 
the  longitude  and  latitude  errors  were  replaced  by  east  and  north  position 
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Table  5.2- 1  F undamental  Matrix  of  the  INS  General  Error  Differential  Equations 
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-p  /R  tanL/R 


Table  5.2-2  Notation  Used  In  Fundamental  Matrix 


V 

2 

^n' 

f 

2 

Q 

n 

=  Q 

COS  : 

n 

2 

=  Q 

sin  : 

"e 

-  - 

V  /R 
n 

Pn 

=  V 

e 

/R 

=  V 

e 

tan 

'^e 

=  ^e 

w 

n 

+  n 

CO 

2 

+  n 

Latitude  of  vehicle 
Earth  rotation  rate 
Earth  radius 

Magnitude  of  gravity  vector 

Components  of  vehicle  velocity 
with  respect  to  earth 

Components  of  specific  force 
Components  of  earth  rate 


Components  of  angular  velocity 
of  E-N-Z  frame  with  respect  to 
earth 


Components  of  angular  velocity 
of  E-N~Z  frame  with  respect  to 
inertial  space 


V  /R 
z 


42 

2  (Q  V 
n  n 

+  Q 

43 

P  + 
2  ^e 

%  ■ 

44 

~  p  tan  L 

52 

-  2  n 

n 

V 

e 
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OP" 
^n  ^2 

“e 

63 

2g/R  - 

'^n 

92 

^2 

tan 

2  2 


n  n' 


2 

)  V  / 

n  e' 
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Table  5  2-3 


E.-Tor  Forcing  FuncUons  of  the  INS  General  Error 
Error  DtfTercntlal  Equations 


w  =0 
— r  — 

Position  error  components 

w  =  C  5r‘  +  6q 

—V  P  “ 

Velocity  error  components 

X  P 

w  =  C  6u)  . 

-€  p  -ip 

Attitude  error  components 

where 

6fP  =  i p  -  f ° 

Specific  force  measurement 
error  (measured  minus  actual 
platform  components) 

,  .n  "n  n 

62_  =  2.  -  i 

Gravity  model  error  (computed 
minus  actual)  in  geographic 
coordinates,  not  including 
position  error  effect 

-  P  P  "ri  p 

60)^  =  -w‘.  -  u 

— ip  — ip  -ip  — c 

Platform  angular  velocity 
error  (actual  minus  esti¬ 
mated  platform  components 
of  the  angular  velocity  of 
the  platform  in  inertial 
space,  less  any  alignment 
correction  rate) 

Transformation  from  plat- 

P 

form  to  geographic  coordi- 
na  tes 
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errors,  respec  c  ive  Lv ,  which  resulted  in  cerms  involving  the  inverse  of  chi- 
radius  of  che  earch  being  removed  in  the  coupling  to  velocity  and  altitude 
errors,  and  added  in  che  attitude  error  rate  coupling;  in  addition,  ciu*  Ce-rn! 
in  the  6,3  element  of  the  matrix  (rg/R)  which  results  in  the  exponeticial 
growth  in  uncalibraced  IN'S  altitude  error  (caused  by  che  positive  feedback  oi 
altitude  error  into  gravity  computation  error)  has  been  dropped,  sitsce  it  is 
assumed  chat  gravity  compensation  will  be  based  upon  a  constant  .ship  altitud.-. 

In  addition  to  these  basic  9  error  states,  a  total  of  15  gyro  error 
sources,  6  accelerometer  error  sources  and  3  gravity  computation  errors  are 
inodelled,  resulting  in  a  total  of  33  INS  erro->'  states.  The  gyro  error.s 
include  3  g  -  insens  i  cive  drifts,  3  g-sensitive  drifts.  3  gyro  scale  factor- 
errors  and  6  input  axis  misalignments.  Each  g - insens i t ive  drift  is  modeiltd 
as  a  random  walk,  with  a  power  spectral  density  selected  to  produce  an  error 
growth  which  is  10%  of  che  initial  value  over  a  100  second  period.  The 
initial  one  sigma  %'alue  is  a  function  of  the  INS  quality  considered  (see  Table 
5.2-4,  which  suimaar izes  che  gyro  error  model  parameters).  Note  that,  in 
Table  5.2-4  (as  well  as  irr  Table  5.2-5),  the  high  quality  INS  numbers  are 
derived  from  an  error  budget  for  .MAPS  DRU  Ring  Laser  Gyro  (RLG)  based  INS, 
while  Che  medluai  and  low  quality/  numbers  are  representative  of  generic  INSs  of 
io'wer  quality.  Scale  factor  errors  are  modelled  as  Markov  processes  with  1800 
second  correlation  times.  The  relatively  long  time  constant  reflects  the  near 
bias-like  nature  of  che  scale  factor  error,  and  is  consistent  with  models 
presented  elsewhere  [27].  The  process  noise  variance  associated  with  each 
Markov  process  is  computed  from  che  initial  state  error  variance  and  che 
correlation  time  based  upon  steady  state  conditions,  a.s  described  in  [10]. 

The  g-sensitive  drifts  are  expected  to  be  significant  only  for  mechanical 
gyros;  i.e.,  RLGs  are  modelled  with  zero  g-sensicive  drifts.  Thus,  che  three 
qualities  of  INSs  represented  in  Table  5.2-4  are  assumed  to  utilize  RLGs. 

Gyro  input  axis  misalignments  are  represented  as  random  biases. 

Modelled  accelerometer  error  sources  include  accelerometer  biases  and 
scale  factor  errors  (modelled  as  slowly  varying  Markov  processes)  ,  and  input 
axis  misalignments,  represented  as  random  biases.  Table  5.2-5  summarizes  the 
accelerometer  error  source  magnitudes  as  a  function  of  INS  quality.  The 
horizontal  components  of  gravity  computation  error  (often  referred  to  as 
deflections  of  the  vertical)  are  modelled  as  Markov  processes  with  one  sigma 
value  of  10  arc-seconds,  and  correlation  distances  of  25  run,  while  the 
vertical  component  (often  referred  to  as  gravity  anomaly)  has  a  one  sigma 
value  of  50  micro-g's  and  a  correlation  distance  of  60  run.  These  values  are 
generally  consistent  with  levels  which  are  discussed  in  [39).  The  magnitudes 
selected  for  both  the  horizontal  and  vertical  components  may  appear  larger 
than  values  published  elsewhere  [32,23];  this  is  partially  due  to  a  desire  to 
make  che  simulation  performance  predictions  somewhat  conservative,  but  also  to 
include  (approximately)  che  effect  of  bias  components  in  each  gravity  error 
component  (as  referenced  in  [39])  without  the  inclusion  of  additional  states. 

5. 2. 2. 3  Barometric  Altimeter  Error  Model 

Barometric  altimeter  updates  are  processed  each  second  by  Che  filter  to 
assist  in  stabilizing  che  INS  vertical  channel  during  GPS  outages.  Wlien  ful ' 
GPS  coverage  exists,  the  bias  error  of  the  barometric  altimeter  can  be 
calibrated.  The  baro  error  model  consists  of  a  single  Markov  process,  with  a 
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one  sigma  value  of  10  raeCers.  cind  a  correlation  time  of  100  secojids  .  In 
addition,  a  measureitenC  noise  variance  of  0.09  meter  squared  is  as-signed  t.u 
the  baro  measurement  when  it  is  processed  by  the  optimal  filter. 

5. 2. 2. 9  OPS  Measurement  Geometry 

Due  CO  the  relatively  short  time  periods  of  interest,  GPS  satellite 
geometry  is  assumed  fixed,  with  the  four  satellite  azimuths  and  elevations 
indicated  in  Table  5.2-6,  below.  These  satellite  Lines  of  sight  are  derived 
from  real  data  obtained  on  30  April  1983.  from  GPS  satellites  9,  11,  13  and 
19.  This  sec  of  four  satellites  (at  0600  UTC)  results  in  a  PDOP  of  2.5,  and 
HDOP  of  1.5  (the  Vertical  Dllusion  of  Precision  (VDOP)  is  therefore  2.0}. 
Thus,  a  good  nominal  GPS  geometry  is  assumed. 


Table  5.2-6:  Assumed  GPS  Satellite  Geometry 


Simulation  No. 

GPS  Satellite  No. 

Azimuth  Angle  <°) 

Elevation  Ang'e(°) 

SVl 

9 

170 

40 

SV2 

11 

40 

30 

SV3 

13 

340 

66 

SV4 

19 

310 

29 

5. 2. 2. 5  GPS  Error  Model 

The  GPS  error  model  includes  receiver  clock  phase  and  frequency  error, 
and  a  maximum  of  5  error  states  for  each  satellite;  these  errors  include 
carrier  multipath  error,  residual  (i.e.,  post-differentially  corrected) 
tropospheric,  ionospheric  and  S/A  errors,  and  integer  ambiguity  error.  Thus, 
a  maximum  number  of  22  states  can  be  modelled  for  4  satellites.  GPS  carrier 
pseudorange  measurements  are  processed  each  second  with  a  one  sigma  noise 
value  of  1  cm. 

The  receiver  clock  is  assumed  to  be  a  rubidium  frequency  standard  (a 
sensitivity  analysis  is  reported  on  later  in  which  a  non  atomic  clock  is 
considered).  Its  error  model  is  derived  from  [28],  and  consists  of  a  random 
walk  frequency  error  with  an  initial  one  sigma  value  of  lxlO‘^\  driven  by  a 
white  noise  with  power  spectral  density  computed  from  an  Allen  variance  of 
5x10'^^  over  2000  seconds.  In  addition,  a  phase  noise  with  a  power  spectral 
density  l.xlO'^^  sec^sec  is  included.  Computation  of  the  appropriate  power 
spectral  densities  from  the  Allen  variance  statistic  is  discussed  in  [35]. 

The  selection  of  a  random  walk  in  lieu  of  a  Markov  process  is  based  upon  the 
fact  that  frequency  correlation  times  on  the  order  of  several  days 
(corresponding  to  aging  effects)  are  typical  for  rubidium  standards,  which  is 
not  a  significant  effect  for  the  simulation  durations  considered.  Due  to  the 
relative  magnitudes  of  the  phase  and  frequency  power  spectral  densities  and 
the  simulation  duration,  the  phase  noise  is  expected  to  be  the  dominant  error 
contributor . 

Carrier  multipath  error  is  modeled  by  a  Markov  process  with  a  5  minute 
correlation  time  and  a  one  sigma  value  of  1  cm;  these  values  are  consistent 
with  roof-top  tests  performed  using  TNL  4000  (i.e.,  survey  quality)  receivers 


i'22].  Residual  tropospheric,  ionospheric,  and  S/A  errors  are  computed  as  a 
function  of  the  separation  distance  between  the  ship  and  the  (master; 
reference  station  using  the  equations  below: 


Residual  Tropospheric  Error 


(j^ 

a2 


tropo 

tropo 

tropo 

tropoo 


^ropo 

~tropo 

a 

'‘^tropo 

tropo 


'^tropoo  (l-e‘2‘^^‘*tropo) 

’^‘^tropo  ^^tropo  ^  *^tropo 

calculated  tropospheric  error  variance 

non  differentially  corrected  tropospheric  error  variance 
distance  from  master  reference  station 
correlation  distance  for  tropospheric  error 
time  rate  of  change  of  tropospheric  error  variance 
inverse  of  tropospheric  error  correlation  time 
power  spectral  density  for  tropospheric  error  model 
(8x10.2  craVsec) 


Residual  Ionospheric  Error 
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^ono 
^lono 
^lono 
\  onoo 
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^ono 

^iono 

'^iono 


=  (7.  •  d  /  d- 
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=»  -25.  cr.  +  q. 

-  calculated  ionospheric  error  variance 

”  non  differentially  corrected  ionospheric  error  variance 
“  distance  from  master  reference  station 
“  correlation  distance  for  ionospheric  error 

-  time  rate  of  change  of  ionospheric  error  variance 

-  inverse  of  ionospheric  error  correlation  time 

-  power  spectral  density  for  ionospheric  error  model 

(1.8x10*^  mVsec) 


Residual  S/A  Error 
"^sao 

”  '^^sa^sa  ^sa 

=  calculated  sa  error  variance 

”  non  differentially  corrected  sa  error  variance 
=  distance  from  master  reference  station 
“  correlation  distance  for  sa 

-  time  rate  of  change  of  sa  error  variance 

-  inverse  of  sa  correlation  time 

-  power  spectral  density  for  sa  (8  x  10’^  m^/sec) 

Table  5.2-7  summarizes  the  parameters  of  these  error  models.  The  values 
utilized  for  tropospheric  and  ionospheric  error  are  consistent  with  parameters 
published  in  the  open  literature  [3,13,19,18],  The  parameters  used  for  the  SA 
error  model  are  representative  of  only  the  so-called  epsilon  term,  since  the 
dither  term  (resulting  from  the  perturbation  of  each  satellite's  oscillator) 
is  common  to  both  the  ship  and  the  reference  station,  and  will  drop  out  at  all 
separation  distances.  The  residual  error  due  to  the  epsilon  term,  since  it 
represents  a  satellite  position  error,  is  strictly  a  function 
of  ship/station  relative  geometry,  which  behaves  linearly  with  separation 
distance  for  distances  which  are  small  compared  with  the  distance  to  the  GPS 
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satellite.  The  constant  relating  the  two  is  roughly  the  inverse  of  the 
distance  to  the  satellites  (i.e.,  20000  km).  The  very  long  correlation  time 
associated  with  the  S/A  error  reflects  the  slowly  varying  nature  of  the 
induced  orbital  errors. 

Integer  ambiguity  error  is  modelled  as  a  random  bias,  with  a  one  sigma 
value  dependent  upon  the  simulation  conditions.  It  is  initialized  to  zero, 
assuming  that  a  correct  set  of  integers  has  been  found;  during  a  period  of 
continuous  satellite  coverage,  it  remains  zero,  assuming  that  carrier  lock  is 
maintained  without  cycle  slips  occurring.  Following  a  simulated  measurement 
gap,  it  is  set  to  a  one  sigma  value  of  100  meters,  reflecting  no  knowledge  of 
the  integer  count  (other  than  through  the  pseudorange  measurement,  assumed  to 
be  S/A-corrupted) .  As  the  carrier  pseudoranges  (some  with  unresolved 
integers)  are  processed  and  the  one  sigma  error  associated  with  the  integers 
decrease,  a  test  is  made  to  see  if  the  integer  can  be  resolved.  As  a  function 
of  the  receiver  utilized  (i.e.,  wide,  half  wide  or  normal  laning) ,  the  integer 
is  considered  to  be  resolved  if  2.5  times  the  one  sigma  error  is  less  than  one 
half  of  the  lane  width.  This  test  corresponds  to  a  0.99  probability  of 
correct  determination  of  the  integer.  It  is  difficult  to  compare  this 
criterion  directly  with  the  criteria  used  in  the  ambiguity  resolution 
simulation,  since  that  was  a  Monte  Carlo,  rather  than  a  covariance  simulation. 
Certainly,  the  covariance  test  utilized  here  is  more  conservative  than  the 
covariance  test  used  in  the  ambiguity  resolution  simulation.  However,  the 
additional  tests  performed  there,  which  do  not  have  a  counterpart  in  the 
covariance  simulation,  may  make  the  two  sets  of  tests  for  lane  convergence 
more  nearly  equivalent.  Once  an  integer  can  be  resolved,  an  artificial 
perfect  measurement  of  the  integer  error  is  formed,  which  results  in  the  one 
sigma  error  associated  with  that  state  being  set  to  zero,  and  all  errors 
correlated  with  that  state  being  correspondingly  reduced.  Note  Chat, 
alternatively,  integer  error  could  be  modelled  fractionally.  Given  that  a 
conservative  measure  of  correct  integer  determination  is  utilized  (which  is 
the  case  here),  integer  modelling  is  preferred,  since  filter  estimation  errors 
will  generally  be  reduced.  Once  any  integer  has  been  resolved,  the  carrier 
phase  measurements  are  reprocessed,  making  use  of  this  information.  The 
process  (at  each  time  instant)  is  concluded  when  no  new  integers  are  resolved 
in  the  current  measurement  epoch.  Note  that  a  separate  integer  ambiguity  is 
modelled  for  each  satellite,  while  double  differencing  algorithms  consider 
pair-wise  differences  of  these  integers.  Since  a  rubidium  frequency  standard 
is  part  of  the  recommended  system,  the  use  of  double  differencing  (which 
removes  the  effect  of  user  clock  error)  will  not  change  the  performance 
predictions  significantly.  If  operating  without  an  atomic  clock,  the  use  of 
additional  satellites  (beyond  the  minimum  4  utilized  in  the  simulations 
reported  in  this  section)  will  enable,  through  double  differencing  or 
processing  by  the  optimal  filter,  removal  of  receiver  clock  errors.  The 
Kalman  filter  will  effectively  make  use  of  double  differencing  when  presented 
with  redundant  satellites  and  a  clock  with  excessive  drift  characteristics. 
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5.3  Simulation  Results 


5.3.1  Scenario  Deiinition 

Given  the  error  models  just  described,  measurement  gaps  of  different 
durations  were  simulated  following  an  assumed  60  second  calibration  period 
with  full  GPS  coverage  (i.e.,  carrier  pseudorange  measurements  available  r rom 
4  satellites  each  second) .  This  period  was  previously  determined  to  be 
sufficiently  long  to  enable  the  filter  to  achieve  a  near  steady  state  error 
condition.  Except  when  specifically  noted,  the  barometric  altimeter  was 
assumed  available  in  all  cases. 

5.3.2  Cases  Considered 

Of  primary  interest  in  running  the  gap  filling  simulation  is  the 
determination  of  the  maximum  gap  length  which  can  be  tolerated  without 
permanent  loss  of  the  integer  ambiguities  over  the  ^i.ll  range  of  operating 
conditions.  This  includes  separation  distances  ranging  from  0  to  100  km, 
high,  medium  and  low  quality  INSs,  wide-laning,  half -wide- laning  and  normal 
laning  receivers,  and  loss  of  from  1  to  4  GPS  satellites  during  the  simulated 
gap.  The  zero  separation  distance  was  included  as  a  limiting  case  for  which 
the  residual  correlated  errors  could  be  completely  removed.  Due  to  the  large 
number  of  runs  involved,  the  simulation  was  configured  to  automatically 
(without  operator  intervention)  select  the  maximum  tolerable  gap  length  by 
iterating  from  an  initial  supplied  length.  It  tests  the  current  length,  then 
either  doubles  the  gap  length  or  halves  the  previous  change  in  gap  length, 
based  upon  whether  the  integers  were  maintained  or  not,  respectively.  To 
better  understand  the  summary  curves  presented  in  this  section,  consider 
Figures  5,3-la.  and  b.,  which  plot  the  one  sigma  errors  along  each  satellite 
line  of  sight  as  a  function  of  time,  starting  with  the  simulated  measurement 
gap.  As  illustrated  in  Figure  5.3-la.,  the  initial  errors  are  very  small 
follow...ng  the  60  second  calibration  period  with  full  GPS  coverage.  A  gap  of 
50  seconds  with  complete  loss  of  GPS  is  considered,  using  a  high  quality  INS 
and  wide-laning  GPS  receiver.  The  separation  distance  for  this  case  was 
100  km.  Note  that  each  of  the  one  sigma  errors  along  the  4  lines  of  sight 
converges  to  less  than  0.2  of  a  wide-lane  wavelength  (86  cm),  indicating 
reacquisition  of  the  correct  set  of  integers  for  all  4  satellites.  Figure 
5. 3- lb.  is  a  blowup  of  the  integer  reacquisition  sequence  from  Figure  5. 3- la,, 
illustrating  how  the  second  carrier  measurement  enables  initial  separation  of 
INS  position  errors  and  integer  errors,  and  Che  integer  "snap  in"  below  the  1% 
confidence  level. 

Of  secondary  interest  in  the  utilization  of  this  simulation  is  the 
resultant  navigation  performance  of  the  integrated  GPS/INS  navigation  filter 
(navigation  performance  is  addressed  in  Section  5) .  To  illustrate  the 
differences  between  navigation  performance  and  integer  convergence,  consider 
Fig.  5,3-2,  which  plots  the  position  error  components  corresponding  to  the 
line  of  sight  error  components  presented  in  the  previous  figure.  Note  the 
significant  differences  in  the  error  level  magnitudes. 

Position  error  components  prior  to  the  gap  and  following  integer 
reconvergence  range  from  about  0,2  meter  to  more  chan  0.5  meter,  while  the 
integer  ambiguity  errors  along  each  line  of  sight  are  a  fraction  of  a 
centimeter.  This  apparent  inconsistency  (which  is  more  than  an  order  of 
magnitude  larger  chan  the  scaling  attributable  to  PDOP) ,  is  caused  by  the  near 
unity  negative  correlation  coefficient  which  develops  between  position  errors 
and  correlated  errors  in  the  measurements  (i.e.,  residual  tropospheric  and 


68 


Figure  5.3- la,  ReacqulsiUon  of  Widelane  Integers  with  INS  after  50  Second  Gap 


Time  (secs) 


Figure  5.3- lb.  RcacqulslUon  ofWidclanc  Integers  with  INS  after  50  Second  Gap 


ime  (secs' 


Figure 


ionospheric  errors  and  S/A).  In  simpler  terms,  the  sum  of  these  residual 
errors  and  the  position  errors  along  each  line  of  sight  can  be  estimated  at  a 
level  below  the  one  sigma  measurement  noise,  but  it  is  not  possible  to 
separate  the  effects  of  the  two.  This  is  thus  representative  of  a  classical 
observability  problem  [10],  Nevertheless,  the  filter  is  able  to  correctly 
resolve  the  integer  ambiguities  following  the  gap,  even  with  relatively  large 
uncertainties  in  position;  a  correct  integer  does  not  reflect  absolute  range 
information  along  a  given  line  of  sight,  but  includes  the  effects  of  residual 
correlated  ranging  errors,  which  induce  the  positioning  error  in  the  resulting 
solution. 

5.3.2. 1  Complete  Loss  of  GPS 

Results  for  complete  loss  of  GPS  are  summarized  by  Figure  5.3-3.  The 
maximum  tolerable  gap  length  (in  seconds)  is  plotted  on  the  vertical  axis, 
while  the  separation  distance  (in  kilometers)  is  the  horizontal  scale.  A 
total  of  9  individual  curves  appear  on  the  plot,  corresponding  to  all  possible 
combinations  of  high,  medium  and  low  quality  INSs  (designated  by  the  use  of  a 
solid  line,  a  dashed  line  and  a  dotted  line,  respectively)  and  wide,  half -wide 
and  normal  laning  receivers  (designated  by  squares,  diamonds  and  triangles, 
respectively).  Note  the  very  significant  variations  in  performance  with  INS 
quality;  a  high  quality  INS  is  able  to  tolerate  gaps  between  60  and  70  seconds 
when  operating  with  a  wide-laning  receiver,  while  a  low  quality  INS  can  at 
best  tolerate  a  gap  of  less  than  20  seconds.  Also  note  that  the  variation  of 
gap  length  with  separation  distance  is  relatively  minor;  a  maximum  variation 
over  the  full  range  of  separation  distances  of  roughly  12  seconds  occurs  for 
the  high  quality  INS  operating  with  the  normal  laning  GPS  receiver.  This  is 
not  surprising,  based  on  the  discussion  at  the  end  of  the  previous  section. 

The  larger  separation  distances  significantly  degrade  positioning  performance, 
but  the  ability  to  resolve  the  integers  is  based  upon  the  total  error  along 
each  line  of  sight,  which  is  dominated  by  the  INS  position  error  growth  during 
Che  gap  (this  error  is  of  course,  uncor related  with  the  residual  correlated 
ranging  errors). 

5. 3. 2. 2  Three  GPS  Satellites  Tracked  During  Gap 

Results  when  three  satellites  are  still  being  tracked  during  the  gap  are 
summarized  in  Fig.  Note  the  significantly  different  behavior  from  the 

previous  case;  much  longer  measurement  gaps  can  generally  be  tolerated, 
especially  for  the  wide-laning  receiver;  significant  variation  of  gap  length 
with  separation  distance  occurs;  and,  for  a  given  receiver  type,  results  are 
relatively  insensitive  to  INS  quality.  The  longer  gap  lengths  can  be 
attributed  to  the  filter  reliance  upon  the  rubidium  frequency  standard  as  a 
"pseudo  satellite"  to  maintain  an  accurate  clock  during  the  gap.  The  maximum 
tolerable  gap  length  is  therefore  no  longer  strongly  dependent  upon  INS 
quality.  In  fact,  maximum  tolerable  gap  lengths  of  more  than  200  seconds 
should  be  attainable  without  an  INS.  The  more  significant  variation  with 
separation  distance  can  be  attributed  to  the  temporal  decorrelation  of  the 
residual  correlated  errors,  since  the  gap  lengths  for  this  case  (especially 
for  the  wide-laning  receiver,  which  exhibits  the  more  pronounced  variation 
with  separation  distance)  are  now  significant  relative  to  the  correlation 
times  of  the  errors. 
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Figure  5.3-3  Gap  Length  v.  Separation  Distance  Loss  of  Four  GPS  Satellites 
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5. 3, 2. 3  Two  GPS  Satellites  Tracked  During  Gap 

Results  for  2  GPS  satellites  being  cracked  during  the  gap  are  presented 
in  Figure  5.3-5.  Virtually  all  the  trends  which  appeared  in  Figure  5.3-3  for 
complete  loss  of  GPS  are  apparent:  the  results  are  relatively  insensitive  to 
variations  in  separation  distance,  but  differ  significantly  as  a  function  of 
INS  quality.  The  only  difference  from  Figure  5.3-3  is  the  small  improvement 
in  tolerable  gap  length  for  all  cases,  e.g.,  about  10%  for  the  high  quality 
INS  operating  with  wide-laning. 

5. 3. 2. A  One  GPS  Satellite  Tracked  During  Gap 

Results  for  this  case  are  presented  in  Figure  5.3-6,  and  are 
approximately  midway  between  the  results  for  loss  of  2  satellites  and  4 
satellites.  All  the  trends  appearing  in  Figures  5.3-3  and  5.3-5  are  therefore 
repeated  in  Figure  5.3-6. 

5. 3. 2. 5  Operation  Without  the  Barometric  Altimeter 

Although  the  barometric  altimeter  is  expected  to  be  a  relatively  low  cost 
item,  it  nevertheless  represents  an  additional  hardware  component  which  must 
be  integrated  into  the  system,  which  can  be  costly  and/or  troublesome.  Note, 
however,  Chat  even  chough  the  barometric  altimeter  may  prove  to  have  an 
insignificant  role  in  coasting  during  GPS  outages,  it  is  probably  required  as 
a  meteorological  sensor  for  generation  of  the  tropospheric  model  corrections. 
An  additional  sec  of  simulations  were  performed  to  determine  the  tolerable  gap 
length  versus  separation  distance  without  barometric  altimeter  updates  for  a 
complete  loss  of  GPS  during  the  gap.  Comparison  of  these  results  with  the 
corresponding  results  from  Figure  5,3-3  indicated  that  there  was  virtually  no 
difference  between  the  two  (differences  in  gap  length  were  generally  less  chan 
1  second) .  This  implies  that  Che  INS  vertical  drift  has  been  sufficiently 
calibrated  during  the  assumed  60  second  calibration  period  prior  to  the  gap. 
Since  the  case  corresponding  to  a  complete  loss  of  GPS  was  expected  to  produce 
Che  most  significant  differences  when  operating  without  the  baro ,  there  does 
not  seem  to  be  any  motivation  for  including  the  baro  altimeter  in  the  system. 
Only  when  the  INS  vertical  channel  has  not  been  accurately  calibrated  should 
the  baro  input  be  expected  to  affect  performance.  In  this  case,  as  discussed 
in  the  next  paragraph,  the  use  of  a  zero  mean  sea  level  input  to  an 
appropriately  designed  complementary  filter  should  produce  the  same  effect  as 
the  baro  input. 

In  place  of  the  barometric  altimeter,  altitude  error  growth  can  be 
bounded  by  use  of  a  zero  mean  sea  level  input,  which  can  be  used  in  a  manner 
similar  to  the  conventional  damping  of  an  INS  vertical  channel  using  an 
externally  supplied  altitude.  Figure  5.3-7  illustrates  the  concept;  INS 
vertical  acceleration,  corrected  by  the  last  available  update  from  the  Kalman 
filter,  is  compensated  by  the  feedback  of  the  difference  between  computed  ano 
reference  altitude,  and  then  integrated  to  produce  a  vertical  velocity 
estimate.  Vertical  velocity  is  similarly  corrected  by  feedback  of  the 
difference  between  indicated  and  reference  altitude,  and  then  integrated  to 
produce  indicated  altitude.  Since  the  loop  (as  indicated)  is  second  order, 
the  loop  gains  can  be  expressed  directly  in  terms  of  the  resultant  damping 
ratio  and  natural  frequency,  as  indicated  in  the  Figure.  (A  third  order  loop 
can  be  used).  The  loop  represents  what  is  generally  referred  to  as  a 
complementary  filter:  the  INS  is  used  as  a  source  of  high  frequency 
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Figure  5.3-5  Gap  Length  v.  Separation  Distance  Loss  of  Two  of  Four  GPS  Satellites 
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Figure  5  3-7  Zero  Mean  Sea  Level  Damping  of  Altitude  Error  Growth 
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information,  and  so  will  be  used  to  crack  the  deviations  from  steady  state 
caused  by  ocean  wave  motion;  the  steady  state  solution  for  altitude  will  be 
that  supplied  by  the  reference.  The  performance  of  the  resultant  scheme 
should  therefore  very  closely  follow  the  performance  of  a  system  equipped  with 
a  barometric  altimeter. 

5.4  Summary  and  Conclusions 

When  more  than  a  single  GPS  satellite  is  lost,  the  use  of  a  high  quality 
INS  can  significantly  increase  the  maximum  tolerable  gap  length,  permitting 
rapid  reacquisition  of  the  correct  sec  of  integer  ambiguities.  Maximum  gap 
lengths  range  from  10  to  80  seconds  for  the  high  quality  INS,  as  a  function  of 
separation  distance  and  receiver  type.  Of  these,  receiver  type  plays  a  more 
important  role  chan  separation  distance.  Variation  with  separation  distance 
for  a  given  INS  quality  and  receiver  type  is  relatively  minor;  at  most,  15 
seconds  over  the  range  from  0  to  100  km  is  observed.  On  the  other  hand, 
variations  with  receiver  type  are  significant  for  all  INS  qualities,  e.g..  a 
high  quality  INS  can  withstand  a  gap  of  70  seconds  when  operating  with  a 
wide-laning  receiver,  but  the  length  is  reduced  to  roughly  40  seconds  when 
using  a  half -wide - laning  receiver. 

Significantly  different  trends  are  observed  when  only  a  single  satellite 
is  lost.  Here,  only  minor  differences  can  be  attributed  to  INS  quality  and 
greater  variations  with  separation  distance  occur  (especially  for  the  longer 
gap  lengths) .  Maximum  gap  lengths  of  more  than  250  seconds  can  be 
accommodated,  due  to  the  assumed  use  of  a  rubidium  frequency  standard. 

Use  of  a  barometric  altimeter,  as  an  aid  for  gap  filling,  is  not 
necessary.  When  the  INS  vertical  channel  has  been  given  sufficient  time  with 
full  GPS  coverage  (i.e.,  60  seconds),  no  significant  performance  benefits  are 
derived  from  the  incorporation  of  the  baro  inputs.  Periods  followirg 
inaccurate  calibration  of  the  baro  can  be  handled  using  a  zero  mean  sea  level 
damping  input  to  a  complementary  filter  (Figure  5.3-10). 
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6.0  Navigation  Performance  Prediction 

6.1  Discussion  of  Differential  GPS  Operation 

In  this  section,  the  navigation  performance  of  the  kinematic  phase 
positioning  system  is  predicted.  Tlie.  error  budget  proposed  here  is 
essentially  the  budget  for  high-precision  differential  GPS,  with  identical 
error  sources:  differencial  troposphere,  differencial  ionosphere,  computed 
satellite  position,  multipath  and  tracking  loop  error.  Resolution  of  the 
integer  or  lane  ambiguity,  using  techniques  previously  described,  is  assumed 
to  have  been  successfully  performed.  Although  knowledge  of  the  lane  integers 
implies  a  pseudorange  accuracy  of  better  than  a  half -wavelength ,  this  does  not 
necessarily  translate  directly  into  differencial  position  accuracy  of  the  same 
order.  When  the  mobile  receiver  is  far  from  the  reference  receiver,  the 
differencial  position  accuracy  is  affected  by  differential  tropospheric  and 
ionospheric  retraction  and  miscompucacion  of  the  satellite  position,  as  well 
as  carrier  phase  noise  and  multipath  at  the  individual  receivers.  Although 
most  of  the  GPS  error  due  to  these  error  sources  is  common  to  all  receivers  in. 
a  wide  area,  over  tens  of  kilometers  Che  effect  of  each  error  changes 
slightly.  This  is  often  referred  to  as  spatial  decorrelation  of  the 
pseudorange  error. 

Since  much  of  the  change  in  the  pseudorange  error  is  linear  over 
distance,  spatial  decorrelation  effects  can  to  a  certain  extent  be  removed  by 
a  network  of  reference  receivers.  The  efficacy  of  networked  corrections  is 
addressed  in  this  portion  of  the  study  as  well  as  the  accuracy  of  differencial 
corrections.  To  illustrate  the  concept  of  networked  corrections,  consider  a 
master  reference  station  and  two  additional  "monitor"  reference  stations,  one 
directly  to  the  north  of  the  master  and  another  directly  to  Che  east.  Using 
Che  two  master/monitor  pairs  of  reference  stations,  Che  change  in  pseudorange 
error  with  respect  to  distance  in  each  direction  can  be  measured  in  real  time 
and  broadcast  with  the  pseudorange  corrections. 

The  user  adjusts  the  received  differential  corrections  for  his  own 
location  by  applying  a  simple  linear  function  using  the  measured  gradients  as 
coefficients  to  the  difference  between  the  user  position  and  the  master 

reference  position  • 

5PRC/5N  -  gradient  of  PRC  in  northerly  direction 
5PRC/5E  -  gradient  of  PRC  in  easterly  direction 

PRC^ser  -  PRC''®^  +  (dFRC/dli,  dV^C/dE.  0)  *  (  -  x''®^’) 

The  gradients  3PRC/3N  and  3PRC/3E  (the  vertical  gradient  is  assumed  to 
have  negligible  importance)  are  obtained  by  best  linear  fit  of  measurements 
from  the  monitor  and  master  stations.  In  general  the  stations  will  not  be 
placed  directly  north  and  east  of  the  master  reference  stations,  and  the 
pseudorange  corrections  can  be  computed  directly  as  a  linear  function  of  the 
PRCs  at  the  reference  stations  in  the  network: 

PRC^Jser  _  PRC''®^’  +  aj  PRC''®^2  PRC''®^^, 
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where  the  coefficients  are  obtained  by  solving  for  the  (horizontal) 
position  of  the  user,  ,  in  terms  of  the  positions  of  three  reference 

stations,  and  minimizing  the  error  e: 


X 


+  aj  X 


The  term  £  is  an  error  term  that  absorbs  altitude  errors  or,  in  the  case 
of  two  reference  stations,  lateral  error  perpendicular  to  the  line  between  the 
stations.  By  picking  an  imaginary  directly  north  or  east  of  the  master 

station,  one  can  obtain  the  north  and  east  gradients  described  above. 

The  gradients  are  broadcast  as  a  relatively  short  message  at  relatively 
low  rates  from  the  master  station.  The  impact  of  the  message  is  lessened 
because  the  spatial  decorrelation  errors  are  small  compared  to  the  total  error 
carried  by  the  high-rate  differential  correction  message.  Precision 
requirements  for  reporting  the  gradients  is  only  a  few  digits,  all  that  is 
required  to  maintain  ranging  accuracy  to  a  level  of  well  under  one  part  per 
million.  The  message  update  rate  depends  primarily  on  how  quickly  the 
gradients  change  over  time.  Note  that  the  gradients  change  much  more  slowly 
than  the  total  error,  which  is  dominated  by  the  clock  dither  error  function  of 
S/A.  The  highest  expected  contributor  to  gradient  dynamics  is  ionospheric 
activity,  which  probably  produces  less  than  1  ppm  variation  in  100  seconds 
under  operational  GPS  conditions.  Preliminary  study  indicates  a  100  second 
update  rate  is  more  than  adequate.  At  this  rate,  the  gradients  add  only  a  few 
bits  per  second  to  the  master  station  message,  and  so  are  not  included  in 
Table  3.4-1. 

The  gradients  are  caused  by  three  basic  effects.  Examining  these  effects 
in  detail  assists  in  determining  the  efficacy  and  form  of  a  multiple  reference 
station  network.  The  magnitudes,  decorrelation  distances  and  the  decor¬ 
relation  times  of  the  differential  errors  are  of  particular  importance  in 
designing  the  system, 

6.2  Simulation  Description 

The  simulation  considers  a  number  of  scenarios,  with  separation  distances 
from  0  to  120  kra  and  using  a  single  reference  station  or  a  network  of  two  or 
three  reference  stations. 

A  block  diagram  of  the  simulation  used  in  predicting  navigation 
performance  is  provided  in  Figure  6.2-1.  The  simulation  is  a  hybrid  .Honte 
Carlo/linear  error  covariance  program.  The  Monte  Carlo  is  provided  by 
choosing  random  sa  ;ellite  geometry  cases  (a  random  time  for  the  so-called  21  +  3 
constellation).  For  each  case,  a  linear  error  covariance  model  was  used  to 
establish  accuracy  performance.  Results  were  averaged  over  1000  random  cases 

For  each  satellite  geometry,  models  are  used  to  compute  S/A  and 
ionospheric  and  tropospheric  delay  error  statistics.  These  are  combined  iiilo 
differential  range  error  statistics,  and  a  least  squares  solution  a  Igor  i  t.h.rn  is 
used  to  predict  the  one  sigma  positioning  error.  If  monitor  stations  art;  used 
as  well  as  a  master  station,  the  error  statistics  of  the  gradients  are 
included.  When  operating  with  a  3  station  network,  the  relative  geometry 
depicted  in  Figure  6.2-2,  representative  of  operations  off  the  coast  of 
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Figure  6.2-1  Navigation  Error  Simulation  Overview 
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Massachusetts,  is  assumed.  The  master  station  is  at  Logan  Airport,  and  two 
monitor  stations  at  Gloucester  to  the  northeast  and  Plymouth  to  the  southeast 
A  two  station  mode  was  also  studied,  with  two  reference  stations  140  km  apart 
along  the  Texas  coast  (Figure  6.2-3).  In  the  Massachusetts  Bay  case,  user 
separation  distances  of  0  km  to  120  km  from  the  master  reference  station  at 
Logan  were  considered;  for  the  Texas  coast,  user  separation  distances  of  0  k.m 
Co  120  km  from  the  midpoint  of  the  line  between  the  reference  stations. 

6.2.1  Error  Modelling 

6.2. 1.1  Computed  Orbit  Error 

The  major  component  here  is  the  intentional  degradation  in  the  broadcast 
orbit  parameters  to  unauthorized  users,  termed  Denial  of  Accuracy  or  S/A.  S/A 
degradation  is  a  combination  of  two  processes,  epsilon  and  clock  dither.  Only 
epsilon  causes  perceptible  differential  range  error;  clock  dither  is  usually 
completely  removed  by  differential  techniques.  The  epsilon  process  comprises 
the  broadcast  of  slightly  incorrect  orbit  parameters,  resulting  in  range  error 
on  the  order  of  30  meters  under  normal  conditions.  Geometry  considerations 
produce  different  range  error  for  different  locations,  which  can  be  explained 
as  follows.  Given  the  baseline  B  between  reference  and  user,  the  s"‘  ' '' ite 

position  error  dS,  and  the  range  vector  to  the  satellite  of  R,  the 
differential  range  error  due  to  epsilon  is  roughly: 

(k  X  (B  X  R))^  dS  /  r5. 

The  maximum  differencial  range  error  is  roughly  (BdS)/20 , 000 ,  where  B  is 
baseline  distance  in  kilometers  and  dS  is  the  radial  satellite  position  error 
in  meters.  Uith  typical  values  for  the  standard  deviation  of  dS  around  30  m, 
Che  range  error  increases  with  baseline  length  at  a  ratio  of  0.15  cm/km,  or 
1.5  ppm. 

The  epsilon  component  of  differencial  range  error  changes  very  slowly, 
and  very  predictably,  over  time.  It  is  not  a  driver  for  the  update  rate  of 
gradient  messages. 

6. 2. 1.2  Ionospheric  Delay 

If  the  user  and  the  reference  station  both  have  an  ionospheric  error,  but 
it  is  the  same  at  the  two  stations,  it  will  drop  out  under  differencing. 

Under  normal  continental  United  States  (CONUS)  ionospheric  conditions,  one  can 
expect  the  ionospheric  effect  on  range  error  to  vary  at  the  rate  of  less  than 
1  mm/km.  Infrequently  this  will  rise  as  high  as  5  nun/km,  and  on  rare 
occasions  into  the  10  mm/kra  range.  Outside  CONUS,  3  mm/kjn  is  not  uncommon  and 
one  instance  of  50  ram/km  position  error  has  been  reported  in  Antarctica. 
Ionospheric  variations  under  Che  effects  of  the  auroras  in  the  polar  and 
subpolar  regions  can  be  quite  severe.  The  simulation  used  an  ionospheric  rate 
standard  deviation  of  1  mm/kra.  It  is  reasonable  to  assume  chat  large  scale 
ionospheric  effects,  such  as  moving  ionospheric  waves  and  diurnal  accic'itv, 
are  to  a  large  extent  linear  with  distance  on  the  100  km  scale.  Ocher 
ionospheric  activity,  such  as  scintillations  or  storms,  may  be  very  local 
Because  the  latter  effects  are  not  linear  over  a  wide  area,  they  will  not  l>e 
mitigated  by  networked  corrections,  and  are  modeled  in  rhe  simulation  as  a 
correlated  process  over  distance. 
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Ionospheric  delay  affects  the  user  only  if  he  is  relying  on  a 
differential  ionospheric  model  rather  than  direct  dual- frequency  measurement 
of  the  ionosphere.  The  model  becomes  important  if  only  one  carrier  frequency 
is  being  used  or  if  linear  combinations  of  both  frequencies  are  used  which 
preclude  use  for  ionospheric  correction  (e.g.  wide - lane ) . 

6. 2. 1.3  Tropospheric  Delay 

Like  ionospheric  delay,  most  tropospheric  error  will  drop  out  under 
differencing.  Compared  to  ionospheric  delay,  tropospheric  delay  tends  to 
decorrelate  over  shorter  distances.  The  wet  part  of  troposphere  effect  seems 
to  decorrelate  particularly  quickly.  As  a  general  rule,  the  more  accurate  the 
tropospheric  modelling,  the  shorter  the  modelling  error  correlation  distance. 
For  a  good  tropospheric  model  (e.g.,  the  Hopfield  model),  field  experience  has 
indicated  differential  tropospheric  error  on  the  order  of  2-3  cm  decorrelating 
over  tens  of  kilometers.  At  worst,  unusual  tropospheric  events,  such  as 
severe  electrical  storms,  can  cause  short-term  errors  on  the  order  of  many 
centimeters;  these  are  considered  rare  enough  to  not  be  included  in  the 
statistics.  Decorrelation  distances  are  short  enough  on  differential 
tropospheric  ranging  error  that  it  acts  similar  to  receiver  noise  at  longer 
distances . 

Table  6.1  shows  the  major  differential  range  error  sources  and  the  range 
of  error  standard  deviation  values  attributed  to  them,  as  well  as  the  values 
ised  in  the  performance  simulations. 


Table  6.1:  Differential  Range  Error  Contributors 


std  dev- range 

std  dev-sim  value 

tropsphere 

0.1  -  0.5  cm/km’/2 

0.3  cm/km^^^ 

inonsphere  (trend) 

0.0  •  1.0  cm/km 

0 . 1  cm/km 

ionosphere  (random) 

0.0  -  0.2  cm/km^^^ 

0.07  cm/km^^^ 

1  orbit  error 

0.0  -  0.2  cm/km 

0.15  cm/km  jj 

raultipath/noise 

0.5  -  3.0  cm 

1.0  cm 

6.2.2  Results  Summary 

The  simulation  results  are  summarized  in  Figures  6. 2 -A  through  6.2-6. 

Tho  results  show  that  decimeter- level  positioning  is  indeed  achievable, 
although  in  the  vertical  it  is  somewhat  dicey  at  longer  separation  distances. 
The  accuracy  figures  represent  an  average  over  1000  random  satellite 
geometries  in  Massachusetts  Bay  for  both  a  single  reference  receiver  and  a 
three  station  network,  and  both  dual  and  single  frequency  receivers. 
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Accuracies  are  presented  as  sCc*iidard  deviations  in  horizontal  (RSS  of  east  and 
north)  and  vertical  directions,  in  centimeters. 


As  expected,  dual  frequency  receivers  (used  in  ionosphere-corrected  mode) 
have  an  "error  floor"  roughly  three  times  higher  than  single  frequency  (LI) 
receivers  at  zero  separation.  This  is  because  receiver  noise  and  multipath  is 
misinterpreted  in  the  ionosphere  correction  process  as  ionospheric  delay, 
resulting  rn  a  range  error.  In  a  well-designed  system,  the  two  frequencies 
would  be  averaged  at  short  ranges  to  reduce  noise  and  multipath  effects,  with 
the  averaging  gradually  turning  into  ionospheric  correction  at  longer 
distances.  This  would  result  in  better  performance  than  a  single  frequency 
system  at  zero  baseline.  The  restalts  shown  below  are  strictly  ionospheric- 
corrected  dual  frequency  (ICDF)  and  do  not  reflect  such  a  system  design.  In 
fact,  the  results  indicate  chat  ionospheric  correction  is  a  hindrance  rather 
than  a  help  in  the  network.  This  is  expected  under  the  assumption  that 
ionosphere  effects  are  roughly  linear  with  distance  over  the  area  under 
consideration  and  thus  completely  observable  by  the  network. 

A  note  is  appropriate  about  the  general  performance  of  the  network. 

First,  performance  of  the  network  and  the  differential  (single  station)  are 
identici' 1  at  zero  separation,  because  the  gradients  are  not  used 
(equivalently,  the  weighting  coefficients  a^  and  a^  are  zero).  Second, 
network  performance  degrades  at  distances  much  larger  than  30  km  for  the 
network  in  this  study.  Up  to  this  distance,  the  user  lies  inside  the  triangle 
formed  by  the  three  reference  stations  and  the  PRCs  are  interpolated,  which 
tends  to  average  the  noise  effects.  The  optimal  point  is  at  the  centroid  of 
the  triangle  (at  roughly  20  km)  ,  where  all  weighting  coefficients  a,  -  aj 
aj  -  0.33  and  the  noise  is  diminished  by  a  factor  of  (3)^^^.  At  longer 
distances  the  PRCs  must  be  extrapolated,  the  weighting  coefficients  a^  become 
quite  large,  and  their  effect  on  the  noise  increases  proportionately  with 
distance.  This  especially  hurts  the  ICDF  system,  whose  higher  measurement 
noise  is  effectively  multiplied  by  the  RSS  of  the  weighting  coeff icien  :s .  The 
weighting  coefficients  for  the  different  distances  are  listed  in  Table  6.2. 

Table  6.2:  Massachusetts  Bay  Network  Weighting  Coefficients 


separation 

_ ^0 

a2 

RSS 

0  km 

1.00 

0.00 

0.00 

1.41 

10  km 

0.64 

0,17 

0.17 

1.21 

20  km 

0.29 

0.35 

0.35 

1.15 

40  km 

-0.40 

0.70 

0.70 

1.46 

60  km 

-1,10 

1.05 

1.05 

2.10 

80  km 

-1.80 

1.40 

1.40 

2  .  86 

LOO  km 

-2.51 

1.75 

1.75 

3,67 

120  tun 

-3.21 

2.10 

2  . 10 

4 . 49 

The  horizontal  accuracy  (Figure  6, 2 -4a)  indicates  that  a  single  ireq-unicv 
network  is  the  most  effective  in  providing  high-quality  positioning, 
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Dl/Terenlial  Horizontal  Accuracy:  Network  vs.  Sltiglc  Reference 


performing  better  than  a  single  reference  differential  or  ICDF  networks.  Tlte 
ICDF  network  degrades  badly  at  longer  distances  due  to  the  extra  noise  in  the 
carrier  measurements  due  to  the  ionospheric  correction  process  and  the  size  of 
the  weighting  coefficients  a-. 

The  vertical  accuracy  (Figure  6. 2 -4b)  is  roughly  three  times  worse  than 
the  horizontal  accuracy,  with  maximum  error  standard  deviations  on  the  order 
of  20-25  cm  at  120  km.  Note  again  the  relatively  good  performance  of  the 
single  frequency  network  and  also  the  ICDF  single  -  station  differential.  This 
indicates  that  ionospheric  error  is  perhaps  the  largest  contributor  to 
vertical  error  at  the  larger  separation  distances. 

Table  6.3  shows  a  breakdown  of  the  accuracy  results  both  for  the  two 
cases  illustrated  in  the  Figure  6.2-4  and  a  third  case,  the  two- reference 
station  scenario  off  the  Texas  coast,  illustrated  in  Figure  6.2-5.  In 
addition,  the  results  corresponding  to  Figure  6.2-7  are  summarized. 

To  test  the  performance  of  the  network  on  single  frequency  under 
worsening  ionospheric  conditions,  the  same  simulation  was  run  but  with  the 
ionospheric  trend  values  doubled  to  2  ppm  (Is).  These  values  are  typical  of 
bad  ionospheric  conditions  in  CONUS,  and  the  results  are  shown  in  Figure 
6.2-6.  The  higher  ionospheric  values  had  no  effect  on  the  network  solutions, 
as  the  network  corrections  effectively  estimate  and  remove  the  ionospheric 
trends.  Single  differential  stations  suffered  badly,  barely  able  to  meet  the 
vertical  decimeter- level  accuracy  at  10  km. 

The  last  study  examined  the  performance  of  wide -lane.  With  a  two 
frequency  receiver,  especially  in  codeless  mode,  it  is  quite  possible  that  the 
wide  lane  integer  will  be  resolved  but  the  LI  and  L2  lanes  will  remain 
ambiguous.  This  is  equivalent  to  uncertainty  in  the  differential  ionosphere. 
In  this  case,  navigation  will  be  done  using  the  wide -lane  phase,  which  is  the 
difference  between  LI  and  L2  phases.  The  sensitivity  of  this  to  troposphere 
and  orbit  error  are  identical  to  the  previous  cases,  as  these  are  independent 
of  frequency.  However,  the  effect  of  differential  ionosphere  is  increased  by 
a  factor  of  30%  over  single  frequency,  and  the  effect  of  multipath  and  noise 
is  increased  by  nearly  six-fold  over  single  frequency  and  doubled  over  ICDF. 
Because  of  these,  response  is  much  worse  than  either  single  frequency  or  ICDF. 
A  comparison  of  the  wide- lane  and  single  frequency  performances  is  shown  in 
Figure  6.2-7a  and  6.2-7b  for  the  Massachusetts  Bay  three  station  network. 

Although  the  results  for  wide-lane  are  not  by  themselves  impressive,  the 
accuracy  from  the  wide-lane  positioning  technique  can  conceivably  be  used  to 
acquire  the  LI  and  L2  lanes.  The  leap  from  wide-lane  to  the  base  frequencies 
is  equivalent  to  measuring,  or  knowing  a  priori,  the  differential  ionospheric 
effect  on  the  LI  and  L2  frequencies  fairly  accurately.  This  is  probably 
possible  over  a  few  kilometers,  and  probably  can  be  resolved  over  longer 
distances  with  a  network.  The  dependence  on  good  modeling  of  ionospheric 
behavior  over  these  distances  is  so  critical  that  we  feel  an  effective  study 
can  not  be  done  in  a  simulation  with  the  current  knowledge  of  ionospheric 
models . 
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Figure  6.2-4b  DUTerentlal  Vertical  Accuracy:  Network  vs.  Single  Reference 
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Separation  Distance,  User  to  Midpoint  between  References  (km) 


Figure  6.2-6  Single  Frequency  DlfTerentlal  Accuracy  Under  Severe  Ionosphere  Condition 
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Figure  6.2-7a  -  Comparison  of  Widelane  and  Single  Frequency  Accuracl 
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Figure  6.2-7b  -  Comparison  of  Wldelane  and  Single  Frequency  Amiraries 


6.3  Summary  and  Conclusions 


The  results  of  the  simulations  indicate  chat  the  performance  of  the 
differencial  system  is  sensitive  to  networking  in  the  environments  expected  to 
be  encountered  in  normal  dredging  operations.  A  decimeter- level  vertical 
accuracy  can  be  sustained  out  Co  20  km  with  single  frequency  differential  lud 
40  km  or  perhaps  further  with  a  network  of  three  reference  scacior\s 
surrounding  the  user.  Lastly,  if  dual  frequency  is  used,  the  system  must  not 
perform  a  simple  ionospheric  correction  on  the  carrier  data,  as  this  will 
hinder  performance  considerably. 

The  subject  of  wide-area  networks,  those  with  several  hundred  kilometers 
between  stations,  must  also  be  viewed  with  suspicion  in  the  light  of  the  above 
simulation  results.  Wide-area  networks  postulate  that  the  gradients  are 
relatively  constant  over  wide  areas.  Specifically,  these  gradients  could  be 
measured  and  broadcast  from  a  network  of  perhaps  ten  to  thirty  monitor 
stations  over  CONUS,  as  might  be  established  by  the  US  Coast  Guard  (USCG) , 
National  Oceanographic  and  Atmospheric  Administration  (NOAA)  or  USAGE.  This 
would  work  quite  well  for  dual  frequency  receivers,  which  can  measure 
ionosphere;  in  fact,  such  a  system  has  been  proposed  to  the  National 
Aeronautics  and  Space  Administration  (NASA)  for  shuttle  emergency  landing 
navigation.  For  single  frequency  and  wide-laning  situations,  however,  the 
wide-area  concept  is  considerably  higher  risk,  mainly  because  of  the  dominance 
of  ionospheric  effect.  The  wide-area  concept  depends  on  the  linearity  of  the 
differential  range  error  over  large  areas,  and  such  linearity  for  ionospheric 
behavior  has  not  been  tested  or  documented.  If,  as  expected,  wide-laning  is  a 
critical  technique  in  performing  the  dredging  surveys,  the  wide-area  network 
might  be  the  only  method  to  guarantee  decimeter  level  vertical  accuracy. 

Field  testing  on  large  scale  fiducial  networks  focusing  on  this  specific 
problem  is  required. 
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7.0  Cost  and  Risk  Analysis 
7.L  GPS  Receiver  Costs 

Table  7.1-1,  abstracted  from  [26],  lists  current  off-the-shelf,  .survey 
quality  GPS  receivers.  Summary  cost  data,  also  abstracted  from  [261,  appears 
as  Figure  7.1-1.  As  can  be  seen  from  the  figure,  costs  range  from  $35K  to 
$50K  for  high  quality  receivers.  Excluded  from  Che  cost  data  is  the  Rogue 
receiver,  designed  by  Jet  Propulsion  Laboratories  (JPL)  and  built  by  Allen 
Osborne,  which  sells  for  $80K.  A  lower  size,  weight  and  power  version  ot  the 
Rogue ,  named  Turbo -Rogue,  which  should  also  have  improved  cracking 
performance,  is  being  developed  and  is  expected  to  cost  in  the  neighborhood  of 
$40K,  which  should  be  competitive  with  other  survey  quality  receiver  designs. 
In  addition,  the  Mini-Rogue,  also  expected  to  sell  for  $40K.,  is  being 
developed  by  Allen  Osborne;  its  projected  availability  date  is  December  1990. 
Mote  that,  of  the  off-the-shelf  GPS  receiver  designs  listed  in  Table  7.1-1, 
only  the  Rogue  is  able  to  utilize  a  full  wide-lane  capability  without  use  of 
the  P-code.  Of  course,  all  designs  which  have  a  P-code  tracking  capability 
support  wide-laning,  e.g.,  the  TI4100  and  the  WM102 .  Half -wide - laning  is 
possible  with  MM2816,  MM2816AT,  and  Macrometer  II,  in  addition  to  the  LD-Xll, 
and  the  4000ST  and  4000SLD. 

Limited  design  efforts  are  underway  for  the  development  of  an  integrated 
GPS/GLONASS  receiver:  Magnavox  has  financed  an  International  Research  and 

Development  (IRAD)  project  for  about  2  years  [7],  Ashtech  is  beginning  a  joint 
development  with  the  former  USSR,  and  the  Federal  Aviation  Administration 
(FAA)  has  sponsored  a  development  at  Lincoln  Laboratory.  Due  to  the 
relatively  immature  state  of  these  developments,  no  reliable  cost  data  are 
available . 

7.2  Reference  Station  Costs 

As  indicated  in  Figure  3.1-1,  the  major  components  of  the  reference 
stations  are  the  GPS  Receiver,  a  processor  and  display  and  a  communications 
subsystem.  The  receiver  costs  have  already  been  addressed;  the  processor  and 
display  costs  are  expected  to  be  in  the  $3K-$5K  range,  assuming  a  10  mHZ  AT 
clone,  with  hard  disk,  and  a  monachrorae  graphics  capability.  The 
communications  subsystem  costs  vary  significantly  between  the  master  and 
monitor  reference  stations.  For  the  master  station: 

•  $4K  for  VHF/UHF  radio. 

•  $300  for  each  of  2  modems  (reception  of  monitor  station  data) . 

For  each  monitor  station,  a  modem  costing  $300  is  required. 

7.3  INS  Costs 

Reliable  cost  data  has  only  been  obtained  for  the  high  quality  INS. 

Costs  are  generally  expected  to  be  in  the  $80K  to  $100K  range;  specifically, 
the  MAPS  DRU  unit,  available  to  the  Array  for  approximately  $86K,  meets  or 
exceeds  the  navigation  accuracy  requirements  for  the  high  quality  INS  as 
assumed  in  this  study.  Since  the  MAPS  DRU  unit  was  not  designed  for  this 
particular  application,  however,  some  adaptation  (specifically  of  the 
interface  to  the  outside  world)  will  be  needed. 
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Table  7. 11  'OfT-The- Shelf  Code  GPS  Survey  Receivers  17-11 
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Year  of  Receiver  Inlroduclioti 


7.4  Barometric  Altimeter  Costs 


Two  barometric  altimeter  units  adequate  for  use  on  this  system  are  the 
following; 

•  A  digital  barometer,  built  by  Atmospheric  Instrumentation  Research, 

Inc.,  which  sells  for  $800. 

•  The  Rosetta  Microsystems  Falcon  Altitude  Encoder,  which  sells  for 

$100-$2C0. 

7  5  Rubidium  Frequency  Standard 

As  shown  in  Section  5,  the  use  of  a  rubidium  frequency  standard  can 
extend  maximum  tolerable  gap  lengths  significantly  when  operating  with  3 
satellites,  so  its  use  in  the  proposed  system  is  highly  recommended.  It  is 
expected  to  cost  in  the  neighborhood  of  $10K. 

7.6  Computer  Costs 

As  previously  referenced  in  Section  3,  the  major  software  functions  which 
execute  in  the  ship-based  computer  include  interface  software,  differential 
correction  calculation  software,  GPS/INS  integration  software,  and  ambiguity 
resolution  software.  The  major  loading  on  computer  resources  (i.e., 
throughput)  is  expected  to  be  the  GPS/INS  Kalman  filter  and  the  ambiguity 
resolution  algorithm,  which  run  at  1  Hz  and  process  the  available  (carrier) 
pseudorange  measurements.  Of  course,  the  ambiguity  resolution  software  will 
only  execute  when  integers  need  to  be  resolved,  so  it  will  execute  whenever 
lock  is  lost  on  a  satellite  in  view,  or  whenever  a  satellite  is  added  to  the 
active  tracking  set.  On  a  COMPAQ  386/25  computer  at  Trimble,  the  56  state 
covariance  simulation  runs  approximately  20  times  faster  than  real  time,  so 
the  17  state  Kalman  filter  can  be  expected  to  run  roughly  100  times  faster 
than  real  time.  The  ambiguity  resolution  software  can  be  made  to  run  roughly 
5  timess  faster  than  real  time  on  a  Macintosh  II  computer.  On  a  386/25 
computer  this  can  be  expected  to  run  twice  as  fast  as  on  the  Macintosh.  Thus, 
a  386/25  seems  adequate  to  execute  the  real  time  software.  Its  cost, 
including  a  monochrome  display  and  hard  disk,  is  conservatively  estimated  to 
be  $7.5K. 

7.7  Software  Development  Costs 

It  is  difficult  to  estimate  software  development  costs,  due  to  the 
uncertainty  associated  with  hardware  integration,  the  level  of  documentation 
required  and  the  extent  to  which  existing  code  can  be  utilized.  The  software 
development  costs  reported  herein  will  assume  chat  there  are  no  significant 
hardware  integration  problems  (e.g.,  the  need  to  build  a  custom  interface 
between  the  INS  and  the  computer) ,  a  minimum  level  of  documentation  and  no 
portability  of  code.  The  assumed  documentation  level  corresponds  to  an 
equation- level  description  of  the  algorithms,  with  heavily  commented  source 
code.  Given  these  assumptions,  the  following  estimates  are  supplied  for  each 
of  the  major  software  modules: 

•  Ambiguity  resolution  software-  18  mm. 

•  Kalman  filter  integration  software-  15  mm. 

•  Reference  station  software-  9  mm  (Master  Station) ,  6  mm  (Monitor 

Station) . 

•  Interface  software-  6  mm. 

•  Software  installation  and  test-  6  mm. 
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Thus,  a  total  of  60  mm  are  required.  Assuming  that  a  man  month  costs  $12K, 
the  total  software  development  costs  are  roughly  $720K. 

7.8  Estimated  Total  System  Costs 

The  total  hardware  costs,  excluding  the  Rubidium  frequency  standard  and 
the  barometric  altimeter,  and  operating  with  a  single  (master)  station  and  for 
all  hardware  elements  included  are  summarized  below. 

ESTIMATED  TOTAL  HARDWARE  COSTS 

(single  reference  station) 

$207K 

(all  hardware  elements  included) 

$294K 

ESTIMATED  TOTAL  SOFTWARE  DEVELOPMENT  COSTS 

$720K 


7.9  Identified  Risk  Areas 

Three  areas  of  potential  risk  have  been  identified  for  this  system: 

•  The  use  of  non  "off- the  shelf"  equipment. 

t  Reliance  upon  GLONASS  in  the  recommended  receiver  design. 

♦  Use  of  multiple  reference  stations. 

Recommended  use  of  the  turbo-Rogue  receiver  for  the  kinematic  phase 
positioning  system  represents  a  risk  factor,  since  it  is  not  yet  commercially 
available.  This  risk  can  be  reduced  significantly  if  the  mini-Rogue  is 
utilized,  which  should  be  available  soon  and  include  all  the  necessary 
features  of  the  turbo-Rogue.  Use  of  GLONASS  is  considered  a  risk  area  since 
available  information  is  not  complete;  e.g.,  can  the  signals  be  encrypted  or 
turned  off?  Multiple  reference  stations  are  required,  especially  for  longer 
distance  separations  (e.g.,  approaching  100  km).  This  is  an  undesirable 
feature  of  the  system,  since  there  is  a  maintenance  cost  associated  with  each 
reference  station,  and  there  can  be  a  potential  for  vandalism.  The  vandalism 
can  perhaps  be  reduced  or  eliminated  by  placing  the  reference  stations  at 
(secure)  USAGE  facility  locations.  In  addition,  the  networking  concept,  as 
analyzed  in  the  previous  section,  has  not  been  proven  or  demonstrated  using 
real  hardware.  Since  it  is  needed  to  maintain  the  required  accuracy,  it 
would  have  to  be  demonstrated  and  tested  as  part  of  the  kinematic  positioning 
system 
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8.0  Evaluation  of  Growth  Potential 

8.1  Alternative  System  Uses 

Several  other  possible  uses  of  a  kinematic  phase  positioning  system  are 
referenced  in  the  open  literature,  in  addition  to  the  possible  use  on  board  a 
helicopter  of  immediate  interest  to  USAGE.  Its  use  as  part  of  a  precision 
range  tracking  system  is  reported  in  [25],  This  reference  documents  flight 
test  results  performed  at  NASA/Wallops  Flight  Facility  using  TI4100  GPS 
receivers  and  a  laser  altimeter  for  a  vertical  reference.  Agreement  at  the 
decimeter  level  is  reported.  In  [21],  an  integrated  GPS/INS  system  is 
proposed  for  use  in  highway  surveying  applications.  A  Litton  LTN  90-100 
strapdown  INS  is  integrated  with  2  Trimble  4000SX  GPS  receivers;  an  accuracy 
of  10-15  cm  is  expected  from  the  tests.  In  [20],  flight  tests  are  performed 
using  a  kinematic  phase  positioning  system  on  board  an  Orion  P-3C  aircraft: 
accuracies  of  12  cm  are  reported  for  relatively  poor  GPS  geometries; 
analytical  projections  for  good  GPS  geometries  are  1-2  cm.. 

8.2  Problems  in  Helicopter  Integration 

Two  possible  problems  exist  in  the  use  of  the  system  on  board  a 
helicopter;  the  possible  increased  frequency  of  cycle  slips,  increased  INS 
error  excitation  caused  by  the  more  severe  dynamic  environment  and  satellite 
obscuration  induced  by  the  helicopter  blades.  The  increased  frequency  of 
cycle  slipping  can,  of  course,  be  reduced  by  increasing  the  bandwidth  of  the 
carrier  tracking  loop;  this  unfortunately  brings  an  inevitable  increase  in  the 
noisiness  of  the  derived  carrier  phase  measurements.  This  tradeoff  generally 
results  in  an  optimum  bandwidth  selection  for  the  expected  djrnamics  and 
acceptable  noise  level  in  the  measurements.  This  optimization  requires  a 
precise  specification  of  the  worst  case  dynamics  and  noise  environment  which 
may  be  experienced.  Given  this  data  for  the  helicopter  environment,  then,  the 
probability  of  cycle  slips  can  be  reduced  to  an  acceptably  low  level.  Th  ■* 
possible  obscuration  by  the  helicopter  blades  has  2  potential  solutions;  the 
GPS  antenna  can  simply  be  installed  above  the  blades  or  the  signal  processing 
software  within  the  receiver  can  be  designed  to  "look"  for  the  signal  in 
between  the  blade  crossings.  This  latter  approach  has  been  demonstrated  using 
Trimble  receivers. 
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9.0  Conclusions  and  Recommendations 


The  major  conclusions  of  this  study  are  summarized  in  the  decision 
matrix  appearing  as  Table  9.0-1.  Note  that  certain  design  options  have  been 
deleted  from  the  list  of  possible  system  configurations.  Only  a  high  quality 
INS  is  considered,  due  to  its  improved  gap  filling  performance  in  the  absence 
of  GPS,  as  demonstrated  in  the  simulations  of  Section  5.  The  possible  use  of 
a  single  frequency  C/A-code  receiver  and  a  C/A-code  receiver  which  derives  L2 
information  by  signal  squaring  is  not  given  serious  consideration,  since  their 
ability  to  resolve  integers  is  significantly  degraded  relative  to  the  other 
receiver  types  (as  demonstrated  in  the  simulations  of  Section  4).  Finally, 
only  rubidium  standard  clocks  are  considered,  based  upon  their  enhancement  of 
measurement  gap  filling  capabilities,  as  reported  in  Section  5  (i.e.,  no 
consideration  is  given  to  a  very  low  cost  quartz  clock) . 

The  acronyms  used  to  describe  the  GPS  receiver  types  are;  P,  for  P  code 
receiver,  and  C/A  XC ,  for  a  dual  frequency  C/A-code  receiver  which  derives  L2 
information  through  cross -correlation.  In  the  description  of  the  possible 
aiding  sensors,  BARO  denotes  barometric  altimeter.  Finally,  S  denotes  a 
single  differential  reference  station,  while  N  refers  to  a  differential 
network  comprised  of  3  stations.  The  navigation  performance  measures  are 
expressed  as:  one  sigma  horizontal  position  error  and  one  sigma  vertical 
position  error,  with  both  numbers  expressed  in  centimeters.  Note  that  the 
navigation  performance  appearing  in  the  Cable  corresponds  to  results  for  which 
the  LI  integers  have  been  resolved;  they  are  therefore  somewhat  optimistic 
projections,  since  LI  convergence  conditions  were  not  established  as  part  of 
this  study.  More  pessimistic  projections  can  be  derived  from  the  wide -lane 
performance  numbers  appearing  in  Table  6,3.  Separate  secs  of  numbers  appear 
for  0  and  100  km  separation  distances.  The  performance  of  the  ambiguity 
resolution  algorithm  is  expressed  in  terms  of  the  maximum  time  to  converge  to 
a  correct  set  of  wide- lane  integers  on  all  satellites  in  view  (expre.ssed  in 
minutes)  at  the  maximum  indicated  separation  distance  (in  km).  Constraints  on 
the  magnitude  of  the  modelled  multipath  error  (i.e.,  M-1)  necessary  to  achieve 
the  indicated  convergence  time  are  noted.  The  maximum  gap  length  numbers  (in 
seconds)  are  for  0  and  100  km  separation  distances,  respectively,  with 
complete  loss  of  GPS  during  the  gap. 

Based  upon  this  decision  matrix,  several  possible  systems  emerge  as 
strong  candidates  for  USAGE  applications.  The  clear  choice  for  a  GPS  receiver 
is  either  the  Turbo  Rogue,  under  development  at  JPL  or  the  Mini-Rogue,  soon  to 
be  made  available  by  Allen  Osborne;  they  represent  the  ideal  combination  of  a 
P-code  receiver  with  a  backup  C/A-code  tracking  capability  which  can  support 
wide-laning.  It  should  be  competitively  priced  at  $40K.  If  the  risk 
associated  with  its  availability  is  determined  to  be  too  great,  use  of  the 
original  Rogue  (built  by  Allen  Osborne)  should  yield  equivalent  performance, 
but  with  a  cost  penalty  of  an  additional  $40K  (and  a  weight  and  size  penalty). 
The  (optional)  use  of  a  receiver  which  tracks  GLONASS,  in  addition  to  GPS, 
increases  the  speed  at  which  integers  can  be  resolved  (by  roughly  40%).  Its 
use  by  USAGE  is  not  strongly  recommended,  however,  since  no  such  receiver  has 
been  built,  and  information  about  the  GLONASS  system  is  not  complete. 
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A  high  quality  INS  is  essential  for  maintenance  of  navigation  accuracy 
during  (complete)  GPS  outages.  Use  of  a  rubidium  clock  is  warranted  for 
fairly  long  periods  of  loss  of  a  single  satellite.  Use  of  the  barometric 
altimeter,  even  though  it  is  relatively  inexpensive,  is  not  warranted,  since 
it  will  not  enhance  performance  significantly,  and  a  zero  mean  sea  level 
damping  approach  should  yield  equivalent  performance.  However,  it  may  be 
required  as  a  metrological  sensor  for  performing  local  tropospheric 
corrections,  and  may  prove  to  be  useful  if  this  system  is  used  in  helicopter 
operations,  as  described  in  the  previous  section. 

The  use  of  multiple  reference  stations,  which  is  required  for  separation 
distances  in  excess  of  20  km  (when  working  with  LI  resolved  integers), 
represents  significant  risk  (note  that  system  performance  even  with  the 
network  does  not  meet  the  10  cm  requirement  at  100  km) .  No  differential 
network  has  ever  been  built,  so  would  have  to  be  undertaken  as  part  of  this 
system  development.  In  addition,  there  could  be  significant  costs  associated 
with  maintaining  each  reference  station  and  securing  them  from  vandalism. 
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APPENDIX  1;  GPS  EPHEMERIDES 
[4-3] 


Summary  of  nominal  constellation: 

a)  Semi  major  axis  at  orbit:  26,609.000  km  (same  for  all  satellites) 

b)  Eccentricity:  0  (same  for  all  satellites). 

c)  Inclination:  55°  (same  for  all  satellites). 

d)  Argument  of  perigee:  0  (same  for  all  satellites). 

The  following  variables  are  for  the  reference  time  26  November  1989,  00:00:00 

e)  Right  ascension  at  the  ascending  mode: 


325.730284° 

(SVl  - 

SV4) 

25.730284° 

(SV5  - 

SV8) 

85.730284° 

(SV9  - 

SV12) 

145.730284° 

(SVl  3 

-  SV16) 

205.730284° 

(SVl  7 

-  SV20) 

265.730284° 

(SV21 

-  SV24) 

f)  Mean  anomaly 


190.96 

SVl 

220.48 

SV2 

330.17 

SV3 

83.58 

SV4 

249.90 

S75 

352.12 

SV6 

25.25 

SV7 

124.10 

SV8 

286.20 

SV9 

48.94 

SVIO 

155.08 

SVll 

183.71 

SV12 

312.30 

SV13 

340.93 

SV14 

87.06 

SV15 

209.81 

SV16 

11.90 

SV17 

110.76 

SVI8 

143.88 

SV19 

246.11 

SV20 

52.42 

SV21 

165.83 

SV22 

275.52 

SV23 

305.04 

SV24 

7.914  .lO'’  rad 

sec 


g)  Rate  of  right  ascension; 
(same  for  all  satellites) 
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APPENDIX  1:  GPS  EPHEMERIDES  (CONTINUED) 


1  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-00 

O.OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-00  3.33288  lOOOOOOOOOD-00 
O.OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-00  5. 158391222000000D+03 
O.OOOOOOOOOOOOOOOD-00  9.599310890000000D-01  O.OOOOOOOOOOOOOOOD-00 
3.69947 1  OOOOOOOOOD-00-7.9 1405  lOOOOOOOOOD-09  O.OOOOOOOOOOOOOOOD-00 
O.OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-00 
O.OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-00 


2  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-00 

O.OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-OO  3.848102000000000D-00 
O.OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-OO  5.158391222000000D403 
O.OOOOOOOOOOOOOOOD-OO  9-5993 10890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
3.69947 lOOOOOOOOOD-00-7.9 1405  lOOOOOOOOOD-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


3  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.762554000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5. 158391222000000D+03 
O.OOOOOOOOOOOOOOOD-OO  9.599310890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
3.69947  lOOOOOOOOOD-00-7.91 405  lOOOOOOOOOD-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


4  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  1.458746000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.1 5839 1222000000D+03 
O.OOOOOOOOOOOOOOOD-OO  9.599310890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
3.69947 1  OOOOOOOOOD-OO-7.9 1405 1000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


5  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  4.361578000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.1583912220000000+03 
O.OOOOOOOOOOOOOOOD-OO  9. 599310890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
4.746669000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  0  OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 

6  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  6. 145653000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.158391222000000D+03 
O.OOOOOOOOOOOOOOOD-OO  9. 5993 10890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
4.746669000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
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7  573  0  0  0  0.00  000000.00  O.QOOOOOOOOOOOOOOD-00 

O.OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-00  0.440696000000000D-00 
O.OCKXKX)OOOOOOOOOD-O0  O.OOOOOOOOOOOOOOOD-OO  5. 158391222000000D+03 
O.OOOOOOOOOOOOOOOD-OO  9.5993 10890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
4.746669000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


8  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  2. 165954000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.158391222000000D+03 
O.OOOOOOOOOOOOOOOD-OO  9. 599310890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
4.746669000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


9  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  4.995 132000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.1 5839 1222000000D-^03 
O.OOOOOOOOOOOOOOOD-OO  9.5993 10890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
5.793867000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


10  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  0.854 164000000000D -00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5. 158391222000000D-t-03 
O.OOOOOOOOOOOOOOOD-OO  9.599310890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
5.793867000000000D-00-7.9 1 405 lOOOOOOOOOD-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


11  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  2.706657000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5. 15839 1222000000D -1-03 
O.OOOOOOOOOOOOOOOD-OO  9.5993 10890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
5.793867000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


12  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  3.206344000000000D-(X) 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5. 158391222000000D4(X3 
O.OOOOOOOOOOOOOOOD-OO  9.599310890000000D-01  O.0OOOO(XXX)0OOO0OD-00 
5.793867000000000D  OO-7.9 1 405 lOOOOOOOOOD-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


13  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOQ^)OD-O0  5.450663000000000D-00 
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O.CXKKXXXXX)000(X)OD-00  0.0<XXX)0(X)OO0O00OD-0O  5. 158391222(XXXXX3D+03 
O  OOOOOOOOOOOOOOOD-OO  9.5993 10890000000D-01  0.0CXXXXXXXXXXXX)0D-00 
0.557879(XXX)000<X)D-00-7.914051000000000D-C)9  O.CKXXX)OOOOOOOOOOD-(X) 
O.OOOOOOOOOOOOOOOD-CX)  0.(XXX)0(XXXXXXXXX)D-00  0.00(XXXXXXXXXX)00D-O0 
0.0(XXXX)OCXX)OOOOOD-O0  O.(X>C'X>(XXXXXXX>0OD-0O  0.O(X>0O0000000CKX)D-00 

14  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-00 

O.OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-00  5.95035  lOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-OO  5.158391222000000D+03 
O.OOOOOOOOOOOOOOOD-OO  9.5993 10890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
0.557879000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


15  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  1.519484000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.1 5839 1222000000D+03 
O.OOOOOOOOOOOOOOOD-OO  9.599310890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
0.557879000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
0,00(XXXXXXXX>0000D-00  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 

16  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  3.661875000000000D-00 
0  OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-OO  5.1 5839 1222000000D-I-03 
O.OOOOOOOOOOOOOOOD-OO  9.599310890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
0.557879000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
0  OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 

' »  0O(K)O0OOOOO0OOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 

17  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

0  OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  0.207694000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.1 5839 1222000000D+03 
O.OOOOOOOOOOOOOOOD-OO  9.599310890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
1.605076000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
U. OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 

18  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  1.933127000000000D-00 
OOOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.158391222000000D+03 
O.OOOOOOOOOOOOOOOD-OO  9. 599310890000000D-01  O  OOOOOOOOOOOOOOOD-OO 
1,6050760000000000-00-7.9140510000000000-09  0  OOOOCXXXXXXKXXXID-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOCXXXXXXXXXX)OOD-(X> 


19  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  2.511 180000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD  00  5. 15839 1222000000D+03 
O.OOOOOOOOOOOOOOOD-OO  9.599310890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
1  605076000000000D-00-7.9140510000000PC)f>-09  O.OOOOOOOOOOOOOOOD-OO 
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0.(XXXXXXXXX)OOOOOD-C)0  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-CX) 
0.000(XXXXX)00000(3D  00  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


20  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  4.295430000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.1583912220000000-^03 
O.OOOOOOOOOOOOOOOD-OO  9.599310890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
1.605076000000000D-00-7.914U51000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOGOOOOGOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


21  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOCOOD-00  0.914902000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.158391222000000D-t-03 
O.OOOOOOOOOOOOOOOD-OO  9.599310890000000D.01  O.OOOOOOOOOOOOOOOD-OO 
2.652274000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


22  573  0  0  0  0  00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  2.894280000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.1 5839 1222000000D-f03 
O.OOOOOOOOOOOOOOOD-OO  9.599310890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
2.652274000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


23  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  4.808731000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.158391222000000D+03 
O.OOOOOOOOOOOOOOOD-OO  9.599310890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
2.652274000000000D-00-7.91405100000000CD-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 


24  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.323952000000000D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  5.1583912220000000-^03 
O.OOOOOOOOOOOOOOOD-OO  9.5993 10890000000D-01  O.OOOOOOOOOOOOOOOD-OO 
2.652274000000000D-00-7.914051000000000D-09  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOCXXDD-OO  O.OOOOOOOOOOOOOOOD-OO 
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APPENDIX  2:  GLONASS  EPHEMERIDES 
[4-1] 


Summary  of  nominal  constellation: 

a)  Semi  major  axis  of  orbit:  25.506.783  km  (same  for  all  sateilitesi. 

b)  Eccentriticy :  0  (same  for  all  satellites). 

c)  Inclination:  64.8°  (same  for  all  satellites). 

d)  Argument  of  perigee:  0  (same  for  all  satellites). 

The  following  variables  are  for  the  reference  time  7  September  1987,  00:00:00 

e)  Right  ascension  at  Che  ascending  mode: 

73.000°  (SVl  -  SV8) 

193.000°  (SV9  -  SV16) 

313.000°  (SV17  -  SV24) 

f)  Mean  anomaly 

0.00°  SVl 
45.00°  SV2 
90.00°  SV3 
135.00°  SV4 
180,00°  SV5 
225.00°  SV6 
270.00°  SV7 
315.00°  SV8 
15.00°  SV9 
60.00°  SVIO 
105.00°  SVll 
150.00°  SV12 
195.00°  SV13 
240.00°  SV14 
285.00°  SV15 
330.00°  SV16 
30.00°  SV17 
75.00°  SV18 
120.00'  SV19 
165.00°  SV20 
210.00°  SV21 
255.00°  SV22 
300.00°  SV23 
345.00°  SV24 


g)  Rate  of  right  ascension: 
(same  for  all  satellites) 


-5.875  .10'’  r^ 

sec 


1  15 


APPKNDIX  2;  GLONASS  EPllEMEkM)P:S  (€0NTLNUP:D) 


25  573  U  0  0  0  CX)  (KXXXX>  00  0  0(XXXKKKK)0(XXXX>D  3K) 

0.00(XXXXXXXXXXXX)D - CX)  0.0(XXXXXXXXXXXXX)D ■  tX)  3.42 1  B92(XXXXXX>(X)I>  (X) 
0.000(XXXXXXXXXXX)DXX)  0.(XXXXXXXXXXXXXX)D-(X)  5,0rX)424u54(XXXKXX>  ttXl 
0  00(XXXXXKXXXXXX)D  (X)  1. 1309734('XXXXXXX)D-{X)  O.OOCXXXXXXXXXXXX)D  tXJ 
5.22963 lOOOOOOOOOD-00-5, 8749550000000000-09  0.0(XXXXXXXXXXXXX)D  (X) 
0.000000000000(XX)D-(X)0.000(XXXXXXXXXXXH)  (K)  O.OOCXXXXXXXXXXXXJD  (X) 
0,CXXXXXXXXXXXXXX)D-O0  0.0(XXXXXXXXXXXXXJl)  00  0.(XXXXXXKXXXXXXX)D  (X> 

26  573  0  0  0  0.00  000000.00  O.OOOOOOOOOCXXXXX)D  00 

0.0CXXXXXXXXXXXXX)D-^0O  0.0000(XXXXXXXXXX)D-(X)  4.20729 1(XXXXXXXX)D  (X) 
0.000(XXXXXXXXXXX)D-00  0.0000000000000(X)D-00  5.050-4240frUXXXXX)D*(X5 
O.OOOOOOOOOOOOOOOD-OO  1  130973400000000D  00  0.00000000(XXXXXX)D  00 
5.22963 1000000000D-00-5.874955000000000D-09  O.CXXXXXXXXXXXXXXiD  OO 
O.OOOOOOOOOOCXXXX)D-O0  0.0(XXXXXXXXXXX)v90r  O*-  i- .Xa.k^  /JOOOO(XXXX)D-CX) 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOf'r.^O  OO  O.OOOOOfXXXXXXKXX)!)  CX) 

27  573  0  0  0  O  CX)  OOOuOO.OOO.OO(XXX)0000(XXXX)D-OCj 

O.OOOCXXXXXXXXXXX)D-CX)  O.V  OOOOOfXWXXXXX)D-{X)  4.992669(XXXXXXX.X)D-(X) 
0.00(XXXXXXXXXXXX)D-(X)  O.CXXXXXXXXXXXy'  CCL-uu  5.050424054000CK)0D4CX5 
O.OOOOOOCXXX)OOOOOD-O0  1.130973400(XXXXX)D-00  O.OOOOOOCX)CXXXXXX)D-(X) 
5.22963 100CX)00(XX)D-CX)-5.874955000000CXX)D-09  O.OOOOOOCXX)0(XXXX)D-O0 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  o.oooooooo(X)oooooD-oo 
O.CXXXXXXXX>OOOOOOD-OOO.OOOOOOOOOOOOOCX)D-(X)  O.OOOOOfXXXXXX)OOOD-(X) 


28  573  0  0  0  O.OOOOCXXX),OOO.OOOOOOOOOOOOOOOD-(X) 
O.OOCXXXXXXXXXX)OOD-00  0.00000000000(XXX)D-00  5.7780870000CXXX)OD-CX) 
O.OOOOOOOOOOOOCXX)D-CX)  0.0o0000(XXXXXXXX)D-00  5.05042405-10{XXXX)D-4a3 
0.0000000(XXXX)0(X)D-(X)  1.13097340CXXXXXX}D-00  O.OOOOOOOOOCXXXX)OD  00 
5.22963 10CXXXXXXX)D-CX)-5.874955000000(XX)D-09  0.(XXX)00000(XXXXX)D-00 
0.00<XXXXXXXXXXX)OD-00  0.000000000(XXXXX)D-00  0.00CXXXXXXXXXXX)0D  fX) 
0.000(XXXXXXXX)CXX)D-CX)  O.OOOOOOOOOOOOOCX)D-O0  0  00CX)CXXXXXXXXXX)D-00 


29  573  0  0  0  O.OOOOOCX)OOOO.OOOOOOCXXXXKXXX)D-O0 

0.0000000(X)00(XXX)D-(X)  0.(XXX)OOOOCX)OOOOOD-O0  6.5634a.5C.KXXXXXX)0D-(X) 
0.000000(XXX)OOOOOD-00  O.OOOOOOOOOOOOOOOD-OO  5,050424054(XXXXX)D-403 
O.OOOOCXXXXXXXXXX)D-O0  1.13097340(X)OOOOOD-''X)  O.OOCKXXKXXXXXXXX)D  CX) 
5.22963 1000(XXXXX)D-00-5.87495500000(X)OOD-09  O.CXXXXXKXXXXXXXX)D-(X) 
O.OOCXXXXXXX-XXXXXDD-OO  O.0O000(X)OOOOOOOOD-OO  O.fXXXXXXXXXXXXXXJD  (X) 
0.fJ(XX)OCX)iXXXKX)CX)D-(X)  0.0(XXX)OOOOOOOOOOD-CX)  0  (XXX>  XXKXXXX-KXX)!)  CX) 

30  573  0  0  0  0  (X)  OOOOOO.CX)  O.OOOCX)0(XXXX)OCXX)D  CX) 

0.0CXXXXXjCXX.XXX)(X>D-CX3  0,0CXXX)0<XXX)CXX)00D-00  1.065698CXX.XXXKKX)D  00 
0.00CXXXXXXX,XXXXX)D-00  0.00CXX)000O(XXX)00D-CX)  5.050424054(X.XXXX)D*(X) 
0.0(XX)00(XXXXXXXX)D-CX)  1.13O9734CK)0(XXKX)D-0O  O.CXXXXXXXXXX)CXXX)D  00 
5.22963 )CXXXXX)0(X)C  CX)-5.8749550CX)(XX)000D-09  0.(XXXX)CXXXXXXXX.X?I)-CX) 
O.CKXXXXXXXXXXXXXjD-OO  0.0000000(X)000<XX)D-00  O.CXX)OOOOOOOCXyXX)D  00 
0.00(XXXX)0000(XXK)D-00  O.OOOOOOOOOOOOOOOD-OO  0.0OOOOCXXXXK30OO0D  00 


31  573  0  0  0  0<X>0<mmm0.(XX)(XHXXXHX)O0O0D-O0 

0  CXXXXKKXXXXXXXXID-CK)  0  (.KXXXXXXKKXXXKXID  (X)  1  .o5 ! O^KXXXXXXKKX)!)  (X> 
O  CXXXKK.XXXXKXXXXiD-tX’l  0  CXXKKXXMXXXKKKXIDXX)  5.05<:)-r2405-l(XXXXX_>I)  +  fV^, 
O.OCXXXXXXXXXXXKXID  CX)  1  -  13097340<XXXXXX)DXX)  O  CXXXXXXXXXKXXXXID  (X) 
5  2296:HC3CaaXXX30D-(X)-;>  H7-}955(XXXXXXXX)D-09  OXXXXXXXXXXXXXXXID  (X) 
0  0000<XXKXXXXXXX)I)-(X.)  0  (XKXXXXXXXXXXXKID  00  O.CXXXXXXXXKxXXXX)!:)  (X) 
0,0000(XXXXXXXXXX)D-CX)  0.(XXXXXXXXKXXXXX)D-(X)  (X(XXXXXXXXXXXXkX)!>  (H) 

32  573  0  0  0  0.00  0(XXKX),00  0.0(XXXXXXXXXXXXX)D  00 

O.OCXXXXXXXXKXXXXID  00  O  (XKXXXXXKXXXKXX)D-OO  2.(130'} 94O0(XXXXXX)I  )-(K) 
0,OOOOOCOOOOOOO(X)D-0<3  0.(XXXXXXXXXXXXXX3D-CX)  5.050424O54iXXXXXn3O)3 
0  0(XXXXXXXXXXXXK3D-00  1 . 1 30973400(XXXXK)D-(X>  O X/OCKXXXXXXl-'XXXX)!^ -(X) 
5.229G3 iOO(XXXXXK)D  (X)-5.87<}95500{XK:KXXX)D-09  0  CKXXXXXXXXXXXKX:)D-(X) 
O.tXXXXXXXXXXXXXXlD  OO  0.0000(XXXXXXXXKX)D-00  0,(XXXXXXXXXXXXXX3D-(X3 
0  0CKXXXXXXXXXXXX)D-00  O.OCXXXXXXXXXXXXXlDXX)  0.0000{XXXXXXXXXX)D-00 

33  573  0  0  0  0.CK3  (XXXXX).00  0.(XXXXXXXXXXXXXX)D  -a') 

0.00000(X.XXXXXXXX)D-00  0.0000000(XXXXXX>0D-{X)  3,OS369200(XX.XKXXJD-('X) 
0.00(XXXXXXXXXXXX)D-(X)  0.00(XXXXXXXXXXXX)D-00  5,05042405'}CXXXXX)D  +  0:5 
0.0000{XXXXXXXXXX)D-00  1.1309734CXXXXXXX3D-(X)  O.OOCXXXXXXXXXXXXID-CXJ 
1.04084 1000(XXXXX)D-00-5.874955(XXXXXXXX)D-09  0.00(XXXXXXX>XKXX)D-(XJ 
O.CXXXKXXXXKXKXXXID-OO  O.CXKXXXXKXXXXXXXDD-CX)  O.CXXXXXXXXXXXXXXDD-OO 
O.OOOOCXXXXXXXXX)OD-O0  O.OOOCXXXXXXXXXXXID-OO  O.CKXXXXXXXXXXXXX)D-00 

3-4  573  0  0  0  0.00  0CKXXX).00  O.CXXXXXXXXXXXXXKID-OO 

0.000000(XXXKXXXX)D-00  0.0(XXXXXXXXXXXXX)D'00  4.469090000000000D-(X) 
O.OOOOOOOOOOOOOOOD-00  O.OOOOOOCXXKXIOOOOD-OO  5.050424054000000D+a3 
0.00(XKXXXXXXXXXX)D-(X)  1. 1309734(XKXX)OO0D-00  0.000000000(XXXKX)D-(X) 
1. 04084 10000000(X)D-(X)-5.87495500(XX3(XXX)D-09  0.0OOO(X)OOO00OOOOD-0O 
O.OOOOCKXXXXXXXXXID-OO  O.CXXXXXKXXXXXXXIOD-CX)  0,0{XXXXXX)OOOOOOOD-(X) 
0.000000000(XXXXK)D-{X)  0  (XXXXXXXXXXXKXXID-OO  0.000(XXXXXX>XXXX)D-00 

35  573  0  0  0  0.00  0(XXXX).00  0.0(X)00(XXXXXXXXK)D-00 

0.000(XXXXXXXXXXX)D-(X)  0.0(XXXXXXXKXXXXX)D-(X)  5.25‘}488000(XXXXX)D-00 
0.(XXXKXXXXXXXXXX)D-(X)  0.0(XXXXXXXXXXXXX)D-00  5.050424054CXXXXX)D4-03 
0.0000(XXXXXXXXXX1D-00  1.13097340(XXXXXX)D-00  O.OOCXXXXXXXXXX)OOD-O0 
1.04084100000(XXX)D-(X)-5.874955000(XXXX)OD-09  0.000000000000<XX)D-00 
O.OCKXXXXXXXXXXXWD-CX)  O.CXXXXXXXXXXXXXXID-OO  O.CKXXXXXXXXXXXXX)D-(X) 
0.(XKXKXXXXXXK3(XK)D-00  0.00(XXXX)OOOOOOOOD-CX)  0.(X)(XXXXXXXXXXXX3D-(X) 

36  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.CXXXKXXXXXXXXXXID  00  0.00000000(XXXXXX)D  00  6.039886O0(XXXXXX)D  CX> 
0.00(X)00(XXXXXXXX)D-(X>  0.0CXKX)0OOO(X)OOOOD-OO  5.05O424O54(XXXXX)D-t.n3 
0.(XXXXX)(XXXX)(XXX)D-00  1.1309734(XXXXXKX)D-00  O.rXXKXXXXXXXXXXXlD-CX) 

1 .04084  l(XXXXXXXX)D-00-5  8749550(XXXXXXX)D-09  O.OOOOCXXXXXXXXXXlD-fX) 
0.(XXXX>0(XXXXKXXX)D-(X3  0.(XXKXX)00000(XKX)D-00  0  OOOCXXXXXXXXXXXJD -00 
O.OCKXXXXXXKXXXXKID-OO  0.0000000000000(X)D-^(X)  0.00(XXXXXXXXXXXX)D  (X) 


37  573  0  0  0  0.00  (XXXX)0.00  0.00(XXXXXXK)OCXXX)D-(X) 

0.0000000(XKXKXXK)D-00  0.(XXXXXX)OOOOOOCK)D-(X)  0  542099000000000D-CK) 
O.OOOOOOOOOOOOOOOD-00  O.OOOOOOOOOOOOOOOD-00  5.050424054000000D+a3 


V  / 


O.OOOOOOOOOOOOOOOD-OO  1.1309734000000000-00  0.000(XXXXXXX)0(X)OD-00 
1 .04084 1000000000D-00-5.874955000000000D-09  O  CKKXXXXXXKXXXXJOD-OO 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  0.0000000000000{X)D-00 
O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  0.00000000(XXXXXX3D-00 

38  573  0  0  0  0.00  000000.00  O.OOOOOOOOOOOOOOOD-OO 

O.OOOOOOOOOOOOOOOD-OO  O.OOOOOOOOOOOOOOOD-OO  1.32749700000(XXX)D-00 
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